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ABSTRACT 



Condensation of R-113 was studied using an evaporator/condenser test platform. 
The condenser section contained four horizontal tubes (nominal outer diameter 15.9 
mm) forming a vertical in-line column with a pitch-to-diameter ratio of 2.25 and a con- 
densing length of 1.2 m. The condenser tubes could be operated either individually (ie. 
as a single tube apparatus) or as a small tube bundle. This allowed investigation of the 
effects of condensate inundation on different types of condenser tubes. Tubes tested were 
smooth copper tubes, copper/nickel KORODENSE roped tubes, KORODENSE tubes 
wrapped with wire, and copper/nickel finned tubes (26 fpi). The outside heat transfer 
coefficient, /i„, was calculated by experimentally determining the overall heat transfer 
coefficient, U, , and then using a modified Wilson Plot procedure. Great care was taken 
to ensure the results were not vitiated by the presence of noncondensibles. 

Results obtained with the smooth copper tubes are in agreement with published data 
and verify satisfactory' operation of the test platform. Furthermore, problems associated 
with the apparatus encountered by previous workers were successfully overcome. In 
comparison to the top smooth tubes, the copper/nickel top KORODENSE tube yielded 
about a 22% increase in /i„. When different diameters of wire were wrapped around the 
copper/nickel KORODENSE tubes, an optimum pitch-to-\\ire diameter of 7 was found 
yielding almost a 90% increase in h„ compared to the top smooth tube. This marked 
increase is presumably due to surface tension effects thinning the condensate film. The 
copper/nickel finned tubes gave the best enhancement (approximately a 7 fold increase). 
Comparison of the data obtained from the top copper/nickel finned tubes agreed well 
with the model of Beatty and Katz. 

The effect of condensate inundation was to reduce h„ for the lower tubes compared 
to the top tube in the bundle. Comparison of h„ for the second tube compared to the top 
tube showed that the effects of condensate inundation are reduced most by wrapping the 
KORODENSE tubes with fine wire. 
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NOMENCLATURE 



Standard Symbols 
A area, m? 

b fin spacing, m 

Cg parameter defined in Eq (2-11) 

C, specific heat, J/kg K 

D diameter, m 

e fin height, m 

g gravitational acceleration, 

h heat transfer coefficient, I'F/m^ A" 

heat transfer coefficient in flooded portion of finned tube, Wjm^ K 
heat transfer coefficient in unflooded region of finned tube, Wjm^ K 
hf heat transfer coefficient for fins in unflooded region, Wlni^ K 

k thermal conductivity, Wjm K 

L length, m 

L average condensing length defined by Eq (2-4) 

m mass flow rate, kg/s 

N number of tubes in a vertical column 

Nu Nusselt number 

Pr Prandtl number 

q heat transfer rate, W 

q" heat flux, WlnO 

R radius, m 

kf fouling resistance, K/W 

R„ wall resistance, K/W 

Re Reynolds number 

RCf film Reynolds number 

T temperature, K 

t fin thickness, m 

U overall heat transfer coefficient, fK/rn^ K 

xii 



V 



velocity, m/s 



Greek Symbols 

6 tape thickness, m 

t] fm efficiency 

^ aspect ratio of condensate surface, Eq (2-13) 

rotation angle normal to fm surface 

n dynamic viscosity, kg/m s 

p density, kglm^ 

a surface tension, N’/m 

4>f flooding angle defined in Eq (2-7) 

Subscripts 

c coolant 

eq equivalent 

ef effective 

f film 

fs fin sides 

ft fin tips 

i inside or inlet 

0 outside or outlet 

r root 

sat saturation 
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I. INTRODUCTION 



A. BACKGROUND 

In the last few years, budgetary considerations have led to a reduction in the funding 
level available for new vessel construction for the US Navy'. In construction of ships, 
weight can often be equated to cost. This has led the Navy to look at various ways by 
which the weight of auxiliary systems can be reduced. 

Of particular interest to the Nav'y has been the issue of reduced weight and therefore 
cost of refrigeration systems. Two major components of refrigeration systems are the 
evaporator and condenser. If the weight of these two components can be reduced, then 
the overall weight and presumably cost, can also be reduced. 

Condenser weight can be significantly reduced by the proper selection of materials 
for condenser construction. A major weight savings comes from changing the material 
which is used for construction of the condenser tubes. As an example, the DDG-51 was 
originally intended to have a refrigeration condenser made of titanium finned tubes. This 
would provide a substantial weight savings and a reduced cost (presumably). However, 
the system was not installed on the current flight of ships. Rather, traditional 
copper/nickel finned tubes were utilized for the refrigeration condenser. i 

In re-designing refrigeration condensers, changes in condenser tube material based 
solely on weight consideration should not be the only consideration. The reason for this 
is that a lightweight material may not have heat transfer characteristics which are com- 
parable to copper or even copper/nickel. Therefore, condenser tube designs have to be 
developed to maintain or even enhance heat transfer or else the condenser will have to 
be enlarged. This enlargement is necessary' in order to provide the requisite surface area 
to maintain the required level of heat transfer. 

Copper/nickel finned condenser tubes have been used in refrigeration condensers 
onboard US Navy ships for many years. However, the use of titanium, as was planned 
for the DDG-51 class ships, holds the promise of providing substantial weight savings 
(not to mention improved corrosion resistance on the tube side). Unfortunately, 
titanium has a thermal conductivity which is substantially less than that of copper or 
copper/nickel. The question which must therefore be addressed is whether titanium 

1 According to Mr. John Goodhue, NAVSEA code SEA 56Y15, the reason the units were 
not installed was that they could not be ready to support proposed construction dates lor fliglit one. 
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condenser tubes can be constructed in such a manner as to maintain heat transfer (due 
to enhancement) despite the lower value of thermal conductivity. 

A major problem along these lines is how to measure the performance of various 
types of condenser tubes. Typically, this is done with single tubes or small tube bundles. 
These tubes tend to be relatively short. Often, the experimental apparatus bares little 
resemblance to a shipboard condenser. Thus, any results become suspect when they are 
extended to the real situation. 

In order to examine heat transfer characteristics of various designs of condenser 
tubes, the best means is to construct a test platform composed of an evaporator and 
condenser with a size comparable to that of a shipboard system. Unfortunately, this is 
seldom feasible due to costs involved and the size of the plants. However, it is possible 
to construct a reasonably large scale test platform in which tube performance can be 
evaluated. This evaluation must include both single tube performance as well as tube 
bundle operation. The reason for this is that condensate forming on one tube falls onto 
lower tubes. This effect, known as condensate inundation, can result in decreased heat 
transfer on lower tubes due to a thickening of the condensate layer. Hence, another 
aspect which must be considered is whether condenser tube geometry can be changed in 
such a manner as to reduce the adverse effects of condensate inundation. This could 
help make up for decreased thermal conductivity associated with lightweight materials 
such as titanium. 

The construction of a suitable condenser/evaporator test platform was begun at the 
Naval Postgraduate School by Zebrowski (Ref 1]. The construction and initial operation 
of the test platform was completed by Mabrey [Ref 2] in 1988. Mabrey conducted ex- 
periments using R-113 and R-114 on smooth copper tubes as well as copper/nickel 
finned tubes (26 fpi). During the course of his experiments, an orange-peel or dimpled 
appearance was observed on the smooth copper tubes. Additionally, the copper/nickel 
finned tubes showed a buildup of "material" between the fins. This material was believed 
to be ice although this was not proved. The cause of the orange-peel appearance was 
also never identified. 

B. OBJECTIVES. 

This thesis concentrates on the condenser section of the test platform with specific 
interest in the condensation of R-113. R-113 was chosen for this study because of 

convienence , ie., it is a liquid at atmospheric pressure. The condenser section consists 
of four horizontal tubes mounted in a vertical in-line column. The nominal diameter of 
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the tubes is 16 mm with a condensing length of 1.2 m. While the condensing medium 
used in this study was R-113, the condensation of many different types of fluids can be 
evaluated. In addition, the condenser can be used to study heat transfer on almost any 
type of condenser tube subject to the nominal outer diameter and condensing length re- 
strictions indicated. 

The main objective of this thesis was to resolve the contamination problem found 
by Mabrey [Ref. 2], and to make the test platform operational. It was felt that such a 
system, once operational, could be used to study heat transfer using condenser tubes 
made of various materials and geometries. In addition, since the tubes are arranged in 
a vertical in-line column, the effect of condensate inundation could also be studied. 

C. SPECIFIC OBJECTIVES 

The specific objectives of this study were as follows: 

1. Resolve the contamination problem reported by Mabrey [Ref 2] 

2. Establish baseline condensation data by examining the condensation of R-113 on 
sm ooth copper tubes. 

3. Examine the heat transfer characteristics of R-113 on various tube geometries 

4. Investigate the effects of condensate inundation 
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II. LITERATURE SURVEY 



A. INTRODUCTION 

Condensation occurs when vapor molecules come into contact with a surface which 
has a temperature below the vapor saturation temperature. This vapor, coming into 
contact with the cooler surface, forms a liquid which can accumulate in two ways. If the 
liquid wets the entire surface, a continuous film of liquid is formed. On the other hand, 
if the liquid fails to wet the entire surface, droplets are formed. These two types of 
condensation are referred to as filmwise and dropwise condensation respectively. 
Dropwisc condensation yields larger heat transfer coefiicicnts but is difficult to maintain 
(Ref. 3]. All condensers today are therefore designed to operate with film condensation 
conditions. 

In 1916, Nusselt [Ref. 4] reasoned that as vapor condenses on a cooler surface, such 
as a horizontal tube, a continuous film of condensate will form. This film grows in 
thickness as the condensate flows around the tube under the influence of gravity. Hence, 
the thickness of the condensate is not uniform around a horizontal tube. Since the 
condensate acts as a resistance to heat transfer, Nusselt argued that the local heat 
transfer coefficient will be a maximum at the top of the tube where the film thickness is 
least and decreases around the tube as the condensate thickness increases. 

B. CONDENSATION OF REFRIGERANTS ON SMOOTH TUBES 

Goto et al. [Ref 5] examined film condensation of various refrigerants, including 
R-113, on smooth copper tubes. Their results showed values for the vapor-side heat 
transfer coefficient on the order of 1.2 to 1.5 kWjm^ K for a wall-film temperature dif- 
ference of 3 to 25 K. Furthermore, their results were about 10% below that predicted 
by Nusselt theory. 

. Masuda and Rose [Ref 6,7] reported that condensation of R-113 on smooth copper 
tubes yielded an outside heat transfer coefficient of about 1500 K for wall-film 
temperature differences of 12 to 20 K. Their results showed good agreement with 
Nusselt theor)’. 

Marto et al. [Ref 8] showed that for condensation of R-113 on smooth copper 
tubes, the vapor-side heat transfer coefficient was on the order of 1000 K for a 
temperature drop across the film of 20 to 25 K. Their values were also in good agreement 
with Nusselt theory'. 
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C. ENHANCEMENT OF FILM CONDENSATION 

The most commonly used procedure for enhancement of film condensation is 
througli the use of integral-fin tubes. With this type of condenser tube, two distinct re- 
gions exist; namely, a flooded and an unflooded region. Typically, as fluid condenses 
on a finned surface, it is pulled from the fin tip to the fin root by surface tension forces. 
Gravity then causes the fluid between fins to drain to the bottom of the tube where it 
drains off. Hence, condensation on a finned tube involves a combination of surface 
tension and gravity forces. It is the added surface tension which tends to enhance heat 
transfer over smooth tubes by thinning the condensate film at the fin tip. 

Since the present study uses R-113 as the working fluid, the remainder of this 
chapter will focus on previous experiments which have used low surface tension fluids 
such as refrigerants. For additional information involving other fluids, the reader is re- 
ferred to recent review articles by Marto [Ref 9,10], Webb [Ref 11] and Sukhatme [Ref 
12 ]. 

1. Condensation on Single Enhanced Tubes: Theoretical Models 

The first theoretical model of condensation on finned tubes was published by 
Beatty and Katz [Ref 13]. Their theoretical model is given by: 

4 Au 

hBK = nfy-J- + K-T ( 2 - 1 ) 



Here rif is the fin efficiency, ^,is the fin surface area, is the unfinned surface area, and 
A,/ is the effective total surface area. The above equation can be simplified to: 



^BK — 0.689 









eq 



( 2 - 2 ) 



where the equivalent diameter is given by: 



1 

^eq 




1.30)/^ 



± 

^ef 






(2-3) 



In this equation, 

- <Dl - D^) 



(2-4) 
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and 



Ae/=fl/^f+^u (2-5) 

The equation for the heat transfer coefTicient is a Nusselt-type condensation equation 
based on the equivalent diameter for a finned tube. This equation predicted the exper- 
imental results obtained by Beatty and Katz to within ± 20%. 

There are several important limitations to the Beatty and Katz model [Ref. 13]. 
It does not distinguish between fin tips and fin flanks, essentially assuming that the fins 
have no thickness. Therefore the model is not suitable for trapezoidal fins. A second, and 
perhaps more significant limitation is that the Beatty and Katz model fails to take into 
account surface tension effects. Instead, the model assumes that the flow of the 
condensate film is simply a function of gravity. While this may not pose a serious limi- 
tation when low surface tension fluids such as refrigerants are examined, this assumption 
can become critical when condensation of high surface tension fluids such as steam are 
investigated. 

Another limitation of the Beatty and Katz model is that it does not take into 
account condensate flooding. Rudy and Webb [Ref 14] argued that the flooded portion 
of the tube was ineffective in condensation. They therefore proposed that the Beatty and 
Katz model be modified by: 
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where R, is the outside radius, \{/ is the fin tip half angle and a is the fluid surface tension. 
It is interesting to note that this modification resulted in an underprediction of Rudy and 
Webb's R-11 condensation data by about 35 /o. 

Webb et al. [Ref 15] next argued that surface tension drainage must be consid- 
ered in the model. They derived the equation: 
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where c is the fin height, t is the fin thickness and b is the width of the space between 
the fins. Substitution into the Beatty and Katz expression yielded agreement with ex- 
perimental data of + 10%. 

Owen et al. [Ref 16] argued that the Rudy and Webb [Ref 15] model had to 
be corrected in order to account for heat transfer in the flooded part of the tube. Owen 
et al. proposed the equation: 



K = + ( 2 - 9 ) 

where is the average outside heat transfer coefficient and is the heat transfer coef- 
ficient in the flooded region of the tube. This value was estimated using a one- 
dimensional radial conduction model. This corrected model of Owen et al. resulted in 
agreement of ± 30% with earlier experimental data. 

In 1985, Webb et al. [Ref 17] proposed a model that included heat transfer in 
the flooded region as well as surface tension drainage on the fin flanks and gravity 
drainage on the interfin surface. Their model equation for the average heat transfer co- 
efficient is: 
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In this expression, is the heat transfer coefficient for the unflooded root surface be- 
tween the fins and hf is the heat transfer coefficient on the fin flanks in the unflooded 
area, and is the heat transfer coefficient in the condensate flooded region. The term 
Q is the fraction of finned tube surface flooded by condensate, or, 

4>f 

Cb=^--W- (2-11) 

where is the flooding angle given in equation (2-7). 

In the Webb et al. [Ref 17] expression for the average heat transfer coefficient, 
/i*, hj and are defined mathematically as follows: 

(2-12) 

where Rcj is the condensate film Reynolds number. 
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where S„ is the length of the convex surface and 0„ is the rotation angle normal to the 
fin surface [Ref 18]. The term (f is a parameter that characterizes the aspect ratio of the 
condensate surface. 

Finally, is defined as; 




This complex model of Webb et al. was able to predict heat transfer coefficients to 
within ± 20% of experimental data obtained with R- 1 1 . 

Honda and Nozu [Ref 19] and Honda et al. [Ref 20] have proposed the most 
comprehensive model to date. This model treats the heat transfer problem as one in- 
volving vapor-to-coolant heat transfer through a finned tube wall. In this model, 
trapezoidal fins were divided into flooded (f) and unflooded (u) regions. Equation (2-7) 
w'as used to determine the flooding angle. Condensate flow in between the fins was as- 
sumed to be driven by gravity whereas for the fin surface, condensate flow was assumed 
to be driven by both surface tension forces and gravity. The condensate film thickness 
on the fin surface was arbitrarily divided into a thick region, where no heat transfer oc- 
curred, and a thin film region. The equation of the condensate film thickness was solved 
numerically. 

Honda and Nozu [Ref 19] showed that the average Nusselt number could be 
written in terms of relevant parameters for the flooded (f) and unflooded (u) regions as: 



Here, and are the dimensionless average wall temperatures at the fin root and 4>f 
is defined as: 
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It should be pointed out that and T^/are determined by solving a circumferential wall 
heat conduction equation. This wall temperature variation must be considered since heat 
transfer rates through flooded and unflooded regions are very' different. 

Honda and Nozu [Ref 19] used their model to predict average heat transfer 
coefficients. These predictions were compared to experimental data on eleven different 
fluids and twenty-two tubes. The predicted values matched most experimental values to 
within ± 20%. 

2. Condensation on Single Enhanced Tubes: Experimental Results 

The first experiments conducted on enhanced tubes were performed by Katz and 
colleagues [Ref 21,22, 13 ]. These authors looked at condensation of R-12 and R-22 on 
smooth and finned tubes. Their results showed a significant enhancement in heat transfer 
for finned tubes. 

In 1971, Karkhu and Borovkov [Ref 23] examined condensation of R-113 on 
five finned tubes. Their results showed that the rate of condensation for these tubes with 
trapezoidal fins could be increased by 50 to 100% over smooth tubes. 

Arai et al. [Ref 24] and Kisaragi et al. [Ref 25] studied the condensation of R-12 
and R-113 respectively on finned tubes with spiral grooves cut along the length of the 
tube. Arai et al. showed that grooves with a pitch of 0.7 mm gave the best enhancement 
in heat transfer. Kisaragi and co-w'orkers also compared grooved finned tubes with tubes 
on which the fins had been covered with a porous coating. In comparison to the finned 
tube, a two fold enhancement was seen v/ith the grooved finned tube while a four fold 
enhancement was obtained with the finned porous coated tube. 

Carnavos [Ref 26] performed one of the first detailed studies of condensation 
of R-11 on enhanced copper tubes. In this study, eleven different tubes were tested. The 
tubes included finned tubes as wcU as finned tubes which had an axial spiral groove cut 
along their length. Various sizes and shapes of fins were included. Enhancements of 4 
to 6 were seen when compared to smooth tubes. 

Rudy [Ref 27] and Webb et al. [Ref 15] reported experimental data obtained 
for R-11 condensing on copper tubes with different fin geometries., Their results showed 
that increasing the fin density from 748 to 1378 fpm resulted in about a 70% increase in 
heat transfer presumably due to an increased surface area. 

Honda et al. [Ref 28] studied condensation of R-113 on copper finned tubes of 
various fin density. Their results showed a 6 to 9 fold enhancement for R-113 
condensation for finned tubes when compared to smooth tubes. A similar study by 
Kabov [Ref. 29] examined the condensation of R-12 and R-22 on tubes with different fin 
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densities His results also showed that the degree of enhancement is determined by fm 
spacing. 

Masuda and Rose (Ref. 6] examined heat transfer for the condensation of R-1 13 
on low integral-fm copper tubes. Fourteen different fin spacings (0.25 to 20 mm) were 
examined. The optimal fm spacing was 0.5 mm which resulted in a 7.3 fold enhancement 
in the vapor-side heat transfer coefficient. 

Lin and Berghmans [Ref 30] examined condensation of R-1 1, R-1 13, and R-1 14 
on an enhanced tube. This enhanced tube, known commercially as "EVERFIN", is made 
of trapezoidal fins (1378 fpm) with a spacing between the fins of 0.064 mm.2 This tube 
gave an 8 to 10 fold increase in the vapor-side heat transfer coefficient when compared 
to a smooth tube. 

Condensation of R-1 1 on a variety of enhanced copper tubes has been reported 
recently by Sukhatme et al. [Ref 31]. All tubes had trapezoidal shaped integral fms. The 
fm density, fm height and fm-tip half angle were systematically varied in order to deter- 
mine the optimal values which would result in maximum enhancement. The maximum 
enhancement obtained by these authors was about 10.3 relative to smooth tubes. This 
enhancement corresponded to a fm spacing of 0.35 mm, a fm height of 1.22 mm and a 
fm-tip half angle of 10 degrees. In addition, a set of specially enhanced tubes were tested 
by these authors. These tubes had axial grooves cut into them w'ith the trapezoidal fms 
present. An enhancement of 12.3 was obtained with these specially enhanced tubes 
when compared to smooth tubes. 

Marto et al. [Ref 8] studied condensation of R-1 13 on a variety of copper 
integral-fm tubes with rectangularly shaped fins. The optimum fm spacing was found 
to lie between 0.25 and 0.5 mm. The vapor-side heat transfer coefficient was enhanced 
on the order of 4 to 7 fold. 

Michael et al. [Ref 32] continued the studies of Marto et al. [Ref 8] by looking 
at the effect of root diameter on the condensation of R-1 13 on integral-fin tubes. They 
found that the root diameter had only a negligible effect on heat transfer enhancement. 
In addition, these investigators found that the optimal fm spacing was about 0.5 mm. 
They obtained an enhancement of 5 to 6. 

2 While the authors quote a fm spacing of 0.064 m, this value seems too small. Tlie schematic 
diagram of the fins together with other dimensions published by these authors tend to suggest that 
the correct fm spacing is 0.64 mm. This value is more in line with values quoted by other authors. 
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In comparing enhancements reported by various authors, a word of caution is 
required. In calculating values for the vapor-side heat transfer coefficients, the surface 
area to be used in this calculation must be defined. This area can be based on outside 
or envelope diameter, root diameter or total fin area. The great majority of authors 
define this area by the outside diameter of the tube as opposed to the root diameter. In 
the various sources cited above, this has been the case. That is, the surface area has been 
defined by the outside tube diameter. 

While it is clear that integral-fin tubes can result in an enhancement of the 
vapor-side heat transfer coefficient, other alternatives exist. As an example, Fujii et al. 
[Ref 33] studied the condensation of R-11 on smooth tubes which were WTapped with 
fine wire. These authors found an optimal pitch-to-wire diameter ratio of 2 which re- 
sulted in an enhancement of 2.5 in the vapor-side heat transfer coefficient when com- 
pared to copper smooth tubes. These authors argued that the enhancement in the heat 
transfer coefficient was due to a thinning of the condensate film. The reason that the 
condensate film thickness is reduced is that the presence of the wire creates a low' pres- 
sure region in the film near the base of the wires. This low pressure region is due to the 
concave shape of the condensate film next to the wire. As a result, condensate film is 
pulled along the tube w^all towards the wire. This results in a thinning of the condensate 
film between the wires. 

3. Summarj' of Single Enhanced Tube Condensation 

It is clear from the preceding discussion that theoretical models, despite their 
degree of sophistication, still only match experimental data to within ± 20%. This 
suggests that the process of condensation on enhanced tubes is extremely complex and 
far from being fully understood. 

In terms of experimental data involving condensation of refrigerants on en- 
hanced tubes, several comments need to be made. The experiments reported here have 
shown that integral-fin tubes can give enhancements of the vapor-side heat transfer co- 
efficient on the order of 7 to 12. The magnitude of this enhancement depends upon such 
factors as fin shape, fin spacing, and fin height. Trapezoidal shaped fins with a spacing 
of 0.5 mm give the highest enhancements. 

Unfortunately, the above experiments w'ere all performed on enhanced tubes 
made of copper. In Naval refrigeration systems, where seawater is used as the coolant, 
corrosion is of major concern. This necessitates the need for copper/nickel or titanium 
condenser tubes in. order to minimize corrosion. However, as can be seen from the 
foregoing summary of the literature, essentially no w'ork has been done on the 
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condensation of refrigerants on copper/nickel or titanium finned tubes. Since both 
copper/nickel and titanium condenser tubes have a lower thermal conductivity than 
copper, one must wonder what effect this will have on the vapor-side heat transfer co- 
efficient. Clearly, experiments using these tubes must be carried out. 

4. Condensation on Enhanced Tube Bundles: Theoretical Studies 

Condensation in tube bundles presents an entirely different situation than 
condensation on single tubes. In a tube bundle, two conflicting factors come into play. 
First, condensate forming on an upper tube drains off this tube onto lower tubes within 
the bundle. This increases the thickness of the condensate layer on lower tubes which 
further increases the resistance of these tubes to heat transfer. This effect is referred to 
as condensate inundation. 

Condensate falling from one tube to another can have a second effect. As 
droplets fall from one tube to the next, ripples and splashing (turbulence) are induced in 
the condensate film. This can actually enhance heat transfer on the lower tubes. 

Vapor velocity can also be a significant factor in tube bundle condensation. As 
the vapor travels into the tube bundle and condenses, the local velocity decreases. This 
decrease in local velocity leads to a decrease in the heat transfer coefficient similar to an 
inundation effect since vapor shear effects that thin the condensate film are reduced. 

In analyzing condensation in horizontal tube bundles, Kusselt [Ref 4] assumed 
that all condensate flowed from one tube to another as a continuous laminar sheet. He 
showed that the ratio of the average heat transfer coefficient for a vertical column of 
horizontal tubes with respect to the top tube could be expressed as: 



Kern [Ref 34] proposed a less conservative approach. He argued that droplets 
of condensate falling from one tube to another causes ripples which disturb the 
condensate film, thereby diminishing the effect of condensate inundation. His math- 
ematical model took the form: 
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where N is the number of tubes in the vertical column. A local coefficient can also be 
calculated for the Nth tube using Nusselt theory'. This coefficient is given by: 
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The local value of h for the Nth tube can be expressed as; 
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In 1972, Eissenberg [Ref. 35] proposed that condensate may not drain only in 
a vertical direction, but rather, may be diverted sideways, presumably due to local vapor 



At the present time, numerous investigations have been carried out in an at- 
tempt to better understand the phenomenon of condensate inundation. However, the 
data are scattered and as a consequence, no successful theoretical model of condensate 
inundation exists [Ref 3j. However, for design purposes, Butterworth [Ref 36] has re- 
commended that the Kern model be used. 

5. Condensation in Enhanced Tube Bundles: Experimental Studies 

Young and Wohlenberg [Ref 37] examined the condensation of R-12 on a 
vertical column of smooth horizontal tubes. These investigators found that the effect of 
the number of tubes in the vertical column is less than that predicted by Nusselt. 

Katz et al. [Ref 38] investigated the condensation of R-12 on a vertical column 
of six horizontal finned tubes. Their results also showed a less severe inundation effect 
than predicted by the model of Nusselt. 

Smirnov and Lukanov [Ref 39] used R-11 to study condensation on 20 hori- 
zontal rows of finned tubes arranged in a triangular fashion. Their results showed a 
dramatic reduction in heat transfer for the first five tubes. However, after this, the heat 
transfer coefficient remained almost constant for the rest of the tubes in the column. 

Honda et al. [Ref 40] have studied the condensation of R-113 on a vertical 
column of horizontal low finned tubes. An inundation tube was used to simulate addi- 
tional rows of tubes. These authors described four modes of condensate flow; namely, 
droplet mode, column mode, column and sheet mode, and sheet mode. Unfortunately, 



flow conditions. This so-called "side drainage" model predicts a less severe effect of 
inundation yields: 
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these authors did not calculate the heat transfer coefficient so that no comparison with 
theoretical models is possible. 

Webb and Murawski [Ref 41] examined the condensation of R-1 1 on a column 
of enhanced condenser tubes. The tubes tested were a standard 1024 fpi tube, a 
GEWA-SC tube, a Turbo-C tube and a TRED-D tube. The GEWA-SC tube is a 1024 
fpi tube but the fins have a Y-shaped cross-section. The Turbo_c and TRED-D tubes 
are similar to the 1024 fpi tube except that the fins have a saw-tooth shape. Webb and 
Murawski found that a standard 1024 fpm tube showed virtually no row (inundation) 
effect. However, when the experiments were conducted on Turbo-C, GEWA-SC and 
Tred-D tubes, row effects (ie., inundation effects) were seen. The inundation effect was 
greatest for the Tred-D tube which also had the lowest single tube performance. These 
investigators found values for the exponent in equations (2-17, 2-19, and 2-21) on the 
order of 0.12 to 0.26. Recall that Nusselt's theory used a value of 0.25. 

More recently, Honda et al. [Ref 42] studied the condensation of R-1 13 on an 
in-line horizontal tube bundle. The tubes tested consisted of either flat-sided (rectangu- 
lar) annular fins or three dimensional trapezoidal fins. These authors found that the ef- 
fects of condensate inundation were greatest for the three dimensional fm tubes. 
Unfortunately, these authors also had vapor velocities on the order of 3.4 m/s or greater. 
Therefore, it is difficult to compare their results directly to the theoretical models of 
Nusselt, Kern, or Eissenberg. 

An interesting facet of the experiments of Honda et al. [Ref 42] deserves men- 
tion. These investigators compared the effects of condensate inundation on the vapor- 
side heat transfer coefficient for smooth and finned copper tubes. Their results showed 
that despite condensate inundation, the heat transfer coefficient for the finned tubes was 
still substantially above that of the unfinned tubes. In addition, the decrease in the heat 
transfer coefficient with increasing row number w’as less with the finned tubes suggesting 
that the effects of condensate inundation were less severe with finned tubes than with 
unfinned tubes. Presumably, this result is due to the fact that condensate draining off 
of a finned tube drains in columns rather than continuous sheets. Lower tubes will 
therefore not be completely covered with condensate but instead will have unflooded 
regions. 

Honda and Nozu [Ref 43] examined the effects of condensate inundation on 
smooth and integral-fin copper tubes using a theoretical model derived recently by 
Honda et al. [Ref 44] This model was a single tube model which was extended to include 
the effects of condensate inundation by realizing that condensate did not flow off tubes 
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as sheets but rather as columns. The model took account of the fact that regions which 
received drainage of condensate had different heat transfer coefficients than regions 
which did not. The results support the more recent experimental results of Honda et al. 
[Ref 42]. The effect of inundation was less severe for the finned tubes than for the 
smooth tubes. The authors suggest that this is due to condensate being charmeled on the 
finned tubes whereas for the smooth tubes, the entire tube surface becomes covered with 
a condensate film of varjdng thickness. 

6. Summary' of Enhanced Tube Bundle Condensation 

Experimental data has shown that the Nusselt model for the effects of 
condensate inundation may be too conservative. Rather, most experimental data seems 
to fall between the Kern and Eissenberg predictions. 

Other experimental data suggests that the effects of condensate inundation are 
less for finned tubes as opposed to smooth tubes. Presumably, this is due to the manner 
in which condensate flows from one finned tube to another. The presence of the seems 
acts to channel the flow of condensate leaving large areas relatively free of flooding from 
inundation. 

As a final comment, the use of wire wrapped tubes provides an interesting al- 
ternative to the more conservative finned tube approach. However, the effects of 
refrigerant condensate inundation on wire wrapped tubes has yet to be examined. 
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III. EXPERIMENTAL APPARATUS 



A. TEST PLATFORM 

Details of the design and construction of the condenser/evaporator test platform are 
provided by Zebrowski [Ref 1] and Mabrey [Ref 2]. Therefore, only a brief summary 
will be presented here. 

The test platform is shown schematically in Figure 1 while a photograph of the unit 
is showm in Figure 2. The unit is made from rolled stainless steel plate (6.35 mm thick) 
and designed to withstand an absolute pressure of 308 kPa (approximately 3 bar). 

I. Condenser Section 

The top section is the condenser. It is 1.3 m in length with an internal diameter 
of 0.61 m. The end plates are made of stainless steel plate (6.35 mm thickness) with a 
diameter of 0.71 m. They are connected to the condenser shell using a bolted flange as- 
sembly. A rubber gasket, 3.22 mm thick, provides a vacuum tight seal. The condenser 
unit has five viewports which allow the condensation process to be viewed. 

Four test condenser tubes are contained in the condenser forming a vertical in- 
line column with a pitch-to-diameter ratio of 2.25. These are held in place at each end 
by a nylon block bolted to the end plate. A rubber gasket provides a vacuum tight seal. 
A stainless steel plate is fitted over the ends of the condenser tubes and bolted to the 
nylon block. O-rings are placed in between the nylon block and steel plate in order to 
provide a vacuum tight seal. 

In order to provide a means of pressure control during testing of the condenser 
tubes, the test platform was fitted with five auxiliary condenser tubes. These auxiliary 
condenser tubes are constructed of 9.53 mm diameter copper tube, wound into coils. 
Each coil is supported by a stainless steel rod welded to one of the end plates. The aux- 
iliary condenser coils enter and exit through this same condenser end plate. 

A stainless steel shroud assembly, shown in Figure 3, is installed in the 
condenser unit. This shroud surrounds the test condenser tubes as w’ell as the auxiliary 
condenser coils. A window, installed in the side of the shroud, allows visualization of the 
test condenser tubes through the viewports. The purpose of this shroud is to funnel va- 
por from the evaporator below, around the outside of the shroud and then vertically 
down past the condenser tubes. The velocity of the vapor past the condenser tubes was 
0.1 m/s. 
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Figure 1. Schematic diagram of test platform. 
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Figure 2. Photograph of the evaporator/condenser test platform. 
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Figure 3. Schematic diagram of condenser shoving shroud. 
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2. Evaporator Section 

The flooded evaporator is formed from rolled stainless steel and has an internal 
diameter of 0.61 m with a width of 0.279 m. Two viewports allow visualization of the 
boiling process. This unit is connected to the condenser by a collar assembly 0.28 m in 
diameter and 0.203 m in length. The volume of liquid contained in the evaporator is 
approximately 0.06 m^, just sufficient to cover the boiler tubes. 

Vapor is produced in the evaporator using a combination of three groups of 
heaters. The first group is the main bundle and is composed of thirty-five tubes in a 
triangular pitch (pitch to diameter ratio of 1.2). Fifteen of these have heaters with a 
nominal rating of 1 kW each. Five of these heated tubes are instrumented with 
thermocouples to enable bundle boiling data to be obtained. In addition, unheated 
dummy tubes are used to guide the two phase flow within the bundle. 

The second group of tubes is composed of five simulation heaters, each with a 
nominal rating of 4 kW. These are placed below the main bundle in order to simulate 
further rows of tubes. 

The final group consists of four auxiliary heaters, each with a nominal rating 
of 4 kVV. Two of these are placed on each side of the evaporator in order to keep the 
liquid pool at saturation conditions; this avoids any subcooling of the liquid. Power to 
all three sets of healers is provided by three separate STACO 240V, 23.5kVA rheostat 
contollers. 

For further information on the design and construction of the test platform, the 
reader is referred to Zebrowski [Ref 1], Mabrey [Ref 2], Murphy [Ref 45], and Chilman 
[Ref 46]. 

3. Coolant System 

The coolant is a 54/46 by weight mixture of ethylene glycol and water. The 
coolant system consists of two main coolant pumps. These pumps are 0.5 hp constant 
speed pumps which take a suction from the main coolant sump. One pump is used to 
supply five rotameters. A ball valve, located at the entrance to each rotameter, is used 
to control the flow rate. Flow to each test condenser tube is supplied by an individual 
rotameter. The fifth rotameter is used to supply one of the auxiliary condenser coils. A 
needle valve is used to provide fine control of flow through the coil. This coil is normally- 
used to control saturation pressure during experimental runs. 

The second coolant pump is used to supply coolant flow to the remaining four 
auxiliary condenser coils. Flow through each coil is controlled by a needle valve located 
at the entrance to each coil. It should be noted that the four auxiliary cooling coils are 
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seldom needed during the condenser experiments; their primary purpose is for use during 
evaporator experiments, (see Chilman [Ref. 46]). 

In addition to supplying the four auxiliary condenser coils, the second coolant 
pump also provides coolant flow to a refrigerant storage tank and to a cold trap located 
at the top of the condenser. The storage tank can be used to store refrigerant during 
condenser or evaporator tube changeout. The cold trap is placed in-line between the 
condenser and a vacuum pump. Vapor from the condenser is drawn through the cold 
trap which is cooled by the ethylene glycol/water mixture. This process is used to remove 
noncondensibles from the condenser wliile preventing loss of refrigerant to the atmos- 
phere. 

Flow of the coolant mixture from the auxiliar)' condenser cooling coils, storage 
tank, and cold trap are returned to the sump tank. Coolant flow from the test condenser 
tubes empties into a common header which then empties into the same coolant sump 
tank. 

In order to keep the coolant mixture temperature approximately constant, it is 
cooled by an eight ton refrigeration system located outside the building. This system 
continually recirculates the coolant using a 0.75 hp pump. Coolant temperature can be 
maintained at any temperature from ambient temperature (approximately 20 °Q down 
to -20 °C. 

4. Instrumentation 

Prior to entering each condenser tube, the coolant flows through a mixing 
chamber. The inlet temperature is measured using a single copper/constantan type-T 
thermocouple. Coolant flows from each condenser tube through an outlet mixing 
chamber. Outlet coolant temperatures are then measured using two thermocouples for 
each tube. An average outlet temperature is then taken for each tube. All thermocouples 
had been previously calibrated by Mabrcy [Ref 2]. 

Refrigerant temperatures were measured at several different locations in the 
evaporator/condenser test platform using similar single thermocouples. Vapor temper- 
ature was measured in two separate locations at the top of the condenser and one lo- 
cation just above the liquid level. Refrigerant liquid temperature was measured at both 
the top and bottom of the tube bundle. The thermocouples were inserted into stainless 
steel housings (with copper tips) which were installed in the condenser/evaporator test 
platform. 
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B. CLEANING AND DECONTAMINATION PROCEDURES 

As mentioned earlier, there was a contamination problem with the test platform 
during the work of Mabrey [Ref. 2]. Therefore, the experimental apparatus was initially 
disassembled and cleaned in accordance with the following procedure. 

Prior to removing the evaporator, all heaters and thermocouples were disconnected. 
The lagging around the mating flange was then removed. A Porta-lift was connected to 
the evaporator mating flange in order to support the weight of the evaporator prior to 
removing the bolts holding the evaporator and condenser together. After supporting the 
evaporator, the bolts wxre removed and the evaporator lowered from the condenser unit. 

Referring to Figure 2, the condenser unit is capped at each end by a steel end plate. 
The left end plate supports the au.xiliary condenser units and is therefore removed after 
the right end plate. The test condenser tubes must be removed prior to the disassembly 
of the condenser. The assembly which holds the condenser tubes in place is composed 
of a nylon plate bolted to the condenser end plate. A gasket insures a tight seal betw'een 
the nylon and condenser end plates. A stainless steel plate is bolted against the nylon 
plate with 0-rings in between w’hich ensure an air tight seal between the environment 
and the inside of the condenser. 

The plates holding the test condenser tubes in place must be removed prior to re- 
moval of the condenser end plates. The stainless steel plates at both ends are removed 
by unbolting the retaining nuts. The nuts holding the nylon plate in place are then re- 
moved. A rubber mallet is used to tap out the test condenser tubes from the left end of 
the platform. This causes the right end nylon plate to move away from the right 
condenser end plate. Screw'drivers can then be used to pry the nylon end plate away from 
the condenser. 

The condenser end plates can be removed after the test condenser tubes have been 
removed. As a word of caution, these end plates weigh in excess of 100 pounds. 
Therefore, they must be supported during the removal process. In these studies, a 
Porta-lift w’as used for this purpose. The right condenser end plate was removed first, 
followed by the left end plate. Prior to removal of the left end plate, the shroud assembly 
had to be removed. This was facilitated by first removing the thermocouple probes which 
penetrate the top of the condenser unit. The shroud W'as then pulled out of the 
condenser. 

If the auxiliary condenser tubes are to be removed, then the the left condenser end 
plate will have to be removed. This is because the five auxiliary condenser tubes are at- 
tached to the left end plate. The location of the test platform in relation to the coolant 






sump is such that this end can not be removed with the evaporator/condenser test plat- 
form in its present location. Rather, if the left head must be removed, the entire test 
platform must be rotated 180 degrees about a vertical axis. 

Initial inspection of the condenser unit showed only one small area of discoloration 
which surrounded the middle viewport. No other obvious signs of contamination were 
evident. The condenser and evaporator were both cleaned using "3-D Supreme Cleaning 
Agent", a commercially available product which leaves no residue. All components were 
rinsed with fresh water, dried, then given a final cleaning with acetone. The auxiliary 
condenser coils were cleaned and decontaminated in the same manner. 

The copper condenser tubes also appeared somewhat discolored. These were cleaned 
using "Shower Power", a commercially available cleaning agent. Prior to installation into 
the condenser unit, they were given a final cleaning with acetone to ensure removal of 
any residue left over from the cleaning procedure. 

During the reassembly, all gaskets were replaced. Prior to replacement of the 
evaporator unit, a test blank was installed on the condenser and a vacuum test success- 
fully performed. After this test was completed, the test blank was removed and the 
evaporator reinstalled. All penetrations were reconnected, the boiling tube bundle 
reinstalled and another vacuum test performed. 

Prior to charging the system with R-113, the ethylene glycol/water coolant pumps 
were energized and flow started through the test condenser tubes. For this initial test, 
smooth copper tubes were used. The temperature of the coolant was -15°C. Since 
Mabrey [Ref. 2] had noted that during condensation, a thickening between the fins of 
copper/nickel finned tubes was occurring which suggested ice formation, it was felt that 
water vapor might be present in the test platform since it had been open to the atmos- 
phere. Therefore, by running the coolant mixture at this temperature through the 
condenser tubes, if no condensation was seen, this would rule out the presence of water 
vapor. During this intitial run, no condensation was seen on any tubes. In addition, the 
orange peel appearance noted by Mabrey on the smooth copper tubes was not in evi- 
dence. 

After checking the evacuated system for possible contamination, the test platform 
was charged with fresh R-113. Coolant flow was again started. This time an orange peel, 
or dimpled appearance could be seen on the top condenser tube. This suggested some 
type of contamination layer. 
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Power (~ 20 kW) was applied to the heaters while the coolant temperature was 
raised to ~ O’C The orange peel appearance on the top tube was still present. Meas- 
urements of the vapor temperature at the top and bottom of the condenser showed ap- 
proximately a 20 ®C differential, signifying the presence of significant amounts of 
noncondensible gases. These gases enter with the refrigerant and are released into the 
condenser when the liquid is boiled. 

A vacuum pump, connected to the top of the condenser was used to remove these 
noncondensible gases. When this was done, the temperature differential between the top 
and bottom of the condenser decreased to less than 0.5° C. The orange peel appearance 
quickly disappeared and all four condenser tubes demonstrated good filmwise 
condensation. Based on these observations, it was standard practice to degas the system 
for at least four hours after the system was open to the atmosphere or fresh R-113 was 
added to the evaporator. 
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IV. EXPERIMENTAL PROCEDURES 



A. GENERAL PROCEDURES 

The experimental procedure was identical for all experiments. The coolant pumps 
were energized and flow started through one of the auxiliary condensers. Flow was also 
started through the four test condenser tubes at a nominal coolant velocity of 0.2 m/s. 
The auxiliary heaters were set at 3 kW, the instrumented tube bundle set at 3.5 kW, and 
the simulation heaters set at 1.5 kW. 

The vapor temperature at the bottom and the top of the condenser was allowed to 
rise to 47.5 which corresponds to the saturation temperature of R-1 13 at atmospheric 
pressure. The rate of temperature rise was controlled to approximately 0.3 degrees per 
minute by adjusting the rate of coolant flow through the auxiliary condenser. During the 
heat up period, a vacuum pump, connected to the top of the condenser via the cold trap, 
was run continuously in order to remove any noncondensibles. The absence of 
noncondensibles was indicated when the vapor temperature difference between the top 
and bottom of the condenser unit was less than 0.5 °C. 

When the vapor temperature was about 45 "C, coolant velocity was adjusted to be- 
tween 1 and 1.2 m/s for the four condenser tubes. The heater power levels were then 
readjusted in order to maintain a relatively constant condenser pressure while leaving 
one auxiliary condenser coil on-line. At lower coolant velocities, condenser pressure was 
maintained using the auxiliary condenser coils. The power settings for the heaters were 
not altered during the course of the experiment. 

During all experimental runs, the vapor temperature at the top and bottom of the 
condenser unit was maintained between 46.5 and 48.5 °C which corresponds to a 
condenser pressure of approximately one atmosphere. Pressure control, was accom- 
plished by adjusting coolant flow through the auxiliary condenser tubes. 

1. Tube Bundle 

Once the vapor temperature was between 46.5 and 48.5 °C , coolant flow was 
adjusted to the desired coolant velocity for the four test condenser tubes. Here, the de- 
sired coolant velocity is an input to the computer program. The program then calculates 
the required flowmeter setting using the calibration curves (Appendix A). The initial 
coolant velocity used was 0.2 m/s. After setting the flow rate, the temperature difference 
between the inlet and outlet of all four condenser tubes was monitored using the com- 
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puter program, DRPCON7. Approximately ten minutes was required for this temper- 
ature differential to stabilize. Once stable, a data set was taken. This procedure was 
repeated in increments of 0.1 m/s up to a maximum of 1.2 m/s. The procedure was then 
repeated by decreasing the flow in 0.1 m/s increments down to 0.2 m/s. Two data sets 
were obtained at each coolant velocity. It should be mentioned that for the 
copper/nickel finned tubes, the range of flow settings was changed. The minimal flow 
setting was 0.6 m/s while the maximum flow setting was 2.0 m's. The reasons for this 
change are discussed in “3. Enhancement of Outside Heat Transfer Coefficient.” on 
page 30. 

2. Single Tube Experiments 

After completion of the tube bundle runs, the experiment was repeated for each 
individual tube starting with the lowest tube (tube D). Flow through the tubes above the 
test tube was stopped and flow through the tube of interest was set for a velocity of 0.2 
m/sec. The vapor pressure was again maintained using an auxiliary condenser coil. 
During this time, the upper tubes were continuously observed for condensate formation. 
No single tube data w^as obtained until condensate formation on the upper tubes had 
ceased. This typically required twenty to thirty minutes. 

Once the coolant temperature rise had stabilized, a data set was obtained. .Fol- 
lowing this, the coolant flow rate was increased in 0.1 m/s increments up to the maxi- 
mum flow velocity. A data set was obtained at each flow velocity. Throughout the 
experiment, tubes above the test tube were visually monitored to ensure that no 
condensate formation occurred. The flow velocity was then decreased in 0.1 m/s incre- 
ments down to the minimum coolant flow velocity, with a data set obtained at each flow 
velocity. This procedure was repeated for each test condenser tube in the bundle. 

B. SPECIFIC EXPERIMENTS 
1. Baseline Experiments 

The purpose of the baseline experiments was to validate the performance of the 
condenser/evaporator test platform by comparing data to that of Mabrey [Ref 2]. These 
initial baseline experiments were performed using smooth copper tubes (inside diameter 
of 13.26 mm , outside diameter of 15.88 mm, thermal conductivity of 386 W/m K). The 
condensing length of the tubes was 1.219 m. The tubes were supplied by Wolverine Tube 
Co. 
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2. Enhancement of Inside Heat Transfer Coefficient 

The main purpose of this study was to obtain accurate values of the outside heat 
transfer coefficient, h,. In order to obtain accurate values for h^, the resistance to heat 
transfer on the inside of the tube - 7 -^ needs to be minimized. The reason that this re- 

h,Ai 

sistance must be minimized is that if the inside resistance is a large part of the overall 
vapor-to-coolant resistance, then will be become inaccurate. In order to minimize the 
inside resistance, one must enhance the inside heat transfer coefficient. 

The inside heat transfer coefficient is influenced by the type of coolant flow 
through the tubes. In the test platform, the coolant is a 54/46 percent (by weight) 
ethylene glycol/water mixture. At temperatures below 0 'C, this fluid is extremely 
viscous with Reynolds numbers in the range of 300 to 2800, suggesting laminar How. 
When the coolant flow is laminar, the inside heat transfer coefficient is low. However, 
if this laminar flow can be mixed, then the inside resistance can be decreased and the 
heat transfer coefficient increased. This can be accomplished through the use of radial 
mixing elements. Radial mi.xing elements have the effect of mixing the thermal boundary 
layer (formed on the inside of the tube) with the ■^older core flow thereby improving heat 
transfer. It should be noted that a penalty in _ .ng such elements is an increased pres- 
sure drop which must always be kept in mind in practical applications. However, for the 
purposes of this study, increased pressure drop was not important. 

In this thesis, two types of mixing elements w^ere investigated. The first was a 
twisted tape insert made of copper [Ref 47,48]. This tape insert had a thickness of 0.559 
mm with a pitch for a 180 degree twist of three times the inner tube diameter. 

The second type of insert used in this study w'as a commercially available prod- 
uct referred to as HEATEX radial mixing elements.3 These elements are manufactured 
from stainless steel and consist of a number of small flexible loops (petals) of steel wire 
attached to a central wire core. These flexible petals are placed in the tube such that they 
face the oncoming flow'. Their outside diameter is slightly larger than the inside diameter 
of .the tube and are therefore kept in position with a "pinch" fit. The petals disturb the 
thermal boundar>’ layer and promote fluid mixing. Figure 4 shows the two types of in- 
serts. It should be noted that two different sizes of HEATEX were used since condenser 
tubes of two different inside diameters were examined. 



3 manufactured by CAL GAVIN Birmingham, England. 
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Figure 4. Twisted tape insert and HEATEX radial mixing element. 
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As pan of the baseline experiments, the smooth copper tubes were also tested 
with no insert installed. Experiments were repeated on the smooth copper tubes using 
the tw'isted tape inserts as well as the HEATEX elements. This data was then used as a 
baseline for comparison with other types of condenser tubes. It should be mentioned 
that w’hen no insert is present, the inside tube resistance to heat transfer is ver}' large 
relative to the outside resistance. In some cases (typically when the coolant flow rate is 
low), this can produce a negative value for the outside heat transfer coefficient. This is 
simply the inaccuracy in the calculation of w’hen the inside resistance dominates and 
clearly demonstrates the need for an insert as mentioned above. 

3. Enhancement of Outside Heat Transfer Coefficient. 

The overall heat transfer coefficient, (4, is influenced by the outside as well as 
the inside heat transfer coeflicient. The means by which the outside heat transfer coeffi- 
cient is enhanced is by altering the outside geometry of the condenser tube. 

Experiments to investigate enhancement of the outside heat transfer coefficient 
utilized three types of condenser tubes. The first type of tube tested was a copper/nickel 
finned tube (26 fins/inch). This tube has an inside diameter of 10.16 mm with a thermal 
conductivity of 43 W/m K. The tubes w'ere manufactured by Wolverine Tube Co. These 
tubes W'ere tested as a tube bundle as well as individually. Experiments were performed 
with no insert as well as with the twisted tape and HEATEX inserts. These tubes are 
show'n in Figure 5. One point must be mentioned with regard to the finned tubes. The 
use of fins on the outside of the tube enhances the heat transfer considerably. Conse- 
quently, the inside resistance again becomes a large fraction of the overall resistance, 
even W'hen inserts are used. 

One additional point must be discussed with regard to the copper/nickel finned 
tubes. Because the size of the holes in the nylon end plates is fixed, the outside diameter 
of any tubes to be tested must be the same. In order to have these finned tubes the same 
outside diameter, the inside diameter had to be reduced. Since the inside tube diameter 
was smaller for the finned tubes, higher coolant velocities were obtained. 

The second series of experiments used a copper/nickel KORODENSE tube also 
manufactured by Wolverine Tube Co. This tube, shown in Figure 5, has a shallow spiral 
groove along the entire condensing length with a pitch of 10 mm. As with the finned 
copper/nickel tubes, experiments w'ere performed on these KORODENSE tubes with no 
insert as v/ell as with the twisted tape and HEATEX inserts. The tubes w'ere tested in a 
tube bundle as well as individually. 
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Fujii and co-workers (Ref. 33] have previously shown that the condensation of 
refrigerants can be enhanced by wrapping smooth condenser tubes with a small diameter 
wire. Therefore, a third scries of experiments were performed in which copper/nickel 
KORODENSE tubes were tightly wound with wire as shown in Figure 5. Since these 
authors showed that there was an optimum wire size wliich gave the greatest enhance- 
ment of heat transfer (0.5 mm at a pitch of 3 mm), the KORODENSE tubes used in this 
study were wrapped with one of three different sizes of wire. The sizes used were 0.029" 
(0.736 mm) diameter stainless steel wire, 0.049" (1.244 mm) diameter stainless steel wire, 
and 0.0675" (1.715 mm) diameter copper wire. The wire lay in the shallow grooves of the 
KORODENSE tubes and was soldered to the tubes at the ends only. Although the wire 
was tightly wrapped around the tubes, it was not physically attached except at the ends. 
Experiments were carried out on each tube individually. Each tube was tested in all four 
locations within the bundle. This was done to see if position within the bundle (i.c. wake 
from tubes above) had any effect. Only the HEATEX inserts were used with these tubes. 

Once an optimal wire diameter was determined, a fourth series of experiments 
was conducted. This series used the same four KORODENSE tubes, only this time each 
was wrapped with the same optimal diameter wire. The tests were done with the 
HEATEX elements installed. Experiments were conducted on the tubes individually in 
order to compare with the previously obtained single tube data. In addition, these wire 
wrapped tubes were tested as a tube bundle in order to investigate whether or not the 
wire wrap had any effect on condensate innundation. 
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Figure 5. Copper/nickel finned tube, Korodense tube, Wire wrapped tube. 
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V. DATA REDUCTION 



A. DATA ACQUISITION 

The data acquisition system consisted of a Hewlett-Packard series 9000 computer 
coupled with a Hewlett-Packard model 3479A Data Acquisition system. The computer 
program used in this study w’as originally written by Mabrey [Ref. 2] but has been sub- 
stantially modified in the course of this thesis. The data acquisition thermocouple 
channel assignments are shown in Table 1. 



Table 1. THERMCOUPLE CHANNEL ASSIGNMENTS. 



Channel 


Temperature Measured 


0 


Vapor (condenser top) 


1 


Vapor (condenser top) 


— 2 


Vapor (top of liquid) 


3 


Liquid (top of bundle) 


4 


Liquid (top of bundle) 


5 


Liquid (bottom of pool) 


6 


Inlet tube 1 


7 


Inlet Tube 2 


8 


Inlet Tube 3 


. 9 


Inlet Tube 4 


10, 11 


Outlet Tube 1 


12, 13 


Outlet Tube 2 


14, 15 


Outlet Tube 3 


16, 17 


Outlet Tube 4 



The computer program allows the user to acquire data, reprocess data, or plot data. 
During data acquisition, the user selects the fluid velocity desired. Based on coolant inlet 
temperature, the coolant fluid properties and mass flow rates are calculated (see Ap- 
pendix C). The program then uses equations derived from the flowmeter calibrations 
(Appendix A) in order to calculate the appropriate flowmeter settings. 

Once the flowmeters are set for the required coolant velocity, the computer system 
displays the vapor temperature at the top and bottom of the condenser as well as the 
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coolant temperature rise for each condenser tube. When these temperature rises are 
stable (approximately 10 minutes), the user can terminate the display and acquire the 
data for that coolant flow velocity. While the temperature rises were stabilizing, coolant 
flow through an auxiliary condenser tube was adjusted in order to maintain constant 
pressure. 

B. DATA REDUCTION 

The main function of the data reduction section of the computer program is to back 
out the outside heat transfer coefficient from an experimentally determined value of U, 
and an assumed inside heat transfer coefficient correlation. This is accomplished in the 
following manner. 

The data acquisition system measures the emf from the thermocouples listed in Ta- 
ble 1. These values of emf are converted to temperatures using correlations obtained 
from the temperature calibrations both contained in Appendix B. The volume flow rate 
of the ethylene glycol/water mixture is calculated based on the desired coolant velocity. 
The density, dynamic viscosity, thermal conductivity, heat capacity, and Prandtl number 
of the coolant are calculated using inlet temperatures (see Appendix C and Appendix 
E). From these quantities, the coolant mass flow rate and fluid velocity are calculated. 
The coolant Reynolds number for each condenser tube is also calculated. The heat 
transferred to the coolant is then found from 



where m is the mass flow rate and AT is the coolant temperature rise. The outside heat 
flux, q'\ is calculated by dividing q by the surface area of the tube based on the outside 
tube diameter. 

The overall heat transfer coefficient, f/'o, is given by: 



q = mCpAT 



(5-1) 



LMTD 



<7 0 



(5-2) 



where LMTD is the log mean temperature difference defined as: 



LMTD 



AT 



(5-3) 
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The inside heat transfer coeflicient (given as a dimensionless Nusselt number, Nu^) 
can then be calculated using one of three correlations. The choice of which correlation 
to use depends upon the type of insert installed inside the condenser tube. If no insert 
is installed, then Nu^ can be calculated from a correlation provided by Hausen [Ref 
49]'* for developing laminar flow in a tube at constant wall temperature. 



Nu, = 3.66 + 



0.0668( Y \ReD,cPrc 



1 + 0.04 



\ 



0.666 



■ (5-4) 



If a twisted tape insert is installed, a correlation given by Hong and Bergles [Ref 
47] can be used to calculate Nu, where 



Nu, = 5.172 



1 + 5.484xlO~Vr?’ 



Re^ 

~y 



1 . 25 ' 



0.5 



(5-5) 



Here, y is the twist ratio (pitch over diameter). Re, is the coolant Reynolds number, 
corrected for tape thickness. Re, is given by: 



Re^ 



Am 

7Zfx{di — AS) 



(5-6) 



where S is the tape thickness. This correlation is valid for Re, up to about 2300. 

A number of correlations for determining the inside heat transfer coeflicient have 
been published when using a wire matrix insert [Ref 50,51,52, 53 ]. However, these 
correlations are specific for the type of mixing element used, and are especially depend- 
ent upon the loop density of the wire matrix. 

The HEATEX mixing elements supplied by CAL GAVIN (nominal diameters of 
13.26 and 10.16 mm) did not match the loop densities used by previous investigators. 
However, CAL GAVIN did supply experimental data for the heat transfer coefficient 
versus fluid velocity at three different temperatures for both sizes of elements. To con- 
vert this data into a usable correlation, the following procedure was developed. 



4 Calculations were performed for coolant flow through the smooth copper tubes at coolant 
velocities of 0.2 and 1.2 m/s. These calculations confirmed the coolant flow was not fuUy devel- 
oped. 
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The correlations of Oliver and Aldington [Ref. 51, 52 ] took the form of 

Nu = cRe'^Pr'' (5 - 7) 

where c, m, and n are arbitrary constants which are dependent upon the type of insert 
used. It was assumed that a similar type of equation could be used in the present study 
by fitting the constants to the data (once non-dimensionalized) provided by CAL 
GAVIN. 

The correlations of Ohver and Aldington [Ref 51, 52 ] had Prandtl number raised 
to a power of 0.46. This w'as primarily to take into account temperature effects. There- 
fore, the same value was assumed in the present study. Values of Nu/Pr®'*‘ were plotted 
as functions of Re on a log-log basis. Figure 6 shows the data obtained for the 13.26 
mm elements while Figure 7shows the data obtained for the 10.16 mm elements. Based 
on the appearance of the non-dimensionalized data, a least squares linear regression was 
performed on each data set. The figures yield the values of the constants m and c (slope 
and intercept respectively) and provide the following correlations for the 13.26 and 10.16 
mm elements respectively: 

Nu, = 0.226Pe°-^^Pr°-^^ (5 - 8) 

and 

Nu, = 0.065Pe° ’Vr° "^ (5 - 9) 



Once Nu, was obtained from equations (5-4, 5-5 5-8 or 5-9), the inside heat transfer 
coefficient was calculated from: 



hi = 



Nu,A, 

~d~ 



(5-10) 
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HEATEX CORRELATION (13.26 MM) 



NU = 0.2E6(RE**.65)(PR**.48) 




Figure 6. HITRAN correlation for 13.26 mm elements 
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HEATEX CORRELATION (10.15 MM) 

NU = 0.067(RE**.76)(PR**.46) 




Figure 7. HITRAN correlation for 10.16 mm elements 
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The overall heat transfer coefficient, U, is found from the summation of resistances 
to heat transfer across the Vi’hole tube: 



1 



1 

^0^0 



+ Rj- + Ry^, + 



1 

hiAi 



(5-11) 



where and R^ are the fouling and wall resistances respectively. In order to find 
equation (5-11) can be rearranged to give: 

K * 7-r-r (5-12) 



In all experiments, the tubes were meticulously cleaned and therefore, the fouling re- 
sistance, Rp was neglected. 

In order to better assess the calculated experimental value of the outside heat 
transfo coefficient, some form of comparison with prediction must be made. In this 
case, the experimentally obtained value of was compared with Nusselt theory. How- 
ever, because the wall temperatures are not known a priori, an iteration must be used 
to determine this predicted value of h„. An initial value of h, was estimated and based 
on this estimate, a wall temperature was calculated using the equation: 



T — T — 
^ wall ^ sat 




(5-13) 



The conden.sate film properties are then calculated using a reference temperature given 
by: 



1 7 

T = — T 4- — T 
^film 2 ^ sat '3 ^ ^ 



wall 



(5 - 14) 



Once these condensate fluid properties are known (Appendix D), the theoretical outside 
heat transfer coefficient can be calculated from Nusselt theory' (Ref 4): 



K 



0.655 



gk]p]hf^ 

o 



1 

3 



(5-15) 



This calculated value is then compared to the initial guess. The error between the initial 
guess and the calculated value are compared. If the error is greater than 0.1%, the 
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"guessed" value for h, is updated and the iteration continued. This process is repeated 
until the error between the calculated and guessed values is less than 0.1%. 

C. MODIFIED WILSON PLOT 

The outside and inside heat transfer coefficients could be determined directly by 
measuring the vapor temperature, the tube wall temperature, and the coolant temper- 
ature rise. Of these, the most difficult temperature to measure directly is the tube wall 
temperature. This would require the use of instrumented tubes with thermocouples em- 
bedded in the tube wall. Unfortunately, the manufacture of such instrumented tubes is 
both costly and time consuming, especially when large numbers of tubes are to be tested. 

An alternative to the use of instrumented tubes is to employ the "Modified Wilson 
Plot" procedure. This procedure has been outlined in detail by Marto [Ref. 10]. The 
modified Wilson plot procedure relies upon the overall heat transfer coefficient being 
reliably measured from experimental data. A summation of heat transfer resistances 
from vapor-to-coolant is then expressed as: 



As mentioned above, the fouling resistance is assumed to be negligible. This summation 
of resistances can be rearranged into a linear equation: 



In looking at the correlations given above for /i, and h„ it can be seen that they all 
employ a leading coefficient. As an example, Nusselt theory’, (equation 5-15) can be ex- 
pressed in the form: 




(5-16) 




(5-17) 



ho = O.F 



(5-18) 



where: 



shj Pfhfg 3 
o 



(5-19) 



In a similar manner, the equations for the inside heat transfer coefficients (equations 5-4, 
5-5, 5-8, and 5-9) can also be transformed into the form: 



h; = CQ (5 - 20) 

where the form of Q would depend on which type of tube insert was being used. As an 
example, using equation 5-5, Q would become; 

0.5 

(5-21) 



Q = -f 



1 -1- 5.4S4A-10"Vr° ’ 



Re^ 

~y~ 



1.25 



The linearized version of the summation equation (equation 5-17) can then be rear- 
ranged to give: 



Tliis is a linear equation with two unknowns, namely, C and a . Letting 



(5-22) 



(5 - 23) 



and 



X = 



AqF 

A^Q 



we can obtain a simplified linear equation: 

r=mX+b 



where 



and 





(5 - 24) 



(5 - 25) 



(5-26) 



(5-27) 



Since £2 and F are both temperature dependent, an iterative solution procedure must be 
used. This is accomplished by applying a least squares fit to the linear equation. 

It should be pointed out that the accuracy of the modified Wilson plot procedure 
relics heavily on the number of data points taken as well as the range of coolant veloci- 
ties used. 
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VI. RESULTS AND DISCUSSION 



A summary of the experimental runs performed in this thesis is presented in 
Table 2. The abbreviation "SM" is for the smooth copper tubes, "KD” is for the 
KORODENSE tubes, and "CF" is for the copper/nickel finned tubes. 



Table 2. SUMMARY OF EXPERIMENTAL RUNS 



Data Set 


T ube type 


Insert 


Comments 


SMOl 


Copper 


Tape 


Test run 


SM02 


Copper 


Tape 


Test run 


SM03 


Copper 


HEATEX 


Test HEATEX correlation 


SM05 


Copper 


HEATEX 


Tube A insert bad 


SM06 


Copper 


HEATEX 




SM07 


Copper 


Tape 




KDOl 


KORODENSE 


Tape 


1st test of KORODENSE 


KD02 


KORODENSE 


None 




KD03 


KORODENSE 


HEATEX 




KD04 


KORODENSE 


Tape 




KD05 


KORODENSE 


HEATEX 


Wire wrapped 


KD06 


KORODENSE 


HEATEX 


Wire wrapped 


KD07 


KORODENSE 


HEATEX 


Wire wTapped 


KD08 


KORODENSE 


HEATEX 


Wire wrapped 


KD09 


KORODENSE 


HEATEX 


0.049" Wire, bad insert 


KDIO 


KORODENSE 


HEATEX 


0.049" Wire 


KDll 


KORODENSE 


HEATEX 


0.049" Wire 


KD12 


KORODENSE 


HEATEX 


0.049" Wire 


CFOl 


Cu/Ni fin 


Tape 


dropwise condensation 


CF02 


Cu/Ni fin 


Tape 


Test run 


CF03 


Cu/Ni fin 


Tape 




CF04 


Cu/Ni fin 


HEATEX 




CF05 


Cu/Ni fin 


None 
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A. DECONTAMINATION 

Prior to charging the system with R-1 13, the system was evacuated to approximately 
15 iiiHg and left for 24 hours. For this initial test, smooth copper tubes were used. The 
ethylene glycol coolant pumps were energized and flow' started through the test 
condenser tubes. The temperature of the ethylene glycol w'as -IS'C. No condensation 
was seen on any tubes. In addition, the orange peel appearance noted by Mabrey [Ref 
2] W'as not observed. 

After checking the evacuated system for possible contamination, the test platform 
was charged with fresh R-1 13. Coolant flow was again started. No orange peel, or 
dimpled appearance could be seen on the top condenser tube. 

Power (~ 20 kW) was applied to the heaters while the coolant temperature was 
raised to ~ 0"C. An orange peel appearance was noted for the top tube. Measurements 
of the vapor temperature at the top and bottom of the condenser showed approximately 
a 20 °C differential. 

A vacuum pump, connected to the top of the condenser, w'as used to remove any 
noncondensible gases. When this was done, the temperature differential decreased to less 
than 0.5° C. The orange peel appearance quickly disappeared from the top condenser 
tube. This suggests that at least part of the contamination problem reported bv Mabrev 
W’as due to the presence of noncondensibles. 

In a subsequent series of experiments conducted on a set of titanium finned tubes, 
a foreign substance appeared to build up on the tubes during the experimental run. This 
occurred despite continuous operation of the vacuum pump. The tubes were removed 
from the condenser and inspected. The foreign material had an oily feel to it and an odor 
similar to that of ethylene glycol. 

The titanium finned tubes had smooth extensions (6" in length) attached to their 
ends using epoxy. These extensions penetrated the nylon block and steel plate and pro- 
vided the smooth surface necessary for the O-rings to give an airtight seal between the 
condenser and the atmosphere. When these tubes w'ere removed from the condenser, it 
was noted that the extensions on two of the tubes had come loose. The condenser was 
opened for inspection. Upon removal of the end plate and shroud assembly, pools of 
ethylene glycol w'ere seen inside the condenser shell and on the auxiliary’ condenser tubes. 
In addition, ethylene glycol w'as found around the shroud assembly. The entire 
condenser and evaporator assembly were therefore dismantled and cleaned in accordance 
with the procedures outlined in "Experimental Apparatus". 
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After reassembly of the condenser and evaporator, the smooth copper tubes were 
cleaned and reinstalled. Operation of the test platform failed to show any orange-peel 
or dimpled appearance on the tubes. In addition, all tubes demonstrated good filmwise 
condensation. 

It should be mentioned that in subsequent tests, copper/nickel finned tubes w'ere 
installed in the condenser. However, unlike the titanium finned tubes, the smooth ex- 
tensions on these tubes were soldered. During system operation, no foreign material 
was noted betw'een the fins. Condensate retention between the fins w'as negligible. The 
tubes demonstrated good film wdse condensation. 

In summary', it is believed that the contamination problems noted bv Mabrev [Ref 
21 were due to two problems. The first problem was the presence of noncondensibles 
during operation of the condenser. This resulted in the relatively low values for overall 
heat transfer as well as the fact that the top tube, which would be most affected bv the 
noncondensibles, had the lowest heat transfer. Secondly, while it can not be proved 
conclusively, the author believes that ethvlene glvcol mav also have been present in the 
system for certain tubes and that this is what was the material noted in between the fins 
bv .Vlabrev. 

B. OVERALL HEAT TRANSFER COEFFICIENT 

In the following sections, the results pertaining to the overall heat transfer coeffi- 
cient are presented. It is important to keep in mind that this value is experimentally de- 
termined and does not rely on any correlations other than for the fluid properties. In 
addition, it should be mentioned that in the following graphs, tube A is the top tube in 
the bundle, tube B the second, tube C the third, and tube D the bottom tube. 

1. Smooth Copper Tubes 

Figure 8 shows the overall heat transfer coefficient, (/„ plotted as a function 
of coolant velocity. The tubes did not contain any inserts. Note the low values for U,. 
These values are almost constant for velocities over the range of 0.2 to 0.7 m/s. In this 
series of experiments, Reynolds numbers were typically of the order of 300 to 2300 in- 
dicating that flow is laminar. Without any insert to promote mixing, one would there- 
fore expect relatively low values for . This is because, in laminar flow, the inside 
resistance to heat transfer is relatively large and therefore leads to a lower value for heat 
transfer. Finally, there is essentially no difference between the tubes at low coolant ve- 
locities. This suggests that with these relatively low' condensation rates, inundation is not 
a significant factor. At higher flow rate, it can be seen that the values for heat transfer 
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start increasing and that difTerences between the tubes become more apparent. As 
coolant velocity is increasing, flow is beginning to become turbulent. This reduces the 
inside resistance to heat transfer which leads to an increase in the overall heat transfer. 
As heat transfer increases, the effects of inundation become more apparent. This is pre- 
sumably due to the outside resistance becoming more important as the inside resistance 
becomes less important. 

The results obtained when these smooth copper tubes were tested individually 
is shown in Figure 9. Here, all the tubes behave similarly. In addition, comparison with 
Figure 8 shows that the values of for the tubes when tested individually are approx- 
imately the same as for the top tube during the bundle run as one would expect. 

The results for the same run with tape inserts installed are shown in Figure 10. 
It can be seen that U, increases as coolant velocity increases. The values are clearly 
greater than when no insen is used (approximately 50% greater) indicating the advan- 
tage gained when enhancing the inside heat transfer coefficient in situations where 
coolant flow is laminar. The effect of condensate inundation on the overall heat transfer 
coefficient is also clearly seen. Here, condensate inundation results in approximately a 
15% decrease in when the top tube is compared to the second tube. Interestingly, 
there appears to be little difference in the values for for tubes B, C, and D indicating 
that the condensate film is of similar thickness for these tubes. The data for a similar 
run performed by Mabrey [Ref. 2] is shown in Figure 11. The data obtained in the 
present study is almost twice the values reported by Mabrey. In addition, the present 
data shows that the top tube has the highest heat transfer which is to be expected as- 
suming noncondensibles are not present. This was not the case in Mabrey's work. 

The results for these same tubes when operated individually are shown in 
Figure 12. All the curves are essentially the same with values for U, corresponding to 
that obtained for the top tube during the bundle run. The maximum variation is about 
12 %. 

The next series of experiments were performed on the smooth copper tubes with 
the HEATEX elements installed. Figure 13 shows the results of these experiments when 
the tubes were operated as a bundle. As with the experiments performed with the 
twisted tape inserts, there is an increase in as coolant velocity is increased. In this 
case the increase is 100% in comparison to the case with no insert! Also, as in the pre- 
vious experiment, condensate inundation appears to result in about a 15% decrease in 
U„ for the second tube compared to the top tube. There appears to be little difference 
in the heat transfer coefficients for tubes B, C, and D, although some effect may be ar- 
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Figure 8. Bundle run on smooth copper tubes without inserts. 
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Figure 9. Single tube run on smooth copper tubes without inserts. 
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Figure 10. Smooth copper tubes with tape insert operated as a bundle. 
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Figure 12. Smooth copper tubes with tape insert operated individually. 
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gued at the higher coolant velocities. This may reflect the fact that, at the higher coolant 
velocities, the inside resistance decreases while the outside resistance increases in im- 
portance. Anything which causes further increases in this resistance, such as condensate 
inundation, may cause noticeable effects in heat transfer. 

As before, these same tubes with the HEATEX inserts installed were tested in- 
dividually; the results are shown in Figure 14. The values of U, for all four tubes are 
virtually identical and match the values for the top tube obtained during the bundle run 
(Figure 13). The values appear to be about 10% greater than for the tape insert ex- 
periments indicating that the HEATEX provides better tube-side enhancement. 

While the preceding figures have shown the effect of inserts and also inundation 
on the heat transfer coefficient, two important questions must be addressed. First, what 
is the uncertainty associated with the calculated values for V\. Secondly, the repeat- 
ability of the data must also be examined. 

An uncertainty analysis was performed on U, as shoAvn in Appendix F based 
on the procedures of Kline and McClintock [Ref 54]. At the low coolant velocities, the 
uncertainty w'as 5.2%. This uncertainty was based primarily on the uncertainty in the 
mass flow measurement. Here, one must include a scale interpolation term as well as a 
term to account for the time-wise jitter for the flowmeter. Since the flow rate is low for 
this coolant velocity, these two terms become more significant than at the higher flow 
rates. This is especially true if the magnitude of the time-wise jitter is the same for all 
flow rates which was the case in this study. 

At high coolant velocities, the uncertainty in the value for (/„ was about 5.6%. 
Here, the major factor driving this value is the uncertainty in the determination of tem- 
perature. The reason for this is that the AT between the tube inlet and outlet is smaller 
than for low coolant velocities. Since the relative uncertainty for the thermocouples re- 
mains the same regardless of coolant velocities, uncertainty in the temperature meas- 
urement, specifically the LMTD, becomes the driving factor in determining the 
uncertainty in 

The uncertainty analysis performed here was only for the overall heat transfer 
coefficient. U„. In this analysis, any uncertainty in the correlations used to determine the 
physical properties of the ethvlene glvcol coolant or R-113 were neglected. Therefore, 
the level of uncertainty reported here is slightly conservative. 

In order to examine the reproducibility of the data, repeat experiments were 
performed on the smooth copper tubes with the HEATEX elements installed. It is im- 
portant to mention that during the time between the two runs, the tubes were removed 
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Figure 13. Smooth copper tubes with HEATEX inserts. Bundle operation. 
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Figure 14. Smooth copper tubes with HEATEX installed. Individual operation. 
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from the condenser, cleaned and reinstalled. The results for the duplicate experiment are 
shown in Figure 15. When these values are compared to Figure 13, differences of the 
order of 3% is seen in the values of (/„ for the top tube. for the lower tubes are slightly 
larger with a more pronounced difference between the tubes. However, these differences 
all lie within the 6% uncertainty band indicating that the system gives good 
reproducibility. 

In a similar manner, the results for a repeat experiment on the smooth copper 
tubes tested individually w'ith the HEATEX installed is shown in Figure 16. When the 
values for showm here are compared with those showm in Figure 14, no significant 
differences are seen. Again, this demonstrates the good reproducibility obtained by this 
system. However, a word of caution is in order. In all the experiments reported here, 
great care has been taken to control as many variables as possible. Of particular im- 
portance is to ensure that the heater power levels are set at the same level for each ex- 
periment. This will tend to keep vapor velocity in the condenser at similar levels for each 
experiment which will result in better repeatability. 

The effects of having no insert as opposed to the twisted tape or HEATEX in- 
sert is shown in Figure 17. In this figure, the values for U, for the top tube only have 
been plotted as a function of coolant velocity for the smooth copper tubes. As previously 
noted, without any insert, coolant flow is laminar and O’ is relatively constant up to 
about 0.7 m/s. When either the twisted tape or HEATEX inserts are installed, U, in- 
creases markedly as coolant velocity increases. Interestingly, the HEATEX results in 
about a 10 to 14% increase in Uq when compared to the twisted tape inserts. Hence, 
from a heat transfer standpoint, the HEATEX element is the better insert. This is 
probably because the wire loops which make up the insert, touch much of the inside 
surface of the tube thereby "scrapping" the thermal boundary layer from the wall. With 
the twisted tape insert, mixing is provided only by inducing swirling w'hich is probably 
not as an effective mechanism. 

2. Copper/Nickel KORODENSE Tubes 

The results obtained for the overall heat transfer coefficient, t/„, for 
KORODENSE tubes operated as a bundle are shown in Figure 18. This series of ex- 
periments was performed without any inserts installed. Note the similarity between these 
results and those obtained under similar conditions with the smooth copper tubes (Fig- 
ure 8). Again, values of appear relatively constant for coolant velocities up to 0.7 
m,'s. As coolant velocity increases further, values of increase dramatically. In addi- 
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tion, an efTect of condensate inundation begins to be apparent at these higher coolant 
velocities. 

Results obtained when these same tubes were operated individually are shown 
in Figure 19. Again, when there is no inundation, all tubes have similar behavior. In 
addition, the values obtained here match those obtained for the top tube during the 
bundle experiment (see Figure 18). Interestingly, the relative magnitude of Ij\ is re- 
markably similar to that obtained for the smooth copper tubes (Figure 9). This is 
somewhat surprising. The KORODENSE tubes are corrugated on both the inside as well 
as the outside. Therefore, one would have expected somewhat higher values for heat 
transfer with the KORODENSE than with the smooth copper tube. 

Values of (/, for the KORODENSE tubes with the twisted tape inserts installed 
are shown in Figure 20. This data was obtained during bundle operation. As coolant 
velocity increased, increased steadily. This was true for all four tubes. The effect of 
inundation reduced by about 15% (tube B relative to tube A). There appeared to be 
more separation in the values of for tubes B, C, and D than was the case for the 
smooth copper tubes (see Figure 10) suggesting a greater inundation effect. In compar- 
ison with the data obtained for the smooth copper tubes, the KORODENSE tubes re- 
sulted in a slight increase (about 6 to 8%) in This difference may be due to the fact 
that the insert reduces the inside heat transfer resistance so that the accuracy of the 
measurement is improved. 

The results obtained for these same tubes operated individually are shown in 
Figure 21. As expected, all four tubes behave similarly when no inundation is present. 
In addition, the values for U, match those obtained for the top tube during the bundle 
operation (Figure 20). In comparison to the smooth copper tube data (Figure 12), the 
data for the KORODENSE tubes appears somewhat higher (about 6 to 8%). However, 
given the known uncertainty, this difference is probably not significant. 

Results obtained for the KORODENSE tubes during bundle operation with the 
HEATEX elements installed are shown in Figure 22. Note that U, again increases 
markedly for all four tube as coolant velocity increases. An inundation effect can be seen 
for tubes B, C, and D which is similar to that seen in the twisted tape runs. Tube B 
shows approximately a 15% decrease in compared to tube A. This effect is magnified 
at higher coolant velocities (condensation rates) as one would expect. When this data is 
compared to that obtained for the smooth copper tubes with HEATEX installed (see 
Figure 13), it can be seen that is greater for the KORODENSE tubes. This suggests 
that roping of the tube provides some enhancement for the outside (and possibly the 
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Figure 18 . KORODENSE tubes with no insert, bundle operation. 
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61 



inside) heat transfer coefficients. While we expect some small difierence, it must be 
pointed out that the differences observed here are at the bounds of the calculated un- 
certainty. 

Results for these tubes tested individually are shown in Figure 23. As before, 
the values of U, are similar for all four tubes and match the values obtained for the top 
tube during the bundle experiments (see Figure 22). In comparison to the data obtained 
for the smooth copper tubes, the values of are somewhat higher as expected, but are 
still at the bounds of the calculated uncertainty. 

The effect of using no insert, twisted tape and HEATEX on U, can be seen in 
Figure 24. The data are for the top tube only. Note that, as with the smooth copper 
tubes (Figure 17), the HEATEX results in about a 10 to 14% increase in U, in com- 
parison to the twisted tape insert. Both provide significant enhancement in U, when 
compared to no insert (almost 100% enhancement). As with the smooth copper tubes, 
when no insert is used, U„ appears relatively constant for coolant velocities over the 
range of 0.2 to 0.7 m/s indicating laminar flow. 

3. Wire Wrapped KORODENSE Tubes 

This series of experiments was performed to determine an optimum wire diam- 
eter as discussed in "Experimental Procedure". All tubes were tested individually and all 
had HEATEX mixing elements installed. Figure 25 shows U, as a function of coolant 
velocity for the case of an ordinary KORODENSE tube together with tubes wrapped 
with 0.029", 0.049" and 0.0675" diameter wire, each tested in the top position (tube A). 
The 0.049" diameter wire gives the biggest enhancement in U, (about 25%). While the 
other two wire diameters also provide enhancement in b\ (about 15%), it is below that 
of the 0.049" wire. The no wire case agrees with the previously obtained data (see 
Figure 22). 

Figure 26 shows similar data but in this case for tube B. Again, the 0.049" di- 
ameter wire provides the greatest enhancement in U, (again about 25%). 

Figure 27 shows what happens when the different diameter wires are used on 
tube C. In this case an interesting fact can be observed. Here, both the 0.049" and 0.029" 
wire provide virtually the same degree of enhancement for U„ (about 20%) with both 
being better than the 0.0675" wire. Figure 28 shows the results obtained for tube D. 
Note that while the 0.049" wire again provides the greatest degree of enhancement for 
b\ (about 16%), the differences between the three different diameter wires are not as 
great. 
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Figure 22. KORODENSE tubes with HEATEX inserts, bundle operation. 
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Figure 23. KORODENSE tubes >vitli HEATEX inserts, individual operation. 
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Figure 26. Effect of >\ire diameter on U, for tube B. 
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Figure 27 and Figure 28 suggest that the position of the wire wrapped tube 
within the bundle may affect the overall heat transfer coefficient. In order to gain further 
insight on this, the data obtained from these experiments was replotted so that for a 
given wire diameter was plotted as a function of coolant velocity for all four tube posi- 
tions. The data for the tube without any wire is shown in Figure 29. The data is all very 
similar and no distinct trends are evident. Note that comparison of this data with 
Figure 23 again demonstrates the reproducibility of the system. 

Figure 30 (0.029" wire diameter) shows that tubes A, B and C behave similarly. 
Tube D produces consistently lower values for U, than the other tubes. Figure 31 
(0.049" wire) demonstrates what appears to be a positional effect on t/,. The data for 
tube A and B is virtually identical and lies above that of tubes C and D which are also 
virtually identical. Figure 32 suggests that there is also a positional effect for the 
0.0675" wire. The top tube seems to provide the greatest enhancement in (/,. The level 
of enhancement for tubes C and D are approximately the same, both being lower than 
tubes A and B. 

Figure 33 shows U, plotted as a function of wire diameter for an arbitrarily 
chosen coolant velocity of 1.0 m/s. Note that the strongest positional affect appears to 
be for the 0.049" wire. Also, it is clearly evident that the 0.049" wire provides the greatest 
enhancement for heat transfer. 

The possible positional effect may be due to vapor velocity. Although the' ve- 
locity of the vapor was very small in this study (approximately 0. 1 m/s) there may still 
be an effect. As the vapor penetrates the tube bundle, it slows down. As its velocity is 
reduced, it can no longer help to strip condensate from the tubes. Hence, heat transfer 
is decreased. This effect looks much like an inundation effect. An interesting point is 
why this shows up in these experiments and not wdth the plain KORODENSE or 
smooth copper tubes. With the wire wrapped tubes, the condensate layer is presumably 
thin. Therefore, anything which could help to thin it funher w’ould be of importance in 
improving heat transfer. The loss of this small amount of vapor shearing may be im- 
portant in this case. 

In this study, the 0.049" diameter wire appears to give the best enhancement in 
heat transfer (see Figure 33). This yields a pitch to wire diameter ratio of 7.2. Fujii et 
al. [Ref 33] found that for condensation of R-12 on wire wrapped smooth copper tubes, 
the optimal pitch-to-wire diameter ratio w'as on the order of 2. This difference could be 
due to our use of the KORODENSE tubes. For a pitch-to-wire diameter ratio of 2, we 
would have had to use wire w'hich was of the order of 5 mm in diameter. This would have 
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Figure 30. Effect of 0.029" wire on U, for all tubes. 
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Figure 32. Effect of 0.0675" wire on U, for all tubes. 
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been unrealistic for our experiments. It may be that the optimal pitch-to-wire diameter 
may in fact be about 2. Additional experiments need to be performed with different 
pitch KORODENSE tubes. 

4. KORODENSE Tubes Wrapped with 0.049" Wire 

Once the optimal wire diameter had been determined, two identical experiments 
were performed on the KORODENSE copper/nickel tubes wrapped with the 0.049" di- 
ameter wire. The reason that two experiments were performed was to check the repeat- 
ability of the data. Figure 34 and Figure 35 show the results from these two experiments 
for bundle operation of the tubes. The two experiments show excellent repeatability. In 
comparing these results with the non-wire wrapped tubes (Figure 22), several interesting 
points can be made. The values of for all four tubes are greater for the wire wrapped 
tubes than for the plain tubes (approximately 16% for the top tube). In addition, the 
inundation effect (seen as a lowering of L\) is reduced for the wire wrapped tubes. 
Finally, note that the differences between the tubes are more apparent than for the 
non-wire wrapped tubes. This suggests that the wire wrapping may have an effect upon 
condensate inundation. This will be addressed in more detail later. 

Figure 36 and Figure 37 show the results obtained when the tubes wrapped 
with the 0.049" diameter wire were operated individually. All the tubes behave similarly 
and demonstrate good reproducibility between the two runs. 

5, Copper/Nickel Finned Tubes 

Figure 38 shows the results obtained when no insert was installed and the tubes 
were operated as a bundle. Comparison of these results with those shown in Figure 8 
for the smooth tubes shows that the finned tubes yield higher heat transfer coefficients 
as expected (approximately 12%). Also, as with the smooth copper tubes, the top tube 
again displays the highest value of U„, due to the fact that the lower tubes suffer from 
condensate inundation. The single tube data for the finned tubes is shown m 
Figure 39. All four tubes are virtually the same as the top tube during bundle operation. 

Figure 40 and Figure 41 show the results obtained for the fumed tubes with the 
tape insert installed during bundle and individual tube operation respectively. In bundle 
operation, the values of U, are nearly double those obtained for both the smooth copper 
tubes (Figure 10) and the KORODENSE tubes (Figure 20) using similar inserts. As 
with the smooth and KORODENSE tubes, inundation again causes a reduction in . 
This reduction for the second tube appears to be about 10%. Comparison of the single 
tube operation (Figure 41) with the bundle operation (Figure 40) shows that all four 
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Figure 34. 0.049" nire wrapped KORODENSE tubes, bundle operation, e.xpt. 1. 



77 



UZERO (W/M**2 K) 

SOO.O 600.0 700.0 800.0 900.0 1000.0 1100.0 



0.049" WIRE WRAPPED KORODENSE 



KD12, R-113, HEATEX, BUNDLE 

















: 




□ ' 


p 












i 


1 ' 


: C 

9 

1 c 

6 




o 

A 1 


? 

[ 




t 


t 

1 i 


T 

j 

C 


? 

i 


t 

b 

^ j 


\ 

L .i 


) i 

i i 

H 


; 1 

I -i 

h 




4 


i- 




c 


> \ 
\ 




h . : 
























LEGEND 
□ = TUBE A 
o = TUBE B 
A = TUBE C 
+ = TUBE D 














i 


i i 


i i 






i i 


: 1 


i i 


i 


f— 


i 


i 


i i 


1 



0.1 0.2 0,3 0.4 0.0 0.6 0.7 0.6 0.9 1 1.1 1.2 1.3 

VELOCITY (M/SEC) 



Figure 35. 0.049" wire wrapped KORODENSE tube.?, bundle operation, e.xpt. 2. 
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Figure 36. 0.049" wire wrapped KORODENSE tubes, individual operation. 
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Figure 37. 0.049" wire wrapped KORODENSE tubes, individual operation. 
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Figure 38. Copper/nickel finned tubes, no insert, bundle operation. 
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Figure 39. Copper/nickel Tinned tubes, no insert, individual operation. 
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tubes operated individually have approximately the same value of (/, as the top tube 
during the bundle operation. 

Data obtained from experiments on the copper/nickel fmned tubes with the 
HEATEX elements installed are shown in Figure 42 and Figure 43 for bundle and in- 
dividual tube operation respectively. Note that in comparison with the tape insert ex- 
periments (Figure 40 and Figure 41), use of the HEATEX inserts results in about a 
20% increase in which again indicates the improvement afforded by using this type 
of insert. Again the inundation effect results in a decrement of 10 to 14%. During in- 
dividual tube operation, all four tubes again behave like the top tube during the bundle 
operation. 

As with the smooth copper tubes and the KORODENSE tubes, the effects of 
using different types of inserts can be summarized. This has been done in Figure 44. 
This data is for the top tube only. The similarity with data from the other tubes (see 
Figure 17 and Figure 24) is good with HEATEX yielding approximately a 20% increase 
in over the twisted tape insert and a remarkable 200% increase over the case with no 
insert! 

C. OUTSIDE HEAT TRANSFER COEFFICIENT 

The following figures deal with the outside heat transfer coefficient, /i„. In all cases, 
except as specifically noted, the outside heat transfer coefficient was determined using 
the modified Wilson Plot procedure. As before, tube A is the top tube with tube D being 
the bottom tube in the bundle. In addition, no values were obtained for /i. for exper- 
iments which had no insert. The reason for this is that, for the coolant velocities used, 
without inserts, the inside resistance was simply too large for any accurate results to be 
obtained for h„. 

1. Smooth Copper Tubes 

Figure 45 and Figure 46 show the values obtained for A, for the smooth copper 
tubes with the tape insert installed during bundle and individual tube operation respec- 
tively. Several features deserve comment. First, values of for the top tube are below 
those predicted by Nusselt theory'. These differences range from about 15% at low 
coolant velocities to about 7% at the higher coolant velocities. In addition, note the 
inundation effect (seen as a reduction in /i„ for tubes B, C, and D). Reductions in range 
from 28% to 19% as coolant velocity increases from 0.2 to 1.2 m/s respectively. As was 
the case for the overall heat transfer coefTicient, when the tubes were operated individ- 
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Figure 40, Copper/nickel finned tubes, tape insert, bundle operation. 
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Figure 41. Copper/nickel finned tubes, tape insert, individual operation. 
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Figure 42. Copper/nickel finned tubes, HEATLX insert, bundle operation. 
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Figure 43. Copper/nickel ruined tubes, HEATEX insert, individual operation. 
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Figure 44. Summary’ of effect of tube inserts on U, 
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ually, all had values for /i„ which were virtually identical to those obtained for the top 
tube during bundle operation. 

Figure 47 and Figure 48 show values for /i„ obtained when the HEATEX ele- 
ments were installed in the copper tubes during bundle and individual tube operation 
respectively. As was the case for the tape insert, values of range from about 10 to 
14% below that predicted from Nusselt theory. Again, there is a substantial inundation 
effect which results in a 15 to 18% decrease in h, when comparing the top tube to the 
second tube. Note that the values for h, obtained for the individual tubes agree with 
those obtained for the top tube during bundle operation. Finally, comparison of the 
value for /?„ obtained when the tape inserts were used as opposed to the HEATEX ele- 
ments reveal no effect on h, which is to be expected since we are examining the outside 
heat transfer coefficient. 

It is interesting to note that Goto et al. [Ref 5] found that values of /i„ for 
R-113 condensing on long, smooth copper tubes were about 10% below that predicted 
by Nusselt theory. Marto et al. [Ref 8 ] found that the vapor-side heat transfer coeffi- 
cient for R-113 condensing on short (0.133m) smooth copper tubes was in relatively 
good agreement with Nusselt theor>'. This result was subsequently confirmed by Michael 
et al. [Ref 32]. However, it is important to point out that both the latter reports used 
very short condenser tubes. The tubes used by Goto et al. were about 0.6 m in length, 
while the tubes used in this study had a condensing length of 1.2 m. Hence tube length 
may be an important consideration. The reason for this is as follows. The coolant en- 
ters the condenser tube with a given inlet temperature. As the coolant flows through the 
tube, the coolant temperature increases. The wall termperature is also changing down 
the length of the tube. At the tube entrance, the temperature difference between the tube 
wal and the vapor saturation temperature causes the formation of a thick condensate 
layer. This layer is thinner at the tube outlet due to the smaller wall-vapor saturation 
temperature difference. This causes condensate to flow in the axial direction. These fac- 
tors taken together imply that the temperature difference is a function of the length of 
the tube. This factor is not taken into account in the Nusselt analysis. As coolant ve- 
locity increase^; the difference betwee the coolant inlet and outlet temperatures is re- 
duced. The wall temperature becomes more constant along the tube length so that values 
for /j„ are in better agreement with those detennined experimentally. 

The effect of using the twisted tape and HEATEX insert on the inside and out- 
side heat transfer coefficients is demonstrated in Figure 49 and Figure 50 respectively. 
In both cases, the values for are virtually identical as expected. However, the 
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Figure 45. Smooth Copper tubes, tape insert, bundle operation. 
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Figure 46. Smooth Copper tubes, tape insert, individual operation. 
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Figure 47. Smooth Copper tubes, HEATEX insert, bundle operation. 
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Figure 48. Smooth Copper tubes, HEATEX insert, individual operation. 
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HEATEX elements essentially double the inside heat transfer coeflicient. This is what 
provides the improvement in overall heat transfer coefficient when the HEATEX ele- 
ments are used. 

The effects of inundation were examined by calculating an average and a local 
heat transfer coefficient. The local h, is found by dividing the value for h, for a given 
tube by the value for the top tube. Values obtained are then compared with Nussclt 
theory (Eq. 2-18) and the Kern model (Eq. 2-20). The average heat transfer coeflicient 
is determined by taking an average value of for the tube of interest plus the tubes 
above it and dividing this value by h, for the top tube. Values obtained are then com- 
pared to Nusselt theory (Eq. 2-17), the Kern model (Eq. 2-19), or the Eissenberg model 
(Eq. 2-12). For both the local and average values, the flatter the curves, the less the ef- 
fect of condensate inundation. 

The effects of inundation on /i, for the smooth copper tubes are shown in 
Figure 51, Figure 52, Figure 53 and Figure 54 at three arbitrarily chosen coolant ve- 
locities. As expected, there is virtually no difference on the local or average values of the 
heat transfer coefficient when the tape insert is used as opposed to the HEATEX ele- 
ment. In addition, the values for the average heat transfer coefficient lie in between the 
values predicted by the Kern and Eissenberg models (equations 2-19, and 2-21). Values 
for the local heat transfer coefficient lie slightly above those predicted by the Kern model 
(equation 2-20). These results suggest that Nusselt theory is too conservative with re- 
gard to the effect of condensate inundation. It is interesting to note that the main 
inundation effect occurs for the second tube with little additional effect for subsequent 
tubes in the bundle. This again supports the notion that the smooth copper tubes can 
only support a limited amount of condensate 
2. KORODENSE Tubes 

Figure 55 and Figure 56 show values for h, for the KORODENSE tubes with 
the tape inserts installed. The figures show data obtained when the tubes were operated 
as a bundle and individually respectively. The values of h, for the top tube are above that 
predicted from Nusselt theory which is due to the small enhancement provided by the 
roping of the tubes. This enhancement is approximately 15 to 20% in comparison to the 
smooth copper tubes. As with the smooth copper tubes, there is again a substantial ef- 
fect of inundation (on the order of 20% for the second tube compared to the top tube). 
There is also more separation between the values of for tubes B, C, and D than for 
the smooth copper tubes indicating that the lower tubes may not be "saturated" with 
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condensate. Finally, note that the relative values for for the KORODEN'SE tubes are 
about 10 to 15% higher than the corresponding values for the smooth copper tubes. 

The data for the KORODEN'SE tubes operated individually show much more 
scatter than was seen for the smooth copper tubes especially at low coolant flow rates. 
Also, rather than being virtually identical with those obtained for the top tube during 
bundle operation, the values of lay scattered between those for the top tube (bundle 
operation) and the values predicted by Nusselt theory. Presumably, this is due to either 
inaccuracies or uncertainty in the inside correlation especially at low coolant flow rates. 
At higher coolant velocities, the scatter appears to be within about 10%. 

Values for for the KORODEN’SE tubes obtained when the HEATEX ele- 
ments were installed are shown in Figure 57 and Figure 58 for bundle and individual 
tube operation respectively. As wth the experiments conducted with the tape inserts 
installed, values of /i, again lie above those predicted by Nusselt theory. The enhance- 
ment appears to be approximately 15 to 20%. Again there is a substantial inundation 
effect which is of the order of 20% for the second tube in comparison to the top tube. 
The data obtained for the individual tube operation shows less scatter. However, values 
for all the tubes again lie between those predicted by Nusselt theory and those obtained 
for the top tube during bundle operation. 

As with the smooth copper tubes, the effects of changing inserts on the inside 
and outside heat transfer coefficients are shown in Figure 59 and Figure 60. The inside 
heat transfer coefficients are identical to those of the smooth copper tubes (Figure 49 
and Figure 50). It was felt that the enhancement provided by the roping would be much 
less than the enhancement provided by the inserts. Therefore, the same correlations 
were used to determine the inside heat transfer coefficient. The values for h, are slightly 
greater for the KORODEN'SE tubes than for the smooth copper tubes (approximately 
5 to 10%). 

Figure 61 and Figure 62 show the local heat transfer coefficient for the four 
tubes for the tape insert and the HEATEX elements. There is virtually no difference 
between the two. The values lie above the Kern estimate. Figure 63 and Figure 64 show 
the average heat transfer coefficient for the four tubes with both inserts. The differences 
for values using the tape inserts as opposed to the HEATEX elements are again negli- 
gible. The average values for h, also lie just above the Kern model. This is in contrast to 
the smooth copper tubes where the average values lay near the Nusselt estimate. This 
suggests that the effect of condensate inundation is indeed less for the KORODEN’SE 
tubes. 
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Figure 55. KORODENSE tubes, tape insert, bundle operation. 
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Figure 56. KORODENSE tubes, tape insert, individual operation. 
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Figure 57. KORODENSE tubes, HEATEX insert, bundle operation. 
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Figure 58. KORODENSE tubes, HEATEX insert, individual operation. 
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Figure 60. Effect of HEATEX insert on /r, and A, 
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Figure 61 . Local h, for tape insert. 
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Figure 62. Local //, for HEATEX element. 
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3. Wire Wrapped KORODENSE Tubes 

The efiect on of WTapping wire around the KORODENSE tubes is shown in 
Figure 65, Figure 66, Figure 67 and Figure 68. In viewing these figures, one should 
remember that the values for the unwrapped KORODENSE tubes were about 10% 
greater than predicted from Nusselt theory- and about 15% greater than for the smooth 
copper tubes. From the figures, it appears that the 0.049" diameter wire results in the 
greatest enhancement in for all four tube positions. However, the relative magnitude 
of the enhancement decreases as lower tubes are examined. Since this data was obtained 
during individual tube operation, the results are not due to condensate inundation. 

Values for are shown for each wire diameter in all four tube positions in Fig- 
ure 69, Figure 70, Figure 71 and Figure 72. When no wire is used, the results are the 
same as those shown in Figure 58. For the 0.029" wire, positions A and B are essentially 
identical with tube C showing the greatest enhancement and tube D the worst. The 
magnitude of the deviation is about 20%. The 0.049" diameter wire results in the best 
enhancement in /j„ regardless of position. However, note that position D is not over- 
whelmingly greater than for the 0.029" wire. The data show that the level of enhance- 
ment depends upon bundle position. The top tube positions behave similarly while 
positions C and D show reduced levels of enhancement Finally, the 0.0675" diameter 
wire yields values for which are nearly the same although a positional effect appears 
evident. An interesting comparison are the values for h, for the 0.049" and 0.0675" wire 
in the fourth position. Here it appears that the 0.0675" wire gives roughly the same level 
of enhancement especially at high coolant velocities. 

4. KORODENSE Tubes Wrapped with 0.049" Wire 

The effects on of wrapping 0.049" diameter wire on KORODENSE tubes on 
are shown in Figure 73 and Figure 74. These two experiments are identical runs. In 
both cases, there is a substantial improvement in for all tubes when compared to val- 
ues predicted from Nusselt theory’. It is interesting to note that even though tubes B, 
C, and D are all being inundated, the values for h, are still well above Nusselt theory 
(32% for tube A, 28% for tube B, 22% for tube C and 16% for tube D). In comparison 
with the non-wire wTapped KORODENSE tubes, wrapping with the 0.049" diameter 
wire yields an increase of about 70% for tubes A and B and about a 30 to 50% increase 
for tubes C and D. Comparison w-ith the top smooth tube yields almost a 90% increase 
in h,. 

The results for the single tube runs are shown in Figure 75 and Figure 76. Both 
runs demonstrate considerable scatter in the data of about + 22%. It is also interesting 
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Figure 65. Effect of wire diameter on A, for tube A. 
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Figure 66. Effect of >\iie diameter on h, for tube B. 
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Figure 67. Effect of >>ire diameter on for tube C. 



115 



HZERO (W/M**2 K) 

800.0 1000.0 IZOO.O 1400.0 IGOO.O IBOO.O ZOOO.O 



HZERO WIRE WRAPPED KORODExNSE 

TUBE D, R-113, HEATEX 



P □ 



LEGEND 
□ = NO TVIRE 
o = 0.029" WIRE 
A = 0.049" WIRE 
+ = 0.0675" WIRE 



+ 

o 



o 



+ 

o 



—I ! 1 ^ 1 

0.5 0.6 0.7 0.8 0.9 

VELOCITY (M/SEC) 



0.1 0.3 0.3 0.4 



1.0 1.1 1.2 1.3 



Figure 68. Effect of wire diameter on for tube D. 
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Figure 69 . Effect of no >\ire on /i, for all tubes. 
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Figure 70. Effect of 0.029" wire on h„ for all tubes. 
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Figure 72. Effect of 0.0675" wire on for all tubes. 
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Figure 73. 0.049" wire wrapped KORODENSE tubes, bundle operation, expt. 1. 
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Figure 74. 0.049" wire wrapped KORODENSE tubes, bundle operation, expt. 2. 
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to note that the data for these single tube runs is somewhat higher than for the bundle 
runs (of the order of 11%). This may be due to uncertainty in the correlations used for 
/i, which are reflected in 

The values for the average and local vapor-side heat transfer coefficients are 
shown in Figure 77, Figure 78, Figure 79 and Figure 80 for the two experiments con- 
ducted with the 0.049" wire wrapped KORODENSE tubes. The average vapor-side heat 
transfer coefficient is equal to that predicted by the Eissenberg model whereas the local 
vapor-side heat transfer coefficient is substantially above the Kem estimate. Recall that 
for the plain KORODENSE tubes, the average heat transfer coefficient was above the 
Kem estimate but below that predicted by the Eissenberg model. The local heat transfer 
coefficient was just above the Kern estimate. The differences between the plain and 
wire-wrapped KORODENSE tubes demonstrate that wrapping the tubes with wire 
substantially reduces the effects of condensate inundation. 

5. Copper/Nickel Finned Tubes 

The values for obtained from experiments on copper/nickel finned tubes using 
the tape inserts are shown in Figure 81. The values obtained for the top tube are ap- 
proximately 7 times greater than predicted by Nusselt theory. Condensate inundation 
caused about a 35% decrease in h, for the second tube. 

The values for obtained using the same tubes but with the HEATEX elements 
are shown in Figure 82. In comparison to the values obtained for the tape inserts fsee 
Figure 81), this data appears very scattered with enliancement over Nusselt theory of 
the order of 10 to 15 fold. This is a large discrepancy. It is possible that the correlation 
for the small diameter HEATEX element may be suspect. The data supplied by CAL 
GAVIN was collected on smooth copper tubes with a diameter similar to that of our 
copper/nickel finned tubes. This data was collected over a lower coolant velocity range 
than that used in this study. Therefore, the correlation for this HEATEX element may 
not be valid for the finned tube over the range of coolant velocities used here. This area 
needs to be studied further, perhaps even to gathering new experimental data which can 
be used to develop a new correlation. 

The degree of inundation can be seen in Figure 83, Figure 84, Figure 85 and 
Figure 86. The data for the average vapor-side heat transfer coefficient for both 
HEATEX as well as tape inserts lie between the Nusselt and Kern models. This is also 
true for the local vapor-side heat transfer coefficients. 

An interesting comparison is to look at the average and local heat transfer val- 
ues for all tubes tested. Regarding the average heat transfer coefficient, for the smooth 
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Figure 75. 0.049"' wire wrapped KORODENSE, individual runs, expt. 1. 
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Figure 76. 0.049* wire wrapped KORODENSE, individual runs, expt. 2. 
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Figure 77. Average h, wire wrapped KORODENSE, expt. 1. 
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Figure 78. Average h, wire wrapped KORODENSE, expt. 2. 
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Figure 79. Local h, wire wrapped KORODENSE, expt. 1. 
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Figure 80. Local h, wire wrapped KORODENSE, expt. 2. 
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Figure 81. A, for copper/nickel finned tubes, tape insert. 
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Figure 82. A, for copper/nickel finned tubes, HEATEX insert. 



131 



copper tubes, the values lay between the Kern and Eissenberg estimates. The plain 
KORODENSE tubes had values similar to those for the smooth copper tubes w’hile the 
0.049“ wire-wrapped KORODENSE had values which matched estimates based on 
Eissenberg's model. The copper/nickel finned tubes, on the other hand, had values which 
were between the Kern and Nusselt estimates. Since the lower the values the more 
condensate effect present, this data would imply that the effect of inundation is worse 
on the finned tubes and least on the wire-wrapped tubes. If the data for the local heat 
transfer coefficients is analyzed in a similar manner, the same conclusions are reached. 

The wire wrapped KORODENSE tubes have the least effect of condensate 
inundation presumably because the ttire serves to create a low pressure area at the base 
of the wire. This causes condensate to be pulled along the tube which enhances heat 
transfer. The values for the finned tubes are not in agreement with those of Honda et 
al. [Ref 42]. These investigators found that the degree of inudation effect was substan- 
tially less than that reported here. There are two apparent reasons for this. The pitch- 
to-diameter ratio used by Honda et al. was 1.37 whereas in the present study it is 2.25. 
This difference may cause differences in the way condensate drains from one tube onto 
another. Alternatively, the correlations used to determine the heat transfer coefficient 
may be suspect. Hong and Bergles [Ref 47] have reported that their correlation for the 
twisted tape has an uncertainty of about 15%. Therefore, it is possible that the large 
effect of condensate inundation is due more to uncertainty in the correlations than to a 
real phenomenon. 

In order to assess the validity of the data obtained on the copper/nickel finned 
tubes, a comparison was made with the model of Beatty and Katz [Ref 13]. In order to 
do this comparison, the fin efficiency had to be determined using equations from 
Incropera and DeWitt [Ref 55]. It was assumed that the fins could be approximated 
as rectangular fins with adiabatic tips. The fin efficiency was then approximated by: 



According to Incropera and DeWitt, the term niL^ can be approximated by [Ref 55]: 



tanh(mL^) 



(6-1) 







(6-2) 
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Figure 83. Average //, for Cu/is'i finned tubes, tape insert. 
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Figure 84. Average h, for Cu/Ni finned tubes, HEATEX insert. 
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Figure 85. Local h, for Cu/Ni finned tubes, tape insert. 
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Figure 86. Local for Cu/Ni finned tubes, HEATEX insert. 
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where is defined as; 



L,= L + ^ (6-3) 

Here, L is defined as the height of the fin and t is the fin thickness. The term A,, is de- 
fined as: 



Ap = L,t (6 - 4) 

In order to calculate a value for the fin efficiency, an iterative procedure was 
used. First, a value for /?„ was guessed. Based on this value, i/y, was calculated using 
equation (6-1). An effective area was then calculated from the following equation (Marto 
(Ref 3]): 



Aef- V/Afs + rifAfi + A 



(6-5) 



where A^ is the surface area of the fins sides, Ap, is the surface area of the fin tips and 
A„ is the interfin surface area. Using these values, the effective diameter, D, was deter- 
mined from (Marto [Ref 3]): 




1.30?7/- 
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where L is defined as: 



( 6 - 6 ) 



n{Dl - d1) 

4Z), 



(6-7) 



The heat transfer coefficient is then calculated from a modified version of the Beatty and 
Katz model (Eq. 2-2) which uses heat flux rather then the film temperature difference, 
AT. 



= .6085 




(6-8) 



The results of these calculations for the copper/nickel finned tubes with the tape 
inserts are shown in Figure 87. In this figure, the values of A. predicted from the Beatty 
and Katz model have been expressed in terms of the smooth tube area, that is ndX . 
Also note that values for /j„ have been calculated for copper, copper/nickel and titanium. 
It is important to note that the raw data lies below the predicted values for copper/nickel 
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finned tubes. The data obtained after reprocessing with the modified Wilson plot proce- 
dure lies somewhat above the prediction. However, both lie within + 20% which is 
within the range of validity for the Beatty and Katz model. Therefore, it would appear 
that the data obtained in the present study for the copper/nickel finned tubes with the 
tape inserts are reasonable. 

The results for the copper/nickel finned tubes with the HEATEX elements in- 
stalled are not as good. These results are shown in Figure 88. This figure shows the 
tremendous amount of scatter in the data irrespective of whether the raw data or the 
data obtained from the modified Wilson plot procedure is used. In either case, it is ap- 
parent that the data do not match values predicted from the model of Beatty and Katz. 
The only difference in these two experiments is the type of insert used. The values for 
are based on an experimentally determined value for U, and a value for /i, which is 
obtained from a correlation. In this study, the correlations for the HEATEX elements 
were based on empirical data supplied from the manufacturer. This data was obtained 
on smooth copper tubes having the same diameter as the copper/nickel firmed tubes used 
in this study. Another point which must be considered is that the range of velocity for 
which the empirical data was supplied was from 0.2 to 1.2 m/s. However, for the 
copper/nickel finned tubes used in this study, a coolant velocity of 0.7 to 2 m/s was used. 
These factors would suggest that the correlation used for the small diameter HEATEX 
element is in error. Further research will have to verify whether or not this is the case. 

D. IMPLICATIONS OF PRESENT WORK 

This study used R-113 as the refrigerant and a refrigerated ethylene glycol/water 
mixture as the coolant. However, Naval refrigeration systems use R-114 for the 
refrigerant and sea water as the condensing medium. Therefore, one must ask whether 
the results obtained from the present work can be extended to Naval apphcations. In 
order to perform the following analysis, the dimensions of the smooth copper tubes and 
copper/nickel finned tubes (26 fpi) used in this study were assumed. The inside heat 
transfer coefficient was calculated using the Dittus-Boelter correlation. For the smooth 
tubes, /i„ was calculated using the Nusselt equation, modified to use heat flux rather than 
the film temperature difference, AT . For the finned tube, /i, was again calculated using 
the Dittus-Boelter correlation. However, was calculated using the modified Beatty and 
Katz equation (Eq. 6-8). The fin efficiency was assumed to be 1.0. The condenser design 
parameters were modeled after the requirements for the DDG-51. Hence, the refriger- 
ation load was fixed at 242.7 tons. The number of condenser tubes was fixed at 340 
tubes. Coolant inlet temperature was assumed to be 65 F (18.3 C) while the saturation 
temperature of R-114 was 107.8 F (42.1 C). The coolant velocity was held constant at 
2 m, s. 
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Figure 87. Experimental data (tape insert) vs. Beatty/Katz model. 
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Figure 88. Experimental data (HEATEX insert) vs. Beatty/Katz model. 
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Ratios of the surface areas for finned to smooth tube and a finned KORODENSE 
tube to smooth tube were calculated. In order to carry out this calculation for the 
KORODENSE tube, some type of correlation must be used to calculate h,. However, 
as shown by Cunningham and Ben Boudinar [Ref 56] in their work on steam 
condensation on roped tubes, the enhancement due to the roping depends upon the type 
of groove cut into the tube. In this analysis, an enhancement of 1.6 for h, was assumed 
for the KORODENSE tube over the finned tube. Once the area ratio is obtained, a 
condenser tube length can be calculated. This is done as follows. The area ratio can be 
expressed in terms of length by: 

A„ Npid„L, 

where is the area of the finned tube, is the area for the smooth tube, and N is the 
number of tubes. The term Uj is given by: 

aj={Aj^A,)N ( 6 - 10 ) 



where A^ and A„ are the finned and unfinned areas respectively. 

Table 3 shows the results of this calculation. For the finned tube, the area ratio of 
1.75 leads to a decrease in the condenser length of about 50%. On the other hand, if a 
roped KORODENSE tube could be manufactured which had fins on the outside surface, 
the reduction in length of the condenser tubes would be 58%. This could lead to even 
greater weight savings for refrigeration systems. 



Table 3. SUMMARY OF AREA RATIO CALCULATIONS 





Smooth 

tube 


Finned 

Tube 


KORODENSE + Fins 




1.0 


1.75 


1.51 












1.0 


0A91 


0.42 


L, 
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VII. CONCLUSIONS 



The following conclusions can be made based on the results obtained in this study. 

A. GENERAL 

1. The condenser/evaporator test platform is fully operational and shows no signs of 
the contamination problems reported by Mabrey [Ref. 2 ]. 

2. Values for outside heat transfer coefficient obtained using smooth copper tubes 
appear to be reasonable and agree with other published data. 

B. OVERALL HEAT TRANSFER COEFFICIENT 

1. KORODENSE tubes provide enhancements of 6 to 8% over smooth copper tubes 
regardless of the type of insert used. 

2. For the KORODENSE tubes used in this study, wire wrapping results in an en- 
hancement in heat transfer. 

3. The optimal pitch-to-w’ire diameter was about 7 for the KORODENSE tubes used 
in the present work. 

4. Wrapping the KORODENSE tubes with 0.049" diameter wire resulted in an en- 
hancement of 16% over the non wire- wrapped KORODENSE tubes and a 22 to 
24% enhancement over the smooth copper tubes. 

5. Copper/nickel finned tubes (26 fpi) yielded a doubling of the heat transfer in com- 
parison to smooth copper tubes. 

6. The use of HEATEX radial mixing elements yielded a 10 to 15% increase in com- 
parison to the use of twisted tape inserts. 

C. OUTSIDE HEAT TRANSFER COEFFICIENT 

1. The smooth tubes yielded values of /i„ which were 7 to 15% below those estimated 
from Nusselt theory. 

2. The KORODENSE tubes gave values of which were about 15% greater than 
those estimated from Nusselt theory and approximately 25% greater than values 
obtained for the smooth copper tubes. 

3. Wrapping the KORODENSE tube with 0.049" diameter wire resulted in a 90% 
increase in the outside heat transfer coefficient in comparison to the smooth copper 
tubes. 

4. The copper/nickel finned tubes yielded a 7 fold increase in over the smooth 
copper tubes. 

D. EFFECTS OF CONDENSATE INUNDATION 

I. For smooth copper tubes, condensate inundation resulted is a 15 to 18% decrease 
in for the second tube compared to the top tube. 
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2. KORODENSE tubes showed a 20% decrease in /i„ for this same comparison. 

3. Wrapping the KORODENSE tubes with 0.049" diameter wire resulted in a decrease 
of only 5 to 10% in for the second tube compared to the top tube. 

4. For the copper, 'nickel finned tubes, condensate inundation resulted in almost a 
35% decrease in for the second tube relative to the top tube. This large decrease 
is presumably due to uncertainty in the correlations used to determine the heat 
transfer coefficient. 
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VIII. RECOMMENDATIONS 



Based on the results of this study, the following recommendations are made: 

A. GENERAL RECOMMENDATIONS 

1. The ball valves which control coolant flow to the flo\vmeters should be replaced 
with needle valves. This will provide better control of coolant flow. 

2. The flowmeters should be replaced with larger capacity meters. This will allow in- 
vestigation of heat transfer at higher coolant velocities and therefore higher 
condensation rates. This change will have to be accompanied by an increase in the 
diameter of the coolant return header in order to accommodate the increased 
coolant flow. 

3. Extra lines from the top ends of the condenser to the cold trap need to be added. 
This will ensure that no isolated pockets of noncondesibles exist in the condenser. 

B. EXPERIMENTAL RECOMMENDATIONS 

1. -The correlation for the small diameter HEATEX element must be re-investigated. 

This may require obtaining additional data from CAL GAVIN or designing an 
experiment at the Naval Postgraduate School which can be used to obtain the re- 
quired data. 

2. The correlation used for the twisted tape insert also needs to be re-evaluated in 
light of the large inundation effects. As with the HEATEX, this may require that 
additional experiments be carried out at the Naval Postgraduate School. 

3. Once a reliable correlation is obtained, the experiments conducted on the 
copper/nickel finned tubes need to be repeated. 

4. Experiments need to be conducted using both titanium finned tubes as well as 
titanium KORODENSE tubes. This data can then be compared to that obtained 
with the copper/nickel finned tubes and the copper/nickel KORODENSE tubes to 
determine whether titanium tubes can deliver the required performance. 

5. Smooth copper tubes should be wrapped with different wire diameters at different 
pitches in order to determine an optimal pitch-to-wire diameter ratio. 

6. Experiments need to be conducted on wire wrapping of KORODENSE tubes which 
have a different pitch. These experiments should be based on results obtained from 
the smooth copper tube experiments. 

7. Experiments need to be repeated using R-114 and possible alternate refrigerants 
that the Navy is interested in such as R-24. 

8. A KORODENSE type tube which is roped on the inside and finned on the outside 
needs to be developed and tested since calculations have shown that this type of 
tube can lead to smaller condensers. 



144 



APPENDIX A. FLOWMETER CALIBRATION 



The flowmeters used to measure flow through the main condenser tubes were cali- 
brated in the following manner. The meters were disconnected from the condenser tubes. 
An empty 55 gallon drum w’as placed on a scale (resolution 0.5 lbs). The condenser 
pump was started with flow through the bypass. 

Flowmeters were calibrated at 0 °C and 24 °C. Each flowmeter was initially opened 
to maximum flow. The flow was then throttled back to 10%. The weight of fluid in the 
drum W’as recorded and the timer started. The w'eight was then recorded at the end of the 
calibration time. This procedure was repeated in 10% increments up to 100% flow. 
Flow'meter number 2 was only calibrated up to 90% flow since a mechanical stop pre- 
vented it from going to 100%. The calibration runs were then repeated from the maxi- 
mum flow back to 10% flow in 10% increments. The mass flow rate was then calculated 
in kg/sec based on an average of the two readings. The data for the flowmeter cali- 
brations as a function of nominal flow’ rate are shown in Table 4. 



Table 4. MASS FLOW RATE (KG/SEC) AT 0 AND 24 C. 





Flow Meter A 


Flow’ Meter B 


Flow Meter C 


Flow Meter D 


Nominal 
Flow Rate 
(%) 


0“C 


24*C 


OX 


24X 


OX 


24X , 


OX 


24X 


10 


0.0149 


0.0172 


0.0127 


0.0181 


0.0141 


0.0179 


0.0253 


0.0318 


20 


0.0363 


0,0403 


0.0308 


0.0408 


0.0374 


0.0391 


0.0513 


0.0613 


30 


0.0544 


0.0621 


0.0531 


0.0626 


0.0566 


0.0616 


0.0798 


0.0894 


40 


0.0794 


0.0851 


0.0739 


0.0844 


0.0793 


0.0899 


0.1129 


0.1216 


50 


0.1003 


0.112 


0.0953 


0.1 102 


0.1000 


0.1103 


0.1406 


0.1555 


60 


0.1179 


0.1339 


0.1 134 


0.1295 


0.1211 


0.1318 


0.1728 


0.1823 


70 


0.1388 


0.1556 


0.1379 


0.1506 


0.1427 


0.1552 


0.1975 


0.2114 


80 


0.1563 


0.1769 


0.1533 


0.1715 


0.1639 


0.1716 


0.2261 


0.2409 


. 90 


0.1846 


0.2028 


0.1799 


0.1912 


0.1813 


0.1975 


0.2461 


0.2667 


100 


0.2019 


0.2151 






0.2034 


0.2183 


0.2788 


0.2924 
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A least squares linear regression was performed for each flowmeter at each of the 
two temperatures for which the calibration runs were carried out. The linear regression 
equation takes the following form; 



m = a + bN 



where a is the y-intercept in kg/s, b is the slope in kg/s/% and N is the flowmeter setting 
in %. The results are shown in Table 5. 



Table 5. FLOWMETER CALIBRATION REGRESSION RESULTS. 





Flow Meter A 


Flow Meter B 


Flow Meter C 


Flow Meter D 




0 


24 *C 


0 ’C 


24 ”C 


0 "C 


24 


0 "C 


24 ”C 


Slope 


0.00206 


0.00225 


0.00204 


0.00216 


0.00208 


0.0022 


0.0028 


0.0029 


Y intercept 


-0.00414 


-0.00329 


-0.00615 


-0.00138 


-0.0036 


-0.00211 


-0.00122 


0.0048 



In order to obtain flowrates at any intermediate temperatures, a simplified interpo- 
lation procedure was used. The results of the calibrations for the four flowmeteres are 
shown in Figure 89, Figure 90, Figure 91, and Figure 92. 
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Figure 89. Calibration for flowmeter A. 
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Figure 90, Calibration for flowmeter B, 
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Figure 91. Calibration for flomiieter C, 
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Figure 92. Calibration for flowmeter D. 
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APPENDIX B. EMF VERSUS TEMPERATURE 



The EMF read from the thermocouples must be converted into a temperature for 
proposes of calculation. This is accomplished using a 7"* order polynomial which has the 
general form: 

r= C(0) + C(/)F^ 

w’here V is the EMF reading in volts, and I has values from 1 to 7. The leading coeffi- 
cients are given in Table 6. 

Figure 93 shows the relationship between EMF and temperature. 



Table 6. LEADING COEFFICIENTS FOR CONVERSION OF 



EMF TO TEMPERATURE. 



Leading Coefficient 


Value 


C(0) 


0.100S6 


C(l) 


25727.944 


C(2) 


-767345.8 


C(3) 


78025595.8 


C(4) 


-9247486589 


C(5) 


6.977.V10" 


C(6) 


-2.66x10'^ 


C(7) 


3.9a:10'" 
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Figure 93. EMF versus temperature. 
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APPENDIX C. PROPERTIES OF ETHYLENE GLYCOL/WATER 

MIXTURE 



The physical and thermodynamic properties of the ethylene glycol/water mixture 
were determined using equations and figures from Cragoc [Ref 56] and Gallant [Ref. 
57). 

The density of the ethylene glycol/water mixture is calculated in lines 2490 to 2530 
of the computer program. Density as a function of temperature for various concen- 
trations of ethylene glycol are showm in Figure 94. 

The kinematic viscosity of the ethylene glycol/water mixture is calculated in lines 
2370 to 2420. Figure 95 shows the kinematic viscosity plotted as a function of temper- 
ature for various ethylene glycol concentrations. 

The heat capacity is calculated in lines 2430 to 2480. Figure 96is a plot of heat ca- 
pacity as a function of temperature for difierent ethylene glycol concentrations. 

The thermal conductivity is calculated in lines 2560 to 2610 of the program. Values 
of the thermal conductivity as a function of temperature are plotted for various ethylene 
glycol concentrations in Figure 97. 
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Figure 94. Density of ethylene glycol/iiater mixture. 
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Figure 95. Kinematic viscosity of ethylene glycol/water mixture. 
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Figure 96. Thermal heat capacity of ethylene glycol/water mixture. 
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Figure 97. Thermal conductivity of ethylene glycol/water mixture. 
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APPENDIX D. PHYSICAL AND THERMODYNAMIC PROPERTIES OF 

R-113 



The physical and thermodynamic properties of R-113 were calculated based on 
equations supplied by Chapman [Ref 58). 

The heat capacity, C„, of R-113 was calculated from the equation; 

Cp = (929.0 -f 1.037) 

where T is in °C and C, is in kJ/kg K. Figure 98 shows C, as a function of temperature 
for R-113. 

The dynamic viscosity for R-113 is calculated from the equation: 

503 

_c (T-2AS) 

M= 1.34x10 10 

where T is in K and n is given in kg/m s. The relationship between n and temperature 
is shown in Figure 99. 

The thermal conductivity of R- 1 1 3 can be calculated from: 
k = 1.7308(4.846x10“^ -Tx 6.57x10“^) 

where T is in °F and K is in W/m K. Figure 100 shows the thermal conductivity as a 
function of temperature for R-113. 

The density of R-113 can be calculated from: 

p= 1.62x10^- Tx (2.219 -I- 7x 2.358x10“^) 

where T is in °C and p is in kgjm^. Density of R-113 is shown as a function of temper- 
ature in Figure 101. 

The latent heat of R-1 13 can be calculated from the equation: 
hfg = 2326.1(70.55 - t x (0.04838 -t- 1.26xl0“^7)) 

where T is in “7 and, is in J/kg. Figure 102 shows the relationship between latent heat 
and temperature for R-113. 
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Figure 98. Heat capacit) versus temperature for R- 113 . 
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Figure 99. Viscosity as a function of temperature for R-113. 
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Figure 100. Thermal Conductivity of R-113 as function of temperature. 
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Figure 101. Density as a function of temperature for R-113. 
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LATENT HEAT OF R-113 




Figure 102. Latent heat versus temperature for R-1 13. 
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APPENDIX E. SAMPLE CALCULATION 



A sample calculation was conducted on run 1 of data set SM06. This run was carried 
out on smooth copper tubes with the HEATEX insert installed. The sample calculation 
was performed for the top tube only. 

Temperatures were recorded as emfs using copper/constantan thermocouples. EiMF 
values were converted to temperatures using the equation showm in Appendix B. 
Table 7 summarizes the results. 



Table 7. EMF AND TEMPERATURES FOR RUN 1 DATA SET SM06. 



Thermocouple 


EMF (microvolts) 


Temperature (C) 


T{0) 


1908 


46.83 


T(l) 


1895 


46.54 


T{2) 


1909 


46.86 


T{3) 


1922 


47.16 


T{4) 


1918 


47.06 


T{5) 


1919 


47.099 


T{6) 


367 


9.44 


T(7) 


369 


9.496 


T(8) 


370 


9.52 


T(9) 


377 


9.7 


T{10) 


893 


22.52 


T{11) 


893 


22.53 


T{12) 


831 


20.99 


T(13) 


831 


20.99 


T(14) 


798 


20.18 


T(15) 


798 


20.18 


T{16) 


853 


21.54 


T{17) 


853 


21.54 
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For tube A, the temperature difTerential is calculated by: 



7'(10) + 7(11) 
2 



- 7(6)= 13.086 C 



The average temperature is then calculated from: 

Tavg = 7(6) + A7x 0.5 = 15.98 C 

Based on the ethylene glycol/water mixture physical properties (Appendix C), and using 
and average temperature of 15.98 C and an ethylene glycol concentration of 54%, the 
following fluid properties are obtained: 



Based on the calibration curves for flowmieter A, the mass flowrate of coolant cor- 
responding to a setting of 15.74 can be calculated. A simplified interpolation procedure 
is used to account for temperature differences. Based on this calculation, we obtain a 
mass flow rate of 0.2979 kg/s for this run. 

The coolant velocity can be calculated from: 



p,= 1071.3 Ag/m' 
fj-c = 0.0051 kgim s 
= 3196.99 Jl kg K 
K = 0.3742 WIm k 
Pr^ = 43.57 







where the inside cross-sectional area is given by: 




Ai = ^ (0.0 13259)^ = 0.000 1 38 



and 



0.02979 



= 0.201 mis 



" (1071.3)(0.000138) 
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The coolant Reynolds number, Re„ is calculated from: 



M/ 



(0.201X0.013259) 

Re. = = 559.89 

4.76jc10~^ 



The heat transferred to the coolant is calculated from: 



q = mCp_^T= (0.02979)(3196.99)(13.085) 
9= 1246.19 W 

and the heat flux is then calculated from : 







where 

= nd^L = k(0.0 15875X1.2 192) = 0.0608 rr? 
where L is the condensing length of the tube. Based on this equation, we obtain; 

" _ . 1^46 .19_ 20494 9 Wjni^ 

In order to calculate the inside and outside heat transfer coefficients, the physical 
and thermodynamic properties of R- 1 1 3 must be calculated (Appendix D). These fluid 
properties are based on a reference or condensate film temperature. In order to obtain 
this film temperature, the wall temperature is needed. To obtain this, an iteration is 
carried out with an initial value of /i„ assumed (in this case 1000 K). The wall 

temperature is then calculated from: 



= 46.925- 



20494.9 

1000 



26.43°C 
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The film temperature is then calculated from: 



Tf=' 



r, 



SQ t I 2 

3“ T ‘ 



wall 



or 



Tf= ^ (26.43) = 33.26“C 

Based on the equations of fluid properties for R-113 (Appendix D), the following values 
were obtained for the condensate: 



/,^ = 147036.5 Jlkg 
k = 0.0734 Wjm K 
p = 1544.34 kglni^ 
p = 6.029j:10'^ kgjm s 

Using these fluid properties together with equation (5-15), we obtain: 

= 1907.06 WjtnK 

Using an iterative procedure, we adjust the initial guess of /i„ by taking the average of 
the calculated and guessed value of h^. The wall and film temperatures are then recalcu- 
lated and new fluid properties for R-113 are obtained. A new value of is then calcu- 
lated. This procedure continues until the newly calculated and previosuly calculated 
values are within 0.1%. For this run, we obtained: 

1251 IVIrn^ K 

The next step is to find a value for the overall heat transfer coefficient, t/„, where: 

r — ^ ° 

° LMTD 



where, for this run, 



LMTD = 30.48° C 



and 
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,, 20494.9 2 

= ' 3o ' 48 ~ = ^ 

The inside heat transfer coefficient, h„ can be obtained from a correlation for the 
Nusselt number for single phase coolant flow in a tube. This Nusselt number is obtained 
from a correlation (in this case) for the HEATEX element where: 

T^u, = 0 . 226 {Ref^\Prf^^ 

Nu, = 0.226(559.89)-^^(43.75)-^® = 78.41 



and 



_ NuA _ (78.41)(0.3742) 

”'■ di 0.013259 

/i;= 2213.07 IVIm^ K 

The outside heat transfer coefficient is calculated from equation (5-12) where: 



ho = - 



C/. 



do 

— — — R 4 

dihi 



The wall resistance, R„, is defined as: 



R 



m 




2Lk 



, f 0.015875 \ 
V 0.013259 ) 
2(1.2192)(386) 



1.91x10'^ KjW 



Therefore, 



ho = 



0.015875 



672.49 (0.0 13259)(22 13.07) 



-(1.91jc10 )(0.0608) 



= 1070.2 lE/m^ K 
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Figure 103 is a copy of the computer printout for this run. The small differences 
seen in the computer calculations versus these sample calculations are the result of 
roundoff error. 
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Data set number = 1 



RESULTS FOR TUBE NO. 1 



Mass flow rate = 
Inlet temperature = 
Saturation temp (Deg C> = 
DELT temp Dif. = 
Log. Mean temp Di f . = 
Heat fluA = 
Resistance of Tube Metai= 
Prandtl number = 
Reynolds number = 
Reynolds number H&B S = 
Ins ide h . t . c . = 
Inside NUSSULT NO. 

OUERALL H.t.c. < Uo ) 



2.980E-02 
9. 443E+00 

4. B93E+01 
1 .307E+01 
3.048E+01 
2.048E+04 
1 .913E-04 
4.359E+01 

5. B08E+02 
B.745E+02 
2.21BE+03 
7.852E+01 
B.720E+02 



TUBE 

t 

1 



FM UE6 
(%) (m/s) 
15.74 0.20 



DELT Uo HO 

(K) (W/m''2.K) 

13.074 6.720E+02 1.0B8E+03 



Figure 103. Computer printout of results for run 1 data set SM06. 



HNUS 

1 .251E+0 
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APPENDIX F. UNCERTAINTY ANALYSIS FOR Uo 



Run number 1 of data set SM06 was also chosen for the uncertainty analysis. The 
analysis was only conducted on l/„. The measured and calculated parameters found in 
the sample calculation (Appendix E) were therefore used in this section. In order to as- 
sess the effects of coolant velocity on the uncertainty in the measurement of U, , an 
uncertainty analysis was also performed on run number 1 1 of data set SM06. This run 
had a coolant velocity of 1.2 m/s as opposed to 0.2 m/s in run number 1. 

The uncertainty analysis utilizes the procedure suggested by Kline and McClintock 
[Ref 54]. This procedure states that if 



where x„ is the measured variable and 8x„ is the uncertainty in the measured variable. 

The uncertainty in the temperature measurements is due to uncertainty in the volt- 
age measured by the thermocouples. For the thermocouples used corresponds to 0.1 
°C. The uncertainty in in this thesis, this value was assumed to be 4 microvolts or 0.1 
C. The uncertainty in AT, and coolant inlet and outlet temperatures, T„ and , 
can be calculated as follows: 



R = .^n) 



then the uncertainty in R, namely, 5R is given by the equation: 




r(10) + 7(11) 
2 



Taking the partial derivatives we obtain: 



= J. 

a 7(10) 2 



likewise. 
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The uncertainty in the measurement of 7„ is thus; 



<5 7 ;, 



57(10) 




+ 



57 



CO 



57(11) 



57 



o.s 



Based on the thermocouple uncertainty, 57 = 0.1 C, we obtain: 

57,^ = 0.0778 C 

The coolant inlet temperature is calculated from a single thermocouple, number 6 for 
tube A. The uncertainty in tliis measurement is calculated as: 

oTci _ 

57(6) 



and. 



^7;,= 



57 



ci 



57(6) 



57 



0.5 



or. 



57„ = 0.1 C 



Finally, the uncertainty in A 7 is calculated from: 

5A7=[(57j' + (57,,)'f 



and. 



5A7= 0.1347 C 

Coolant velocity is calculated from: 



V =-72_ 

^ pA 

In this run, m was calculated to be 0.02979 kg/s. The uncertainty in the reading of the 
flowmeter is defined as the scale interpolation. This value corresponds to one-half the 
value of the smallest marked increment. In this case, the scale interpolation factor, 5,„ 
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is 0.5% which corresponds to 0.00101 kg/s resolution. In addition, a time-wise jitter, 
of 0.5% was noted during the runs. Therefore, the uncertainty in m is given by; 



6m 

rh 







or 



= 0.0479 
m 



The uncertainty in the cross-sectional area of the tubes was estimated from: 



6Ai 

Ai 



sd : 



where 6df is given as 0.1 mm based on tolerances supplied from the manufacturer. 
Therefore, 



6Ai 



= 0.0107 = 1.07% 



The uncertainty in the coolant velocity is now given by: 



Vc 



^ m > ^ A, ^ 



- 10.5 



or 



-^ = 0.049 = 4.9 % 

' c 
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The uncertainty in the calculation of coolant Reynolds number is given by: 



Re, 



[>f 



2 2 



no. 



5 



or substituting appropriate numbers we find 



5 Re, 

= 0.0497 = 4.97 % 
Re, 



The uncertainty in the calculation of the heat transferred to the coolant is related 
to the uncertainty in the measurements of coolant mass flow rate and coolant temper- 
ature rise where: 






[' 



(i^)2_^( 
m 



SAT 

AT 




0.5 



or 



= 0.0489 = 4.89 % 

q 



The uncertainty in the heat flux is then: 

^q’o ^ ^ >,2 , iio 

q'o ^ q ^ A, 



f 



where the uncertainty in is given by: 



Ao 






- 10.5 



Here, the uncertainty in the length of the tube, SL , is assumed to be 5 mm while Sd, is 
assumed to be 0.1 mm (again from given tolerances). Substitution of the appropriate 
numbers yields: 



Sq’’ 

= 0.0496 
q o 
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Finally, since (/„ is given by; 



^ LMTD 

the uncertainty in the LMTD must be calculated. The LMTD is given by: 

LMTD = — 

In 

In order to calculate the uncertainty in the measurement of the LMTD, the following 
equations were used: 



AT 

T — T 

‘ sat ‘ c,o 



dLMTD 

BAT 




1 

~ Tc,i 
T — T 

^ sat '■ c,o 



BLMTD 

dTsat 






In 



T —T- 

^ sat ^ c,i 

~f — T 

^ sat * c,o 



i^sat T'c,l){T'sat T'c,d 



dLMTD AT 



dLMTD 


In 


T — T ~] 

^ sat ^ c,i 

T — T 

' sat c,o 

AT 


"2 

i^sat ~ '^c,o) 




In 


~ T _ 7 . 1 
^ sat ^ c,i 




2 

iTsat-TJ 


T — T 

^ sat ^ c,o 



The uncertainty in the LMTD is therefore given by: 



dLMTD = 






dLMTD 






2 , f dLMTD \2 



C,0 



STcj 



For the sample run with a coolant velocity of 0.2 m/s, the uncertainty in the measure- 
ment of the LMTD is: 



21 2 



SLMTD = [[(2.329)(0.1347)]^ 4- [( -1.0I57)(0.0778)]^ + [( -2.9097)(0.0778)]^ -t- [(1.894)(0.1)]^] 
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or 



SLMTD = 0.446 K 



and 



SLMTD 

LMTD 



0.446 

30.48 



= 0.0146= 1.46% 



The uncertainty in the overall heat transfer coefficient, is then given by; 

-\ 1 



<5o; 



Sq 



+ 



( SLMTD \i 
\ LMTD ) 



or 



^ = [(0.0496)^ + (0.0146)^]i 



and 



SU. 

= 0.052 = 5.2 % 

K 

Table 7 summarizes the results of the uncertainty analysis for the smooth copper tubes 
with the HEATEX insert for coolant velocities of 0.2 and 1.2 m/s respectively. 

The uncertainty analysis neglects any uncertainty in the physical properties of the 
coolant or the R-113. For the most part, this is probably a reasonable assumption. The 
exception to this may be the viscosity due to its sensitivity to temperature. Hence, the 
uncertainty reported here is an underestimate. 

It is curious that the relative magnitudes of the uncertainty for the low and high 
coolant velocities are essentially the same. Two different mechanisms dominate in the 
calculated uncertainty at low and high coolant velocities. In terms of the mass flow rate, 
the time-wise jitter and scale interpolation terms are the same for both coolant velocities. 
However, the actual mass flow rate is less at the lower flow rates. Hence, the uncertainty 
in the mass flow rate calculation is more significant at the lower coolant velocities. At 
higher flow rates the uncertainty in the calculation of the LMTD becomes the dominant 
term. This is because the coolant temperature rise is smaller at high coolant velocities 
but the uncertainty in the temperature measurements remains the same. The mass flow 
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Table 8. SUMMARY OR RESULTS FOR SAMPLE UNCERTAINTY ANALY- 
SIS. 



Variable 


V, =0.2 m/s 


V, = 1.2 m's 




o.rc 


O.rc 


ST,, 


0.078“ C 


0.078“ C 


ST„, 


0.078“C 


0.078“C 


SAT 


0.135“C 


0.135“C 


dm 

m 


0.048 


0.008 


V 1 


0.011 


0.011 


SV, 

V. 


0.049 


0.013 


SRe, 

Re. 


0.05 


0.015 


Sq 

<7 


0.05 


0.015 


SA, 

A. 


0.008 


0.008 


Sq\ 

q". 


0.05 


0.017 


SLMTD 

LMTD 


0.015 


0.054 


SD\ 


0.052 


0.057 



rate goes up so that the uncertainty in the coollant mass flow rate calculation decreases. 
Thus, at the higher coolant velocities, the uncertainty in the calculation of the LMTD 
is the dominant term. 
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APPENDIX G. PROGRAM LISTING 
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1000 I FILE: 
1005 ! PURPOSE: 
ube bundle. 



DRPC0N7 

This program collects and processes condensation data for the t 



1015 

1020 

1021 

1022 

1023 

1024 
1026 
1030 
1035 
1040 
1045 
1050 
1055 
1060 
1065 
1070 
1075 
1080 
1085 
1090 
1095 
1100 
1105 
1110 
1115 
1120 
1125 
1130 
1135 
1140 
1145 
1150 
1155 
1160 

1215 

1216 

1217 

1218 
1220 
1225 
1230 
1235 
1240 
1245 



! CREATED: NOUEMBER 2, 1988 
! UPDATEl: APRIL 25 1991 <5. MEMORY) 

! UPDATE2: May 17 1991 <S. MEMORY) 

! UPDATE3: JULY 26. 1991 <R.W, MAZZONE) 

! UPDATE4: 13 AUGUST 199KR.U. MAZZONE) 

! UPDATES: 10 OCTOBER 1991 (R.W. MAZZONE) 

! CHANGE H0D=3 AND USING CALCULATED Cl FROM WILSON 
BEEP 



'5 

•6 



PRINTER IS 1 
PRINT USING ■'4.x, 
PRINT USING "6X,' 
PRINT USING ’'GX,"”1 
PRINT USING "ex, “"2 
PRINT USING "GX, ’”'3 
PRINT USING "GX/'M 
PRINT USING "BX, 
PRINT USING ”GX,' 
PRINTER IS 701 
INPUT Icall 
IF lcall=0 THEN CALL 
Icall=l THEN CALL 
Icall=2 THEN CALL 
Icall=3 THEN 
Icall=4 THEN 
Icall=5 THEN 
IcalI=B THEN 



SELECT OPTION 

0 TAKING DATA OR REPROCESING PREUIOUS DATA 

1 PLOTTING H US DELTA-T 

2 PLOTTING HRAT US N 

PLOTTING WILSON 

PURGE FILES 

XYREAD 

NUSSELT ESTIMATE' 



IF 
IF 
IF 
IF 
IF 
IF 
END 
SUB 
COM 
COM 
COM 
COM 
COM 
DIM 
DATA 



Mai n 
Plot2 
Plotl 
CALL Plot3 
CALL Purge 
CALL Xyread 
CALL Nusselt 



Main 

/Cc/ C(7) 

/Fid/ Ift 

/Nus/ Iin,T5at,Qdpl,Hnus,Kf,Rhof,Hfg,Muf,Do,Itube 
/Will/ Doa<4),Dia<4),Kna<4),Iact,Droot(4 ) 

/Wil2/ Delta, I sat , Nset s , Hod ,Cia( 3 , 3 ) , Alpaa( 3 ) 
Emf<17),Tp<4),T( 17),Ho(4),Qdp(4),Uo(4),Pc<4) 

1 . 0 , 1 , 0 , 1 . 0 , 1.0 



DATA 5.172,5,172,5.172,5.172 
DATA 0.226,0.226,0.226,0.226 
DATA 0,063,0.063,0.063,0,063 
READ Cia(*) 

DATA 0,10086091,25727.94369,-767345.8295,78025595.81 
DATA -9247486589, 6, 97688E1 1 ,-2,661 92E13, 3. 94078E14 
READ C( * ) 

DATA 0,015875,0.014000,0.015798,0,0159,0,0 
DATA 0,013259,0,010160,0.012776,0,01346,0.0 
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1250 DATA 386.0,42.975,42.975,21.9,0.0 

1251 DATA 0.0,0.014,0.0, .0141 ,0.0 

1255 DATA 0.0005588,3 

12B0 READ Doa( •), Dia( *) ,Kma< *), Drool (*), Del ta , Hod ! Hod=H/Di 

1285 L=1.2192 ! Condensing length 

1270 Jset=0 

1271 Okaccpt=l 

1275 BEEP 

1280 INPUT "ENTER MONTH, DATE AND TIME ( MM= DD : HH : MM : S5 ) " ,01 g$ 

1285 OUTPUT 709: “TD" iOtg$ 

1290 BEEP 

1295 INPUT "SELECT OPTION (0=DAQ, l=FILE)",In 

1300 Ihard=l 

1305 BEEP 

1510 INPUT "WANT A HARDCOPY PRINTOUT < 1=DEF=YES,0=N )" , Ihard 
1315 BEEP 

1320 INPUT "SELECT < 0=R-1 14 , 1=STEAM , 2=R-1 13 , 3=EG ) " , I f t 
1325 Iin=l 

1530 Isat=2 

1335 BEEP 

1340 INPUT "SELECT SAT TEMP MODE ( 0=LIQ , 1=UAP , 2=( LIQ+UAP )/2=DEF ) " , I sat 
1345 IF Ihard=l THEN PRINTER IS 701 

1350 — IF ln=0 THEN 



1355 BEEP 

1360 INPUT "GIUE A NAME FOR THE NEW DATA FILE",File$ 

1365 CREATE BOAT File$,20 

1370 ASSIGN ©File TO File$ 

1375 BEEP 

1380 INPUT "ENTER TUBE CODE (0=SMOOTH CU,1=FINNED CU ,2=K0R0DENSE , 3=TIT FIN 



) " , 1 1 ube 



1381 


BEEP 






1382 


INPUT 


"0=NO INSERT, 1=TWISTED 


TAPE, 2=HITRAN1, 3=HITRAN 


1385 


BEEP 






1390 


INPUT 


"ENTER EG CONCENTRATION 


(WT PERCENT )" ,Egrat 


1395 


ENTER 


709;Dtg$ 




1400 


OUTPUT 


©Fi le:Dtg$ 




1405 


OUTPUT 


©File: Itube,Egrat,Ddl, 


Dd2,Dd3,Dd4,Dd5 


1410 


Iact=10 




1415 


BEEP 






1420 


INPUT 


"SELECT (0=TOP, 1=SEC0ND 


10=BUNDLE=DEFAULT )" , 


1425 . 


PRINT 






1450 


PRINT 


USING "10X,"”FILE NAME: 


" " ,12A" :File$ 


1435 


PRINT 






1440 


ELSE 






1445 


BEEP 






1450 


INPUT 


"ENTER NAME OF EXISTING 


FILE" , Files 


1455 


BEEP 






1460 


INPUT 


"ENTER NUMBER OF DATA SETS STORED" ,Nsets 


1465 


I w i 1 = 1 






1470 


BEEP 







2 " , Insert 



I act 
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1475 INPUT "WANT TO CALL WILSON? ( 1=DEFAULT=YES , 0=NO ) " , Iwi 1 

1480 IF Iwil=l THEN 

1485 BEEP 

1490 INPUT "WHICH TUBE ( 0-TOP , l=SECOND 10=BUNDLE lact 

1491 INPUT "0=NO INSERT, 1=TWISTED TAPE, 2-HITRANl, 3-HITRAN2" , Insert 

1493 CALL Wilson 

1494 IF Iact=10 THEN 

1495 FOR lcount=0 TO 3 

1498 Cia(Insert,Icount)=Cia<Insert,0) 

1497 Alpaa( Icount )=Alpaa( 0 ) 

1498 NEXT Icount 

1499 END IF 

1500 IF Ihard=l THEN PRINTER IS 701 

1505 PRINT 

1510 PRINT USING "10X,""FILE NAME-- " " , 12A" ;File$ 

1515 PRINT 

1520 IF Iact>9 THEN 

1521 FOR 11=0 TO 3 

1522 PRINT USING " 10X ," "RESULTS FOR TU8E NO . " " , 3D" : 1 1 + 1 

1523 PRINT 

1525 PRINT USING " 10X ,"" INSIDE COEFFICIENTS: " " , DD . 3D" ; Cia< Insert , 1 1 

) 

1530 PRINT USING "10X, ""ALPHAS: " " , 4(1X ,MZ . 5DE , 2X )” ; Alp 

aa( II ) 



1531 NEXT II 

1535 ELSE 

1540 PRINT USING " 10X INSIDE COEFFICIENT FOR TUBE "",D,"": "",DD.3 

D" ;Iact + l,Cia< Insert ,Iact ) 

1545 PRINT USING "10X, ""ALPHA FOR TUBE "",D,"": "",DD.3 

D" ; Iact + 1 , Alpaa( lact ) 

1550 END IF 

1555 ELSE 

1560 BEEP 

1581 INPUT "0=TOP, 1-SECOND, 10-BUNDLE" , lact 

1582 INPUT "0-SMOOTH CU , 1=FIN CU , 2-KOROD, 3=FIN TITANIUM" , I tube 

1583 INPUT "0=NO INSERT, 1-TWISTED TAPE, 2-HITRANl , 3-HITRAN2 ", Insert 

1585 INPUT "ENTER Cl UALUE " ,Cia< Insert , II ) 

1570 END IF 

1575 ASSIGN ©File TO Fiie$ 

1580 ENTER @File:Dtg$ 

1585 ENTER @Fi le : I tube , Egrat ,Ddl , Dd2 , Dd3 , Dd4 , Dd5 

1588 INPUT "0=SMOOTH CU , 1=CUNI FIN, 2=KORODENSE, 3=FIN TITAN ",Itube 
1590 END IF 

1595 Iout=l 

1600 BEEP 

1605 INPUT "WANT TO CREATE AN OUTPUT FILE? ( 1=DEF=YES , 0=N0 ) " , lout 

1610 IF Iout=l THEN 

1615 BEEP 

1620 INPUT "ENTER A NAME FOR OUTPUT FILE",Fout$ 

1625 CREATE BOAT Fout$,5 



181 



1630 

1635 

1640 

1645 

1646 
1650 
1655 
1660 
1665 

1670 

1671 

1672 

1673 

1674 

1675 

1676 

1677 

1679 

1680 
1682 
1683 
1685 

1690 

1691 
/s )? 

1694 

1695 

1696 

1697 

1698 

1699 

1700 

1701 

1702 

1703 

1704 

1705 

1706 

1707 

1708 

1709 

1711 

1712 

1713 

1714 

1715 

1716 

1717 

1718 

1719 



ASSIGN @Fout TO FoutS 
END IF 

! 

Do=Doa ( I tube ) 

Dol=Doa< 1 1 ube ) 

Di=Dia( Itube ) 

Km=Kna( Itube ) 

Ax=PI*Di ''2/4 ! Cross-sectional area 

Ao=Pl*Do*L 

! Rm=Do«L 0G( Do/Di )/( 2*Kn ) 

IF Itube=l THEN 

Dol=Droot< Itube ) 

END IF 

IF Itube=3 THEN 

Dol=Droot < Itube ) 

END IF 
Ao=PI *Do*L 

Rn=LOG(Dol/Di )/<2*L*Km) 

Al=PI*Dol*L 

I 

IF lm=0 THEN 
PRINTER IS 1 
BEEP 

INPUT "WHAT IS REQUIRED NOMINAL COOLANT UELOCITY IN n/s (MAX. = 1.4 m 

,Uegst 

OUTPUT 709: "AR AF6 AL9 UR5" 

FOR 1=6 TO 9 

OUTPUT 709; "AS SA“ 

Usum=0 

FOR J=1 TO 5 
ENTER 709 :E 
Usupi=Usun+E 
NEXT J 

Enf< I )=Usun/5 

T( I ) = FNTvsv(Enf< I > ) 

NEXT I 

Tavgl = (T(6 ) + T<7 ) + T(8) + T< 9 ) )/4 
IF Tavgl<0 OR Tavgl>24 THEN 

PRINT "COOLANT TEMP. NOT IN RANGE 0 TO 24 C" 

GOTO 1691 
END IF 

Rhoegl=FNRhoeg< Tavgl , Egrat ) 

Mdotl=Rhoegl*Ax«Uegst 
Pca0=(Mdotl+. 00414 )/. 00206 
Pca24=(Mdotl+. 00329)/. 00225 
Pcb0=( Mdotl+. 006 15 )/. 00204 
Pcb24=(Mdotl + . 00138 )/.00216 
Pcc0=(Mdotl + .0036 >/. 00208 
Pcc2A=(Mdotl+. 00301 )/. 00223 
Pcd0=(Mdotl+. 00122 )/.0028 
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1720 Pcd24=(Mdotl-.0048 )/.0029 

1721 Pc<0)=Pca0+<Tavgl*(Pca24-Pca0))/24 

1722 Pc( 1 )=*Pcb0+< Tavgl»< Pcb24-Pcb0 ) )/24 

1723 Pc<2)=Pcc0+(Tavgl*(Pcc24-Pcc0))/24 

1724 Pc(3)*Pcd0+(Tavgl»<Pcd24-Pcd0) )/24 

1725 PRINT "SET FLOWMETER READINGS CORRESPONDING TO:" 

172B PRINT USING "BX,""% OF METER A = " " , DDD . D" ; Pc < 0 ) 

1727 PRINT USING "BX,""% OF METER B = " " ,DDD.D" :Pc< 1 ) 

1728 PRINT USING "BX,""% OF METER C = " “ , DDD . D“ : Pc< 2 ) 

1729 PRINT USING "BX,”"% OF METER D = " " , DDD. D " : Pc < 3 ) 

1730 PRINT 

1731 PRINT "HIT CONTINUE WHEN READY" 

1733 PAUSE 



1734 Repeat = ! 

1735 ! PRINTER IS 701 

173B Dtld=47.5 ! Desired temperature of liquid 

1737 Ido=l 



1738 

1739 

1740 

1741 

1742 

1743 ! 

1744 

1745 
174G ! 

1747 

1748 ! 

1749 

1750 ! 

1751 

1752 

1753 ! 

1754 

1755 
175B 

1757 

1758 
el2 

1759 
17B0 ! 
17B1 
17B2 
17B5 
17B4 
17B5 

1766 

1767 

1768 

1769 



ON KEY 0,15 RECOUER 1734 

PRINT USING "4X, ""SELECT OPTION 

PRINT USING "6X,""0=TAKE DATA 

PRINT USING "6X,""1=SET Tsat (DEFAULT) 

PRINT USING “4X,""N0TE: KEY 0 = ESCAPE 

Ido = desired option 
BEEP 

INPUT Ido 

Set default value for input 
IF Ido>l THEN Ido=l 
Take data option 
IF ldo=0 THEN 1797 

Loop to check saturation temperature 
IF Ido=l THEN 

INPUT "ENTER DESIRED Tsat ( DEFAULT=47 . 5 C 
INPUT "ENTER DESIRED Tsat <DEFAULT=2.2 C 
Nn=l 

Nrs=Nn MOD 15 
Nn=Nn+ 1 
IF Nrs=l THEN 

PRINT USING "4X,“" DTsat Tldl Tvab 
Del3 Del4 
END IF 

Read thermocouple voltages 
OUTPUT 709; "AR AF0 AL17 UR5" 

Nend*17 
FOR 1=0 TO 5 

OUTPUT 709; "AS SA" 

Usum=0 

FOR J=1 TO 5 
ENTER 709;E 
Vsum=Usum+E 
NEXT J 



R-113 ) 
R-114 )" 



Tvat 



’ ,Dtld 
,Dtld 



Dell D 
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1770 

1771 

1772 

1773 

1774 

1775 

1776 

1777 

1778 

1779 

1780 

1781 

1782 

1783 

1784 

1785 

1786 

1787 

1788 

1789 

1790 

1791 

1792 

14 

1793 

1794 

1795 

1796 

1797 

1798 

1799 

1800 

1801 

1802 

1803 

1804 

1805 

1806 

1807 

1808 

1809 

1810 

1811 

1812 

1813 

1814 

1815 

1816 

1820 

1825 



Enf ( I ) = lJsum/5 

T( I )=FNTv5v(Emf< I ) ) 

NEXT I 

FOR 1=6 TO Nend 

OUTPUT 709: "AS SA" 

Usuro=0 

FOR J=1 TO 20 
ENTER 709 ;E 
Usum=Usum+E 
NEXT J 

Emf( I )=Usum/20 
T( I )=FNTvsv<Enf( I ) ) 

NEXT I 

Compute average temperature of liquid 
Tliq = (T<3) + T(4) )/2 
Tvap=<T<0)+T( 1 )+T<2 ) )/3 
Tvat=(T<0)+T< 1 ))/2 
Tliql=T(5 ) 

Dell = ((T<10)+T(ll ))/2)-T<6) 

Del2=((T(12)+T(13))/2)-T(7) 

Del3=( (T< 14 )+T( 15 ) )/2 )-T<8) 

Del4=< (T( 16 ) + T( 17 ) )/2 )-T(9) 

PRINT USING "4X,8(MDD.DD,3X )" ;Dtld,T(3),T(2 ),Tvat ,Dell ,Del2,Del3,De 



WAIT 2 
GOTO 1755 
END IF 

TAKE DATA IF lm=0 LOOP 
PRINTER IS 701 

OUTPUT 709; "AR AF0 AL17 UR5" 

Nend=17 ! INCREASE TO 20 IF FIUE TUBES IN BUNDLE 
FOR 1=0 TO 5 

OUTPUT 709; “AS SA" 

Usum=0 

FOR J=1 TO 5 

ENTER 709 :E 
Usum=Usum+E 
NEXT J 

Emf < I )=Usum/5 
NEXT I 

FOR 1=6 TO Nend 

OUTPUT 709; "AS SA" 

Usum=0 

FOR J=1 TO 20 
ENTER 709;E 
Usum=Usum+E 
NEXT J 

Emf( I )=Usum/20 
NEXT I 
ELSE 
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1830 
1835 
1840 ! 
1845 ! 
1850 I 
1855 
1860 
1865 
1870 
1875 
1880 
1881 
1885 
1890 

1895 

1896 
1900 

1905 

1906 

1907 

1908 

1925 

1926 
1930 
1935 
1940 
1945 
1950 
1955 
1960 
1965 
1970 
1975 

1980 ! 

1981 
1985 

1990 

1991 ! 

1992 ! 

1994 ! 

1995 I 

1996 
2000 
2005 
2010 
2015 
2020 
2025 
2030 
2035 



ENTER ©File:Pc( * ),Enf ( * ),Tp< * ) 

END IF 

DATA ANALYSIS 
Nend=17 

FOR 1=0 TO Nend 

T( I )=FNTvsv(Emf < I ) ) 

NEXT I 

Tvap=(T(0)+T(l)+T(2))/3 
Tliq = (T<3) + T<4) )/2 
Tliql=T(5) 

IF lset=0 THEN 
Tsat =T 1 iq 1 
END IF 

IF Isat=l THEN 
T sal =Tvap 
END IF 

IF Isat-2 THEN 

Tsat = ( T vap+Tl iq )/2 
END IF 
Jset=Jset+l 

IF Okaccpt=0 THEN Jsel=Jset-l 
PRINT 

PRINT USING "10X,"”Dala set number = "".DD'' 

PRINT 

lbeg=0 

Iend=3 

IF Iact<10 THEN 
Ibeg=Iact 
Iend=Iact 
END IF 

FOR I=Ibeg TO lend 

Gred=FNGrad( Enf ( 1+6 ) ) 

Tout=<T( 2*1+10 )+T(2*I+ll ) )/2 
Delt=Tout-T( 1+6 ) 

Tavg=T( I+B )+Delt* .5 
IF I=Iend THEN 
Delt=T( 10)-T(9) 

Tavg=T( 10 )-T( 9)/2 
END IF 

Rhoeg=FNRhoeg< Tavg , Egrat ) 

Nueg=FNNueg( Tavg , Egrat ) 

Mueg=Nueg*Rhoeg 
Cpeg=FNCpeg( Tavg , Egrat ) 

Keg=FNKeg(Tavg, Egrat ) 

Preg=Cpeg*Mueg/Keg 
Mdot=FNFmcal(I ,T( 1+6 ),Pc( I ) ) 
Veg=Mdot/(Rhoeg*A>; ) 

Reeg=Ueg*Di /Nueg 



: J set 
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2040 


Re5=4*t1dot/(PI«Mueg*(Di-4* Delta > ) 




2045 


Qdot=Mdot»Cpeg*Delt 




2050 


Qdp( I )=Qdot/Ao 




2055 


IF 1=0 OR I=Iact THEN 




2060 


Qdp l=Odp ( I ) 




2065 


CALL Nusselt 




2070 


END IF 




2075 


Li*itd=Delt/LOG< (Tsat-T< 1+6 > >/(Tsat-T< I+B )-Delt ) ) 




2080 


Uo< I )=Qdp( I )/Lmtd 




2081 


Cfeg=l 




2085 


IF Reeg<4000 THEN 




2086 


IF lnsert=0 THEN 




2088 


Bbl= . 0668*( Di /L )*Reeg*Preg 




2089 


Bb2=l+.04*< (Di/L )*Reeg*Preg )'• .6666 




2090 


Nueg=3.66+(Bbl/Bb2 ) 




2091 


END IF 




2092 


IF Insert=l THEN 




2094 


Nueg=l ,86*< Preg*Reeg )'' , 333333* < Di/L .33333*1 


[^f eg 


2095 


Nueg=Cia( Insert, I )♦( 1+5. 484E-3^Preg'' . 7* ( Re 5 / Hod 1 . 25 . 5 


2096 


END IF 




2097 


IF Insert=2 THEN 




2098 


Nueg=Cia< Insert , I )*<Reeg'’ .65 >*Preg'' .46 




2099 


END IF 




2100 


IF Insert=3 THEN 




2101 


Nueg=Cia< Insert,! >*<Reeg''.76 )*Preg'‘ . 46 




2102 


END IF 




2104 


ELSE 




2105 


BEEP 




2106 


PRINT USING "lOX," "TURBULENT CORRELATION 




2110 


Nueg= .027*Reeg''.8*Preg''.3333*Cfeg 




2115 


END IF 




2120 


Hi=Nueg*Keg/Di 




2125 


Hod )=!/( l/Uo( I )-Do/( Di*Hi )-Rn*Al ) 




2130 


IF 1=0 OR I=Iact THEN 




2131 


PRINT 




2132 


PRINT 




2134 


PRINT USING "10X," "RESULTS FOR TUBE NO."", 3D 


" : I + l 


2135 


PRINT 




2136 


PRINT USING "10X,""Mass flow rate = 


" " ,MZ.3DE" :Maot 


2137! 

\ 


PRINT USING "10X, ""Inside Tube Dia. (m.) 


" " , MZ . 3DE “ ; Dia( I tube 


f 

2138! 


PRINT USING "10X,"“180 DEG OUER Dia. < HOD ) = 


"" ,MZ.3DE" :Hod 


2140 


PRINT USING "10X, ""Inlet temperature - 


"",MZ.3DE“;T(I+6) 


2145 


PRINT USING " 10X ," "Saturation temp (Deg C) = 


" " ,MZ.3DE" :Tsat 


2146! 


PRINT USING "10X, ""DELTA Tape Thickness = 


"",MZ.3DE"; Delta 


2147 


PRINT USING "10X,""DELT temp Dif. 


" " ,MZ.3DE" ;Delt 


2148 


PRINT USING "10X,""Log. Mean temp Dif. 


" " , MZ . 3DE" ; Lmt d 


2150 


PRINT USING "10X,""Heat flux 


"",MZ.3DE";Qdpd> 


2151 ! 


PRINT USING " 10X ," "Conduct ivity E.G. = 


" " ,MZ.3DE" :Keg 


2152! 


PRINT USING " 10X ," "Conduct i VI ty Tube Metal = 


" ",MZ.3DE":Km 
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2153 


PRINT 


USING 


"10X, " 


"Resistance of Tube Metal 


= 


" " .MZ.3DE" ;Rm 


2155 


PRINT 


USING 


“10X, " 


"Prandtl number 


= 


" ",MZ.3DE";Preg 


2160 


PRINT 


US INC- 


"10X, " 


"Reynolds number 


= 


” " .MZ.3DE” ;Reeg 


2161 


PRINT 


USING 


"10X, " 


"Reynolds number H&B 5 


= 


"“,MZ.3DE";Re5 


2165 


PRINT 


USING 


“10X, " 


"Inside h. t .c . 


= 


” " .MZ.3DE” ;Hi 


2166 


PRINT 


USING 


"10X, “ 


"Inside NU5SULT NO. 


= 


" ",MZ.3DE'';Nueg 


2167 


PRINT 


USING 


“10X, " 


"O^JERALL H.t.c. ( Uo ) 


= 


,liZ.3DE" ;Uo( I ) 


2170 


PRINT 












2171 ! 


PRINT 


USING 


"10X, " 


"KF 


= 


.MZ.3DE" ;Kf 


2172! 


PRINT 


USING 


"10X, “ 


"RHOF 




,MZ.3DE" :Rhof 


2173! 


PRINT 


USING 


"10X, " 


"HF6 


= 


"" .MZ.3DE" :Hfg 


2174! 


PRINT 


USING 


"10X, " 


"MUF 


= 


"“,MZ.3DE" :Muf 


2175 


PRINT 












2177 


PRINT 


USING 


"10X, " 


"TUBE FM UEG 




DELT Uo 


HO 


HNUS 












2180 


PRINT 


USING 


“10X, " 


" t (%) (m/s) 




(K) <U/n''2 



K ) 

2185 PRINT USING " 10X , 3D , 4X , 3D . DD , 3X , Z . DD , 4X , DO . 30 , 2X , MZ . 3DE , 3X , MZ . 3DE 

,3X,IiZ.3DE,3X“ : I + l,Pc<I ),Ueg,Delt ,Uo( I ),Ho( I ),Hnus 
2190 ELSE 

2195 PRINT USING " 10X , 3D , 4X , 3D . DD , 3X , Z . DD , 4X , DD . 3D , 2X , MZ . 5DE , 3X ,MZ , 3DE 

,3X“ ; I + l ,Pc( I ),Ueg,Delt ,Uo< I ),Ho( I ) 

2200 END IF 

2205 NEXT I 

2210 IF lm=0 THEN 

2215 Okaccpt=l 

2220 BEEP 

2225 INPUT "OK TO ACCEPT THIS SET < 1=DEFAULT=YES , 0=NO )?" .Okaccpt 

2230 IF Okaccpt = l THEN OUTPUT ©Fi le ; Pc( * ) , Enf ( ♦ ) , Tp ( * ) 

2235 END IF 

2240 IF (Okaccpi = l OR Im=l ) AND Iout = l THEN 

2245 FOR 1=0 TO 2 

2250 OUTPUT ©Fout ;Ho( I ),Qdp( I ) 

2255 NEXT I 

2260 END IF 

2265 IF Ipi=0 then 

2270 Okrpt=l 

2275 BEEP 

2280 INPUT "UILL THERE BE ANOTHER RUN < 1=YES=DEFAULT ,0=NO )“ .Okrpt 

2285 IF Okrpt=l THEN 1683 

2290 ELSE 

2295 IF Jset<Nsets THEN 1683 

2300 END IF 

2305 ASSIGN ©File TO * 

2310 IF Iout=l THEN ASSIGN ©File TO • 

2315 SUBEND 

2320 DEF FNGrad(T) 

2325 Grad=-3. 8778575-5-2*4. 71 42857E-8«T 

2330 RETURN Grad 

2335 FNEND 
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2340 

2345 

2350 

2355 

2360 

2365 

2370 

2375 

2380 

2385 

2390 

2395 

2400 

2405 

2410 

2415 

2420 

2425 

2430 

2435 

2440 

2445 

2450 

2455 

2460 

2465 

2470 

2475 

2480 

2485 

2490 

2495 

2500 

2505 

2510 

2515 

2520 

2525 

2530 

2535 

2540 

2545 

2550 

2555 

2560 

2565 

2570 

2575 

2580 

2585 



DEF FNKcu(T) 

OFHC COPPER 250 TO 300 K 

Tk=T+273. 15 

K=434-. 112*Tk 

RETURN K 

FNEND 

DEF FNNueg(Tc,Egr ) 

RANGE OF UALIDITY: -20 TO 20 DE6 C 

Tk=Tc+273. 15 

Nul=7. 119G507E-3-Tk*( 7.4863347E-5-Tk*(2.6294943E-7-Tk*3,0833329E-10 ) ) 

Nu2=4.9237G38E-3-Tk*(4.9213912E-5-Tk*( l.G457534E-7-Tk*l .8335331E-10) ) 

Nu3=8.G586293E-3-Tk*(8.8837902E-5-Tk*(3.0495032E-7-Tk*3.499999GE-10) ) 

A2=(Nu3-2*Nu2+Nul >/200 

Al = (Nu2-Nul-940*A2 )/10 

A0=Nul-42*Al-1764*A2 

Nu=A0+Egr*( Al+Egr*A2 > 

RETURN Nu 
FNEND 

DEF FNCpeg( Tc , Egr ) 

RANGE OF UALIDITY: 0 TO 20 DEG C 

Tk=Tc+273.15 

Cpl = l .6701550E+3+Tk*G.3 

Cp2=l .4748125E+3+Tk*6.25 

Cp5=9 . 5800500E+2+Tk*7 . 3 

A2=( Cp3-2*Cp2+Cpl )/200 

Al = (Cp2-Cpl-900«A2 )/10 

A0=Cpl-40*Al-1600*A2 

Cp=A0+Egr*( Al+Egr*A2 ) 

RETURN Cp 
FNEND 

DEF FNRhoeg(T,Egr ) 

Rol=l .0G07093E+3-T*( 3.7031283E-l+T*4.0837183E-3) 
Ro2=1.0748272E+3-T*<4,426G195E-l+T*4.0939706E-5) 

Ro3=l .0885934E+3-T*(5.7355G53E-l+T*6.1281405E-3) 

A2=(Ro3-2*Ro2+Rol )/200 
Al = (Ro2-Rol-900*A2 )/10 
A0=Rol-40*Al-1600*A2 
Ro=A0+Egr*< Al+Egr*A2 ) 

RETURN Ro 
FNEND 

DEF FNPreg( T , Egr ) 

Pr=FNCpeg( T , Egr )*FNNueg(T,Egr )*FNRoeg( T , Egr )/FNKeg< T , Egr > 

RETURN Pr 
FNEND 

DEF FNKeg( Tc , Egr ) 

RANGE OF UALIDITY: -20 TO 20 DEG C 

Tk=Tc+273.15 

Kl=2.2824708E-l+Tk*<5.598928GE-4+Tk*5.571428GE-7 ) 

K2=2 . 5846616E-1+Tk * ( 2 . 3978571E-4+Tk*7 . 1428571E-7 ) 
K3=3,2138932E-l-Tk*(3.0042857E-4-TI *1 .4285714E-6 ) 



188 



2590 

2595 

2600 

2605 

2610 

2615 

2620 

2625 

2650 

2635 

2640 

2645 

2650 

2655 

2660 

2665 

2670 

2675 

2680 

2685 

2690 

2695 

2700 

2705 

2710 

2715 

2720 

2725 

2730 

2735 

2740 

2745 

2750 

2755 

2760 

2765 

2770 

2775 

2850 

2855 

2860 

2865 

2870 

2875 

2880 

2885 

2890 

2895 

2900 

2905 



A2 = <K3-2*K2+K1 )/200 
A1 = (K2-K1-900*A2 )/10 
A0=K1-40*A1-1600*A2 
K=A0+Egr* < Al+Egr*A2 ) 

RETURN K 
F.NEND 

DEF FNTanh(X) 

P=EXP<X ) 

Q=l/P 

Tanh=<P-Q)/<P+Q ) 

RETURN Tanh 
FNEND 

DEF FNTvsv(U) 

COM /Cc/ C<7) 

T=C<0) 

FOR 1=1 TO 7 

T=T+C< I )*U''I 
NEXT I 
RETURN T 
FNEND 

DEF FNBeta(T) 

Rop=FNRho(T+.l ) 

Rom=FNRho(T-,l ) 

Beta=-2/<Rop+Rom)«(Rop-Rom)/.2 

RETURN Beta 

FNEND 

DEF FNPsat ( Tc ) 

0 TO 80 deg F CURUE FIT OF Psat 
Tf=l .8*Tc+32 

Pa=5.945525+Tf«< , 15352082+Tf *< 1 .4840963E-3+Tf*9,6150671E-6 ) ) 
Pg=Pa-14.7 

IF Pg>0 THEN ! +=PSIG,-=in Hg 

Psat=Pg 
ELSE 

Psat=Pg*29.92/14.7 
END IF 
RETURN Psat 
FNEND 

SUB Xyread 
BEEP 

INPUT "ENTER FILE NAME",File$ 

BEEP 

INPUT "ENTER NUMBER OF X,Y PAIRS", N • 

ASSIGN ©File TO File$ 

FOR 1=1 TO N 

ENTER @File:X,Y 
PRINT X,Y 
NEXT I 
SUBEND 
SUB Purge 
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2910 

2915 

2920 

2925 

2930 

2935 

2940 

2945 

2950 

2955 

2960 

2965 

2970 

2975 

2976 

2977 

2980 

2985 

2990 

2995 

3000 

3005 

3010 

3015 

3020 

3025 

3030 

3035 

3040 

3045 

3046 

3050 

3051 

3052 

3053 

3054 

3055 

3056 

3057 

3058 

3059 

3060 

3062 

3063 

3065 

3070 

3075 

3080 

3085 

3090 



BEEP 

INPUT "ENTER FILE NAME TO BE DELETED" ,Fi le$ 

PURGE Files 
GOTO 2910 
SUBEND 
SUB Wilson 

COM /Will/ Doa( 4 ) , Dia( 4 ) , Kna< 4 ) , lact , Drool < 4 ) 

COM /Wil2/ Del t a , I sat , Nset s , Hod ,Cia( 3 , 3 ) , Alpaa< 3 ) 

DIM Emf<17),T(17),Xa(100),Ya(100),Pc<4),Tp<4) 

BEEP 

INPUT "PLEASE RE-ENTER NAME OF FILE",File$ 

ASSIGN ©File TO FileS 

INPUT "0=SMOOTH COPPER TUBE, 1=CUNI FIN, 2=KORODENSE, 3=FIN TITAN", Itube 
BEEP 

INPUT "0=NO INSERT, 1= TWISTED TAPE INSERT, 2=HITRAN1 , 3=HITRAN2" , Insert 
BEEP 

INPUT "GIUE A NAME FOR XY FILE",Xy$ 

CREATE BOAT Xy$,20 
ASSIGN @Xy TO XyS 
L=1 ,2192 
Do=Doa< Itube ) 

Di=Dia( Itube ) 

Km=Kma< Itube ) 

Ax=PI*Di ''2/4 ! Cross-sectional area 

Ao=PI*Do*L 

Rn=Do*L0G< Do/Di )/(2*Kn) 



Ini 


tial va 


L 


ues 


Tf= 


Tsat 






Alp 


a=.G55 






IF 


Insert= 


0 


THEN 


Ci = 


1,0 






END 


IF 






IF 


Insert= 


1 


THEN 


Ci = 


5.172 






END 


IF 






IF 


Insert= 


2 


THEN 


Ci = 


.226 






END 


IF 






IF 


Insert= 


2 


THEN 


Ci = 


.063 






END 


IF 






G=9 


.81 






lbeg=0 






lend=0 'CHANGE ' 


IF 


I act <10 




THEN 



Ibeg=Iact 
Iend=I act 
END IF 



IF FIUE TUBES IN BUNDLE 
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3095 

3100 

3105 

3110 

3115 

3120 

3125 

3130 

3135 

3140 

3145 

3150 

3155 

3160 

3161 
3165 
3170 
3175 
3180 
3185 

3190 ! 

3191 
3195 

3200 

3201 ! 

3202 ! 

3204 ! 

3205 ! 

3206 ! 
3210 ! 
3215 
3220 
3225 
3230 
3235 
3240 
3245 ! 
3250 
3255 
3260 
3265 . 
3270 
3275 
3280 

3285 

3286 

3287 

3288 

3289 

3290 



FOR I=Ibeg TO lend 
Sx=0 
5y=0 
Sx5=0 
Sxy=0 
Jset=0 

ASSIGN @File TO File$ 

ENTER OFi le;Dtg$, I lube, Egret ,Ddl ,Dd2 ,Dd3,Dd4,Dd5 
ENTER 0File:Pc(*),Emf<*),Tp<*) 

FOR J=0 TO 17 

T( J )=FNTvsv(Ernf ( J ) ) 

NEXT J 

Tvap=<T(0 )+T( 1 )+T<2 ) )/3 
Tliq = (T(3)+T<4) )/2 
Tliql=T(5 ) 

IF l5at=0 THEN 
Tsa t =T1 iq 1 
ELSE 

T5at=Tvap 
END IF 

6rad=FN6rad(T< 1+6 ) ) 

Tout = (T<2*I + 10)+T( 2*1 + 11 ) )/2 
Delt = Tout-T( 1+6 ) 

Tavg=T< 1+6 )+Delt*.5 
IF I=Iend THEN 
Delt=T<10)-T(9) 

Tavg=(T< 10 )-T< 9 > )/2 
END IF 

Water/Ethylene Glycol Mixture Froperties 
Rhoeg=FNRhoeg< Tavg , Egrat ) 

Nueg=FNNueg( Tavg,Egrat ) 

Mueg=Nueg*Rhoeg 
Cpeg=FNCpeg( Tavg , Egrat ) 

Keg=FNKeg( Tavg, Egrat ) 

Preg=Cpeg*Mueg/Keg 

Mdot=FNFncal< I ,T( 1+6 ),Pc< I ) ) 

Oeg=Mdot / < Rhoeg*Ax ) 

Reeg=Ueg*Di /Nueg 

Re5=4*Mdot/( Mueg*< PI*Di-4*Del ta ) ) 

Qdot=Mdot *Cpeg*Delt 
Odp=Qdot/Ao 

Lntd=Delt/LOG( (Tsat-T<I+6 ) )/< Tsat -TU+6 )-Delt ) ) 

Uo=Qdp/Lmt d 

IF lnsert=0 THEN 

Bbl=.06G8*< Di/L )*Reeg*Preg 

Bb2=l + .04*( (Di/L)*Reeg*Preg)''.BG66 

0foega=3 . 66+< Bbl/Bb2 ) 

END IF 
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3291 

3293 

3294 

3295 

3296 

3297 

3298 

3299 

3300 

3302 

3303 

5305 

3310 

3315 

3320 

3325 

3330 

3335 

3340 

3345 

3350 

3355 

3360 

3365 

3370 

3375 

3380 

3385 

3390 

3395 

3400 

3405 

3410 

3415 

3420 

3425 

3430 

5435 

3440 

3445 

3450 

5455 

3460 

3465 

3470 

3475 

3480 

5485 

3490 

3495 



IF Inserl=l THEN 

Omega = < 1+5 . 484E-3*Preg'' . 7* ( Res/Hod ) " 1 , 25 I” .5 
END IF 

IF Insert=2 THEN 

Omega = ( Reeg'' . 65 )* < Preg'' . 46 ) 

END IF 

IF Inserl=3 THEN 

Omega=( Reeg'' .76 )*(Preg''.46 ) 

END IF 

R-114/R113 Properties 
Hfg=FNHfg(Tsat ) 

Kf=FNK(Tf ) 

Rhof=FNRho( T f ) 

Muf=FNMu<Tf ) 

F=(Kf''3*Rhof''2*G*Hfg/(Muf*Do*Qdp ) I" .33333 

Ho=Alpa*F 

Two=Tsat-Qdp/Ho 

Tf = T sat / 3+2* Two/ 3 

Y=( 1/Uo-Rm )*F 

X=Do*F/( Keg*0mega ) 

PRINT ”0MEGA=" : Omega; “F=“ ;F; "X=“ ;X: "Y=" ; Y 

Xa(Jset)=X ! INEFFICIENT (MODIFY LATER) 

Ya( Jset )=Y 

S.x=5x+X 

Sy=Sy+Y 

Sxs=Sxs+X*X 

Sxy=Sxy+X*Y 

Jset=Jset+l 

IF Jset<Nsets THEN 3155 
ASSIGN ©File TO * 

Slope=( Nset s*Sxy-Sx*Sy )/ ( Nset s*Sxs-Sx'‘2 ) 

Intcpt=(Sy-Slope*Sx )/Nsets 

Cic=l /Slope 

Alpac=l/Intcpt 

Cerr=ABS( <Ci-Cic )/Cic ) 

Aerr=ABS< <Alpac-Alpa)/Alpac ) 

IF Cerr>.001 OR Aerr>.001 THEN 
Alpa=(Alpa+Alpac )*.5 
Ci = ( Ci+Cic )* .5 

PRINT "CIC=" ;Cic: "ALPA=" : Alpa 
GOTO 3100 
END IF 
BEEP 
BEEP 

PRINTER IS 1 

PRINT "CIC=" ;Cic; "ALPA=" : Alpa 
Cia( Insert , I )=Cic 
Alpaa( I ) = Alpac 
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3500 

3505 

3510 

3515 

3520 

3525 

3530 

3535 

3540 

3545 

3550 

3555 

35B0 

35B5 

3570 

3575 

3580 

3585 

3590 

3595 

3B00 

3B05 

3B10 

3B15 

3B20 

3625 

3630 

3635 

3640 

3645 

3650 

3655 

3660 

3665 

3670 

3675 

3680 

3685 

3690 

3695 

3700 

3705 

3710 

3715 

3720 

3725 

3730 

3735 

3740 

3745 



PRINTER IS 701 
FOR J=0 TO Nsets-l 

OUTPUT @Xy;Xa< J ) , Ya( J ) 

NEXT J 

PRINTER IS 1 
NEXT I 

ASSIGN 0Xy TO * 

SUBEND 
SUB Nusselt 

COM /Nus/ I a n , Tsa t , Qdp , Hoc ,K f , Rhof , Hf g , Muf , Do , I tube 
Do= .0159 
Ho=1000 
IF lin=0 THEN 
BEEP 

INPUT "ENTER TSAT AND HEAT FLUX ", Tsat , Qdp 
END IF 

HfQ=FNHfg<Tsat ) 

Two=Tsat-Qdp/Ho 
Tf=Tsat/3+2*Two/3 
Kf=FNK<Tf ) 

Rhof =FNRho< T f ) 

Muf=FNMu<Tf ) 

Hoc=. 655* < Kf "3*Rhof ''2*9,8 l*Hfg/( Muf *Do*Qdp ))" .333333 
IF ABS< <Ho-Hoc )/Hoc )>.001 THEN 
Ho=( Ho+Hoc )* . 5 
GOTO 3585 
END IF 

IF lin=0 THEN PRINT "H0=":Hoc 

SUBEND 

SUB Plotl 

DIM Yaa<4) 

PRINTER IS 705 

Idv=l 

BEEP 

INPUT "OK TO USE DEFAULT UALUES < 1=DEF-Y , 0=N ) " , I dv 
IF Idv=l THEN 
Itn=2 
Xnin=l 
Xtnax=5 
Xstep=l 
Ytnin=0 
Ynax=2 . 0 
Ystep=.5 
ELSE 
BEEP 

INPUT "ENTER MINIMUM AND MAXIMUM X-UALUES" , Xmi n , Xtnax 
BEEP 

INPUT "ENTER MINIMUM AND MAXIMUM Y-UALUES " , Yni n , Ynax 
BEEP 

INPUT "ENTER STEP SIZE FOR X-AXIS" .Xstep 
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3750 

3755 

37E0 

3765 

3770 

3775 

3780 

3785 

3790 

3795 

3800 

3805 

3810 

3815 

5820 

3825 

3830 

3835 

3840 

3845 

3850 

3855 

3860 

3865 

3870 

3875 

3880 

3885 

3890 

3895 

3900 

3905 

3910 

3915 

3920 

3925 

3930 

3935 

3940 

3945 

3950 

3955 

3960 

3965 

3970 

3975 

3980 

3985 

3990 

3995 



BEEP 

INPUT "ENTER STEP SIZE FOR Y-AXIS" , Ystep 
END IF 

PRINT "IN;SP1;IP 2300,1800,8300,6800;" 

PRINT “SC 0, 100,0, 100;TL 2,0;" 

Sfx = 100/(Xrnax-Xtnin ) 

Sfy = 100/(Ynax-Yrun) 

PRINT "PU 0,0 PD" 

FOR Xa=Xnin TO Xnax STEP Xstep 
X = (Xa-Xrnln)*Sfx 
PRINT "PA":X,",0; XT; " 

NEXT Xa 

PRINT "PA 100,0:PU;" 

PRINT "PU PA 0,0 PD" 

FOR Ya=Ynin TO Ymax STEP Ystep 
Y=( Ya-Ymin )*Sfy 
PRINT "PA 0, " : Y, "YT" 

NEXT Ya 

PRINT "PA 0,100 TL 0 2" 

FOR Xa=Xnin TO Xtnax STEP Xstep 
X=(Xa-Xrun)*Sfx 
PRINT "PA";X,",100; XT” 

NEXT Xa 

PRINT "PA 100,100 PU PA 100,0 PD" 

FOR Ya=Ynin TO Ynax STEP Ystep 
Y=( Ya-Ytnin )*Sfy 
PRINT ”PD PA 100, “ ,Y, "YT" 

NEXT Ya 

PRINT "PA 100,100 PU" 

PRINT "PA 0,-2 SR 1,5,2" 

FOR Xa=Xnin TO Xnax STEP Xstep 
X=< Xa-Xnin )*Sf x 
PRINT "PA";X," ,0; “ 

PRINT “CP -2,-l;LB" :Xa; "" 

NEXT Xa 

PRINT “PU PA 0,0" 

FOR Ya=Ymin TO Ymax STEP Ystep 
Y=(Ya-Yiriin)*Sfy 
PRINT "PA 0,";Y,“" 

PRINT "CP -4,-.25;LB" ;Ya; "" 

NEXT Ya 
BEEP 

INPUT "SELECT MODE ( 0=HN/H1 , 1=HN( avg )/Hl )" , I sn 
I srn = I sn+ 1 
IF Idv=l THEN 

IF Isn=l THEN Ylabel$="HN/Hl" 

IF Isn=2 THEN Ylabel$="HN< avg )/Hl " 
Xlabel$="Tube Number" 

ELSE 

BEEP 
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4000 

4005 

4010 

4015 

4020 

4025 

4030 

4035 

4040 

4045 

4050 

4055 

4060 

4065 

4070 

4075 

4080 

4085 

4090 

4095 

4100 

4105 

4110 

4115 

4120 

4125 

4130 

4135 

4140 

4145 

4150 

4155 

4160 

4165 

4170 

4175 

4180 

4185 

4190 

4195 

4200 

4205 

4210 

4215 

4220 

4225 

4230 

4235 

4240 

4245 



INPUT "ENTER X-LABEL" ,.Xlabel$ 

BEEP 

INPUT "ENTER Y-LABEL" , Ylabel$ 

END IF 

PRINT "SR 1.5,2:PU PA 50,-10 CP " : -LEN( Xlabe 1$ )/2 : " 0 : LB " ; Xlabel $ ; " " 
PRINT "PA -11,50 CP 0, ";-LEN<Ylabel$)/2*5/B; "DI 0 , 1 ; LB" ; Ylabel $ ; " ' 
PRINT "CP 0,0 DI" 

Okp=l 

BEEP 

INPUT "WANT TO PLOT DATA FROM A FILE ( 1=DEF=Y , 0=N )? " , Okp 
IF Okp=l THEN 
BEEP 

INPUT "ENTER THE NAME OF THE DATA FILE",Dfile$ 

ASSI6N ©File TO DfileS 
BEEP 

INPUT "ENTER THE BEGINNING RUN NUMBER", Md 
BEEP 

INPUT "ENTER THE NUMBER OF X-Y PAIRS STORED" ,Nset5 
BEEP 

INPUT "SELECT A SYMBOL FOR THE PLOTTER < 1 = * , 2=+ , 3=c , 4=o , 5=" ) " , Sy 
PRINT "PU DI" 

IF Sy=l THEN PRINT "SM«" 

IF Sy=2 THEN PRINT "SM+" 

IF Sy=3 THEN PRINT "SMc" 

IF Sy=4 THEN PRINT "SMo" 

IF Sy=5 THEN PRINT "SM"" 

FOR 1=1 TO Nsets 
FOR J=0 TO 3 

ENTER @File:Yaa< J ),D 
IF J=0 THEN Ytop=Yaa(0) 

Yaa( J )=Yaa( J )/Ytop 
NEXT J 

FOR J=0 TO 3 

X=( J+l-Xmin)«Sfx 
Y=<Yaa< J )-Ynin)*Sfy 
PRINT "PA",X,Y,"" 

NEXT J 
NEXT I 
BEEP 

ASSIGN ©File TO • 

GOTO 4040 
END IF 
BEEP 

INPUT "LIKE TO PLOT THE NUSSELT RELATION < 1=Y , 0=N )■? " , Oknus 
PRINT "PU;SM" 

IF Oknu5=l THEN 

FOR Xa=Xnin TO Xnax STEP Xstep/50 
X = ( Xa-Xp)in )*Sf x 

IF Isn=l AND Xa>Xrnin THEN Ya=Xa'' . 75-( Xa-1 ) " . 75 
IF Isn=2 AND Xa>Xmin THEN Ya=XaM-,25) 
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4250 

4255 

42B0 

4265 

4270 

4275 

4280 

4285 

4290 

4295 

4300 

4305 

4310 

4315 

4320 

4325 

4330 

4335 

4340 

4345 

4350 

4355 

4360 

4365 

4370 

4375 

4380 

4385 

4390 

4395 

4400 

4405 

4410 

4415 

4420 

4425 

4430 

4435 

44^0 

4445 

4450 

4455 

4460 

4465 

4470 

4475 

4480 

4485 

4490 

4495 



IF Xa=Xnin THEN Ya=l 
Y = ( Ya-Ymin )*Sfy 
PRINT ”PA" ,X,Y, "PD" 

NEXT Xa 
BEEP 

PRINT "PU” 

INPUT “MOUE THE PEN TO LABEL THE NUSSELT LINE", Ok 
PRINT "LBNusseir’ 

END IF 

IF Ism=2 THEN 
BEEP 

INPUT "LIKE TO PLOT EXPTL CURUE ( 1 = Y , 0=N ) " , Okex 
Nq=0 

IF Okex=l THEN 
BEEP 

INPUT "ENTER THE EXPONENT", Ex 
FOR Xa=Xmin TO Xmax STEP Xstep/10 
Nq=Nq+l 
Ya=Xa"(-Ex ) 

X = ( Xa-Xmin )*5f x 
Y=(Ya-Ynin)*Sfy 
IF Nq MOD 2=0 THEN 
PRINT "PA" ,X,Y, "PD" 

ELSE 

PRINT "PA",X,Y, "PU" 

END IF 
NEXT Xa 
PRINT "PU" 

BEEP 

INPUT "MOUE PEN TO LABEL AND HIT ENTER", Ok 
PRINT "LBs=0" 

PRINT "PR -1 0" 

PRINT "LB"; Ex;"" 

GOTO 4300 
END IF 
END IF 
GOTO 4510 
BEEP 

INPUT "LIKE TO PLOT KERN RELATIONSHIP <1=Y , 0=N )? " , Yes 
IF Yes=l THEN 

FOR Xa=Xt«)in TO Xnax STEP Xstep/20 
Ya=Xa"(-l/6 ) 

X = (Xa-Xi«iin )*Sfx 
Y=(Ya-Ymin)*Sfy 
PRINT "PA" ,X,Y, "PD" 

NEXT Xa 
PRINT "PU" 

BEEP 

INPUT "MOUE THE PEN TO LABEL KERN RELATIONSHIP" ,0k 
PRINT "LBKern;PU" 
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4B00 

4505 

4510 

4515 

4520 

4525 

4530 

4535 

4540 

4545 

4550 

4555 

4560 

4565 

4570 

4575 

4580 

4585 

4590 

4595 

4600 

4605 

4610 

4615 

4620 

4625 

4630 

4635 

4640 

4645 

4650 

4655 

4660 

4665 

4670 

4675 

4680 

4685 

4690 

4695 

4700 

4705 

4710 

4715 

4720 

4725 

4730 

4735 

4740 

4745 



END IF 

PRINT "PU PA 0,0" 

BEEP 

INPUT "LIKE TO PLOT EISSENBERG RELATION ( 1=Y,0=N )?" ,Okei 
IF Okei=l THEN 

FOR Xa=Xmin TO Xnax STEP Xstep/10 
Ya=,B+.42*Xa"(-.25) 

X=<Xa-Xrun)*Sf>; 

Y=( Ya-Ynin )*Sfy 
PRINT "PA" ,X,Y, "PD" 

NEXT Xa 
PRINT "PU" 

BEEP 

INPUT "MOUE THE PEN TO LABEL THE EISSENBERG LINE", Ok 
PRINT "LBEissenberg;PU“ 

END IF 

PRINT "PU SP0" 

SUBEND 

DEF FNPvst(Tc ) 

COM /Fid/ Ift 
DIM K(8) 

IF lfl=0 THEN 
BEEP 

PRINT "POST CORRELATION NOT AVAILABLE FOR R-114" 

STOP 
END IF 

IF Ifl=l THEN 

DATA -7.691234564,-26.08023696,-168.1706546,64.23285504,-118.9646225 
DATA 4.16711732,20.9750676,1.69,6 
READ K(») 

T=<Tc+273.15)/647.3 

Sum=0 

FOR N=0 TO 4 

Sun=Sun+K(N)*< 1-T )''(N+1 ) 

NEXT N 

Br=Sum/(T*( 1+K < 5 )* ( 1-T )+K( 6 )* ( 1-T >''2 ) )-( 1-T )/( K ( 7 )* ( 1-T )"2+K(8 ) ) 
Pr=EXP<Br) 

P=22120000*Pr 
END IF 

IF Ift=2 THEN 

Tf=Tc»l .8+32+459.6 

P=10-'(33. 0655-4330. 98/Tf-9.2635*LGT(Tf )+2 . 0539E-3*Tf ) 
P=P*101525/14.696 
END IF 

IF Ift=3 THEN 

A=9. 394685-3066. l/(Tc+273. 15) 

P=133.32«10'-A 
END IF 
RETURN P 
FNEND 
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4750 

4755 

4760 

4765 

4770 

4775 

4780 

4785 

4790 

4795 

4800 

4805 

4810 

4815 

48161 

4820 

4825 

4830 

4855 

4840 

4845 

4850 

4855 

4860 

4865 

4870 

4875 

4880 

4885 

4890 

4895 

4900 

4905 

4910 

4911 
4913 
4915 
4920 
4925 
4930 
4935 
4940 
4945 
4950 
4955 
4960 
4965 
4970 
4975 
4980 



DEF FNHfg(T) 

COH /Fid/ Ift 
IF lft=0 THEN 
Tf=T*l .8+32 

Hfg=6, 1451558E+l-Tf*(6.951079E-2+Tf*( 1 . 3988688E-4+1 . 9607843E-7*Tf > > 
Hfg=Hfg*2326 
END IF 

IF Ift=l THEN 

Hfg=2477200-2450*<T-10 ) 

END IF 

IF Ift=2 THEN 
Tf=T*l .8+32 

Hfg=7.0557857E+l-Tf*(4.838052E-2+l .2619048E-4*Tf ) 

Hfg=Hfg*2326.1 
Hfg=< 1 .611 -.0031 *T >*1 .0E+5 
END IF 

IF Ift=3 THEN 
Tk=T+273. 15 

Hf g=l . 35264E+6-TI< *( 6 . 38263E+2+Tk* . 747462 ) 

END IF 
RETURN Hfg 
FNEND 

— DEF FNMu<T) 

COM /Fid/ Ift 
IF lft=0 THEN 
Tk=T+273. 15 

Mu=EXP( -4. 4636+101 1 .47/Tk )*1 .E-3 
END IF 

IF Ift=l THEN 

A=247.8/(T+133.15) 

Mu=2.4E-5*10''A 
END IF 

IF Ift=2 THEN 

Mu=8.9629819E-4-T*( 1. 1094609E-5-T*5.566829E-8 > 

Tk=T+273. 15 

Mu = l . 345-5*10.0'- (503.0/ (Tk-2. 15 ) ) 

END IF 

IF Ift=3 THEN 

Tk = l/(T+273. 15) 

Mu=EXP(-ll .0179+Tk*( 1 .744E+3-Tk*(2.80335E+5-Tk*l . 1266 lE+8 ) > ) 

END IF 
RETURN Mu 
FNEND 

DEF FNUvsUTt ) 

COM /Fid/ Ift 
IF lft=0 THEN 
BEEP 

PRINT "UUST CORRELATION NOT AUAILABLE FOR R-114" 

STOP 
END IF 
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4985 

4990 

4995 

5000 

5005 

5010 

5015 

5020 

5025 

5030 

5035 

5040 

5045 

5050 

5055 

5060 

5065 

5070 

5075 

5080 

5085 

5090 

5095 

5100 

5105 

5110 

5115 

5120 

5125 

5130 

5135 

5140 

5145! 

5146 

5150 

5155 

5160 

5165 

5170 

5175 

5180 

5185 

5190 

5195 

5200 

5205 

5210 

5215 

5220 

5225 



IF Ift=l THEN 
P=FNPvst<Tl ) 

T=T1+273.15 

X=1500/T 

Fl = l/< l+T*l.E-4 ) 

F2=(1-EXP<-X ))'‘2.5*EXP(X )/X''.5 
B=.0015*F1-.000942*F2-.0004882*X 
K=2*P/(461 .52*T ) 

U=<1+(1+2*B*K)“.5)/K 
END IF 

IF Ift=2 THEN 
Tf=Tl*l .8+32 

U=13.955357-Tf *( . 16127262-Tf *5. 1726190E-4 ) 

V=IJ/16.018 
END IF 

IF Ifl=3 THEN 
Tk=Tt+273. 15 
P=FNPv5t<Tt ) 

'.^=133.95*Tk/P 
END IF 
RETURN U 
FNEND 

DEF FNCp(T) 

COM /Fid/ Ift 
IF lfl=0 THEN 
Tk=T+273.15 

Cp=.40118+Tk*(l .65007E-3+Tk*( 1.51494E-6-Tk*6.67853E-10) ) 

END IF 

IF Ift=l THEN 

Cp=4.21120858-T»(2.26826E-3-T*(4.42361E-5+2.71428E-7*T) ) 

END IF 

IF Ift=2 THEN 

Cp=9.2507273E-l+T*( 9. 3400433E-4+1 .7207792E-6*T ) 
Cp=(929.0+1.03*T)/1000. 

END IF 

IF Ift=5 THEN 
Tk=T+273.l5 

Cp=4. 1868*( 1 .6884E-2+Tk*< 3 . 35083E-3-Tk*( 7 .224E-6-Tk*7.61748E-9 ) ) ) 
END IF 

RETURN Cp*1000 
FNEND 

DEF FNP.ho<T) 

CON /Fid/ Ifl 
IF lft=0 THEN 
Tk=T+273,15 
X=l-( 1 .8*Tk/753.95) 

Ro=36,32+61 . 146414*X''( 1/3 )+16 . 41 8015*X + 17 . 476838*X'' . 5+1 . 119828*X''2 
Ro=Ro/. 062428 
END IF 

IF Ifi=l THEN 



199 



5230 

5235 

5240 

5245 

524B! 

5250 

5255 

52G0 

5265 

5270 

5275 

5280 

5285 

5290 

5295 

5300 

5305 

5310 

5315 

5320 

5325 

5330 

5335 

5340 

5345 

5350! 

5351 1 

5352 I 

5353 

5354 

5355 
5357 
5360 
5365 
5370 
5375 
5380 
5385 
5390 
5395 
5400 
5405 
5410 
5415 
5420 
5425 
5430 
5435 
5440 
5445 



Ro=999.52946+T*< . 01269-T* < 5 . 482513E-3-T* 1 , 234147E-5 ) ) 

END IF 

IF Ift=2 THEN 

Ro= 1.G207479E+3-T«(2 . 2186346+1*2.357829 lE-3 ) 

Ro=l./<< .G17+.000B47*T"1, 1 )*1.0E-3) 

END IF 

IF Ift=3 THEN 

Tk=T+273. 15-338.15 

Uf=9.24848E-4+Tk»(G.279GE-7+Tk*< 9.2444E-10+Tk*3.057E-12 ) ) 
Ro=l/Vf 
END IF 
RETURN Ro 
FNEND 

DEF FNPr(T ) 

Pr=FNCp<T)*FNMu<T )/FNK<T ) 

RETURN Pr 
FNEND 

DEF FNK<T) 

COM /Fid/ Ift 

IF lfl=0 THEN K=.071-.000261*T 
IF Ift=l THEN 

X=<T+273.15)/273.15 

K=-.92247+X*<2.8395-X«< 1.8007-X*< . 52577- .07344*X ) ) ) 

END IF 

IF Ift=2 THEN 

K=8.2095238E-2-T*<2.2214286E-4+T*2.3809524E-8 ) 

K=.0802-.000203*T 

K=.071-.0002G1*T 

Tf=T*1.8+32 

K=4.8461905E-2-Tf*G.5714286E-5 
K=K*1 .7308 
END IF 

IF Ifl=3 THEN 
Tk=T+273. 15 

K=4. 1868E-4*<519.442+.320920*Tk ) 

END IF 
RETURN K 
FNEND 

DEF FNHf(T) 

COM /Fid/ Ift 
IF lft=0 THEN 
BEEP 

PRINT “HF CORRELATION NOT FOR R-114" 

STOP 
END IF 

IF Ift=l THEN 

Hf=T*<4.203849-T*<5.88132E-4-T*4.55160317E-6 ) ) 

END IF 

IF Ift=2 THEN 
Tf=T*l .8+32 
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5450 

5455 

54B0 

5455 

5470 

5475 

5480 

5485 

5490 

5495 

5500 

5505 

5510 

5515 

5520 

5525 

5530 

5535 

5540 

5545 

5550 

5555 

5550 

5555 

5570 

5575 

5580 

5585 

5590 

5595 

5800 

5505 

5510 

5615 

5520 

5525 

5630 

5635 

5640 

5645 

5650 

5655 

5660 

5665 

5670 

5675 

5680 

5685 

5690 

5695 



Hf =8. 2078571 +Tf*< . 19467857+Tf*l .3214286E-4 ) 
Hf=Hf*2.326 
END IF 

IF Ifl=3 THEN 

Hf=250 ! TO BE UERIFIED 

END IF 

RETURN Hf*1000 

FNEND 

SUB Plot2 

COM /Drl/ Star , Syn , Icon 

COM /Fid/ Ift 

DIM C<9),Xya(7),Doa(4) 

DATA 0.0158,0.0158,0.0158,0.0158 
READ Doa( * ) 

Fw=l 

PRINTER IS 1 
BEEP 

PRINT USING "4X," "Select Option X-Y Limits: 

PRINT USING “6X,“"0 Use default values 

PRINT USING "6X,""1 Use new values 

INPUT Okd 
BEEP 

INPUT "ENTER TUBE CODE",Icode 
Do=Doa< Itube > 

Iht=2 

BEEP 

PRINT USING "4X," "Select option: 

PRINT USING "6X,""0 h versus q 

PRINT USING "BX,""1 q versus Delta-T 

PRINT USING "6X,""2 h versus Delta-T (default)" 
INPUT Iht 
PRINTER IS 705 

IF Okd=0 THEN lAXIS DEFAULT UALUES 
IF lht=0 THEN ! ( h vs q ) 

Ynin=0 
YMax=60 
Ystep=10 
Xmin= .2 
Xmax=l . 4 
Xstep=.2 
END IF 

IF Iht = l THEN !(q vs t ) 

Xnin=0 
Ymin=0 
Ynax=.5 
XMax=15 
Xst ep=3 
Ystep=. 1 
END IF 

IF Iht=2 THEN ! ( h vs t ) 
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5700 

5705 

5710 

5715 

5720 

5725 

5730 

5735 

5740 

5745 

5750 

5755 

5760 

5765 

5770 

5775 

5780 

5785 

5790 

5795 

5800 

5805 

5810 

5815 

5820 

5825 

5830 

5835 

5840 

5845 

5850 

5855 

5860 

5865 

5870 

5875 

5880 

5885 

5890 

5895 

5900 

5905 

5910 

5915 

5920 

5925 

5950 

5935 

5940 

5945 



Xmin=0 

Ymin=0 



Xmax=50 
Ymax=6 
X5tep=10 
Ystep=l 
END IF 
END IF 

IF Okd=l THEN 
BEEP 

INPUT “ENTER MINIMUM AND MAXIMUM X-UALUES" ,Xnin,Xmax 
BEEP 

INPUT “ENTER MINIMUM AND MAXIMUM Y-UALUES" ,Ymin,Ymax 
BEEP 

INPUT “ENTER STEP SIZE FOR X-AXIS" ,Xstep 
BEEP 

INPUT "ENTER STEP SIZE FOR Y-AXIS" , Ystep 
END IF 
BEEP 

PRINT “IN;SP1;IP 2300,1800,8300,6800;“ 

PRINT "SC 0, 100,0, 100;TL 2,0;" 

Sfx=100/<Xnax-Xnln) 



Sfy=100/<Ynax-Ymin) 

BEEP 



lcg=0 

INPUT “LIKE TO BY-PASS CAGE < 1=Y , 0=N=DEFAULT )? “ , leg 
IF Icg=l THEN 6175 
PRINT "PU 0,0 PD“ 

FOR Xa=Xnin TO Xmax STEP Xsiep 
X=<Xa-Xmin)*Sfx 
PRINT "PA" ;X, “ ,0; XT: " 

NEXT Xa 

PRINT "PA 100,0;PU; " 

PRINT "PU PA 0,0 PD" 

FOR Ya=Ytnin TO Ymax STEP Ystep 
Y=<Ya-Ymin)*Sfy 
PRINT "PA 0, "; Y, "YT" 

NEXT Ya 

PRINT "PA 0,100 TL 0 2" 

FOR Xa=Xmin TO Xmax STEP Xstep 
X=(Xa-Xmin>*Sfx 
PRINT "PA" :X , " , 100; XT" 

NEXT Xa 

PRINT "PA 100,100 PU PA 100,0 PD" 

FOR Ye=Ymin TO Ymax STEP Ystep 
Y=(Ya-Ymin)«Sfy 
PRINT "PD PA 100,",Y,"YT" 

NEXT Ya 

PRINT "PA 100,100 PU" 

PRINT "PA 0,-2 SR 1,5,2“ 
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B950 
5955 
59B0 
59B5 
5970 
5975 
5980 
5985 
5990 
5995 
B000 
6005 
B010 
6015 
6020 
6025 
6030 
6035 
6040 
6045 
6050 
6055 
6060 
6065 
6070 
6075 
6080 
6085 
6090 
6095 
6100 
6105 
kW/n 
6110 
6115 
6120 
R 1.1 
6125 
6130 
6135 
1.5,2 
6140 
6145 
2,0 
6150 
6155 
6160 
6165 
6170 
6175 



FOR Xa=Xnin TO Xnax STEP Xslep 
X=( Xa-Xnin )*Sfx 
PRINT "PA" ;X, •' ,0; " 

IF Xa<l AND Xa<>0 THEN PRINT "CP -1 . 5 , -1 ; LB0 ; PR -1,0;LB 
IF Xa=0 THEN PRINT "CP -.5,-1 ;LB0" . 

Xin=0 

IF Xa MOD 1=0 THEN Xin=l 

IF Xa>=10 THEN PRINT "CP -2 ,-l : LB " : Xa : " " 

IF Xa>l AND Xa<10 AND Xin=l THEN PRINT "CP -1.25,-1;LB" 
IF Xa>l AND Xin=0 THEN PRINT "CP -2 , -1 ; LB" : Xa i " " 

IF Xa=l THEN PRINT "CP -1,-1:LB1.0" 

NEXT Xa 

PRINT "PU PA 0,0" 

FOR Ya=Ynin TO Ymax STEP Ystep 
Y=< Ya-Ymin)*Sfy 
PRINT "PA 0,":Y,"" 

IF lhl=0 AND Ya>0 THEN PRINT "PR 2,0" 

IF Ya<l AND Ya<>0 THEN PRINT "CP -4 ,- . 25 ; LB0: PR -2,0;LB 

IF Ya=0 THEN PRINT "CP -2,-.25;LB0" 

IF Ya>l AND Iht<2 THEN PRINT "CP -5 ,- . 25 ; LB " ; Ya : “ " 

IF Ya>9 AND Ihl=2 THEN PRINT "CP -4 , - . 25 ; LB " ; Ya ; " " 

IF Ya>0 AND Ya<10 AND Ihl=2 THEN PRINT "CP -3,-.25:LB"! 

! IF Ya=l THEN PRINT "CP -4 , - . 25 ; LBl . 0 " 

NEXT Ya 
IF Okd=l THEN 
BEEP 

INPUT "ENTER X-LABEL" , XlabelS 
BEEP 

INPUT "ENTER Y-LABEL" , YlabelS 
END IF 

IF IhtOl THEN 

PRINT "SR 1.5,2;PU PA -12,35:DI 0,l;LBh;PR l,0.5:LBo:PR 

PRINT "PR -1,0.5:SR 1 , 1 .5 ; LB2 ; SR 1.5,2;PR .5,.5;LB.;PR . 
ELSE 

PRINT "PA -12,39:DI 0 , 1 ; LBq/( MW/m ; SR 1,1.5;PR -1,0.5;LB2 
,0.5: LB)" 

END IF 

IF lht=0 THEN 

PRINT "SR 1.5,2;PU PA 40 , -10 ; DI ; LBq/( MW/m ; PR 0.5,1;SR1, 
:PR .5,-1; LB)" 

ELSE 

PRINT "DI PA 38,-10;LB(T:PR .5,-l;LB5;PR .5,1;LB-T;PR -2 
PU;PR .5,-4;" 

PRINT "LBwo;PR .5,1:LB)/K" 

END IF 

PRINT "CP 0,0 DI" 

Xlg=l.E+6 
Xug=-1 .E+6 
Xal=50 



";Xa; 

;Xa: “ " 

";Ya; "" 

Y a ; " " 

-1 ,0.5;LB/< 
5,0:LBK >" 
;SR 1.5,2:P 

1 .5;LB2:SR 
.4,3 PD PR 
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G180 

6185 

6190 

6195 

6200 

6205 

6210 

6215 

6220 

6225 

6230 

6235 

6240 

6245 

6250 

6255 

6260 

6265 

6270 

6275 

6280 

6285 

6290 

6295 

6300 

6305 

6510 

6315 

6320 

6325 

6530 

6335 

6340 

6345 

6350 

6355 

6560 

6365 

6370 

6375 

6380 

6385 

6390 

6395 

6400 

6405 

6410 

6415 

6420 

6425 



Yal=95 

Nrun=0 

BEEP 

INPUT "WANT TO PLOT DATA FROM A FILE < 1=Y , 0=N )? " , Ok 
Xll = l ,E+6 
Xul=-1 .E+6 
Okp=0 

IF Ok=l THEN 
BEEP 

INPUT "ENTER THE NAME OF THE PLOT DATA FILE" ,D_f i le$ 

ASSI6N ©File TO D_file$ 

IF lcomb<>0 THEN 6265 

Sx=0 

Sy=0 

6x2=0 

Sxy=0 

Md=l 

BEEP 

INPUT "ENTER THE BEGINNING RUN NUMBER (DEF=l)",Md 

Npairs=9 

BEEP 

INPUT "ENTER THE NUMBER OF X-Y PAIRS STORED ( DEF=9 ) " , Npairs 

Nrun=Nrun+Npairs 

PRINTER IS 1 

BEEP 



PRINT 


USING 


"4X, " 


"Select a symbol * 


PRINT 


USING 


"4X," 


" 1 


Star 2 Plus sign 


PRINT 


USING 


"4X, " 


II - 


Circle 4 Square 


PRINT 


USING 


"4X," 


" 5 


Rombus 


PRINT 


USING 


"4X," 


" 6 


Right-side-up triangle" " 


PRINT 

INPUT 

BEEP 


USING 

Syn 


"4X." 


" 7 


Up“Side-down triangle""" 



INPUT “ENTER TUBE NUMBER FOR PLOTTING ( 0=TOP , l=SECOND I tube 

PRINTER IS 705 

IF Syn=l THEN PRINT "5M*" 

IF Sy(n=2 THEN PRINT "SM+" 

IF Syfn=3 THEN PRINT "SMo" 

IF Md>l THEN 

FOR 1=1 TO (Md-1 > 

ENTER ©File;XyaC ) 

NEXT I 
END IF 

FOR 1=1 TO Npairs 

ENTER ©File:Xya( * ) 

Ya=Xya( Itube*2 ) 

Xa=Xya( Itube*2+1 ) 

Yc=LOG(Xa > 

Xc=LOG(Xa/Ya ) 

Sx=Sx+Xc 
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6430 

6435 

6440 

6445 

6450 

6455 

6460 

6465 

6470 

6475 

6480 

6485 

6490 

6495 

6500 

6505 

6510 

6515 

6520 

6525 

6530 

6535 

6540 

6545 

6550 

6555 

6560 

6565 

6570 

6575 

6580 

6585 

6590 

6595 

6600 

6605 

6610 

6615 

6620 

6625 

6630 

6635 

6640 

6645 

6650 

6655 

6660 

6665 

6670 

6675 



Sy=5y+Yc 
Sx2=Sx2 + Xc''2 
Sxy=Sxy+Xc -»Yc 
IF lht=0 THEN 
Xt=Ya 
Ya=Ya/Xa 



Xa=Xt 

IF Xa/l.E+6>Xul THEN Xul=Xa/l.E+6 
IF Xa/1.E+6<X11 THEN Xll=Xa/l,E+6 
END IF 

IF lht=0 THEN 

X=<Xa*l.E-6-Xmin)*Sfx 
Y=<Ya*l.E-3-Ymin>*5fy 
END IF 

IF Iht=l THEN 

X = ( Xa-Xmin >«Sfx 
Y=< Ya*l ,E-6-Ymin )*Sfy 
END IF 



IF Iht=2 THEN 

X=(Xa/Ya-Xmn )*Sfx 



Y=< Ya*l .E-3-Ymin)*Sfy 
END IF 



IF Y>100 OR Y<0 THEN 6585 



IF Sym>3 THEN PRINT 
IF Sym<4 THEN PRINT 
PRINT ''PA",X,Y,"“ 

IF Sym>3 THEN PRINT 
IF Syn=4 THEN PRINT 
IF Sym=5 THEN PRINT 
IF Syn=6 THEN PRINT 
IF Sym=7 THEN PRINT 
NEXT I 



"SM" 

"SR 1,4, 2. 4" 

“SR 1.2, 1.6" 

” UC2 , 4 , 99 , 0 , -8 , -4 , 0 , 0 , 8 , 4 , 0 ; " 
"UC3, 0,99, -3, -6, -3, 6, 3, 6, 3,-6; 
"UC0,5.3,99,3,-8,-6,0,3,8; ■’ 

" UC0 , -5 , 3 , 99 , -3 , 8 , 6 , 0 , -3 , -8 ; " 



BEEP 



INPUT "WANT TO LABEL < 1=Y , 0=N )? " , I Ibl 
IF Ilbl=l THEN 

IF Sym>3 THEN PRINT "SM” 

IF Sym<4 THEN PRINT "SR 1.4, 2. 4" 

PRINT "PA" ,Xal,Yal , "“ 

IF Syn>3 THEN PRINT "SR 1.2, 1,6" 

IF Sym=4 THEN PRINT "UC2 , 4 , 99 , 0 , -8 , -4 , 0 , 0 , 8 , 4 , 0 : " 
IF Sym=5 THEN PRINT "UC3 , 0 , 99 , -3 , -6 , -5 , 6 , 3 , 6 , 3 , -6 ; 
IF Sym=6 THEN PRINT “UC0 , 5 . 3 , 99 , 3 , -8 , -6 , 0 , 3 , 8 ; " 

IF Syn=7 THEN PRINT " UC0 , -5 . 3 , 99, -3 , 8 , 6 , 0 , -3 , -8 ; " 
PRINT "SM" 



IF Sym<4 THEN PRINT "PR 2,0" 

PRINT "PR 2,-1.0:SR 1 . 0 , 1 . 8 ; LB" ; D_f i le$ ; " " 

Yal=Yal-5 

BEEP 

INPUT "WANT TO ADD ANOTHER STRING ( 1=Y,0=N )7" , las 
IF Ias=l THEN 
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G680 

6B85 

GB90 

GB95 

G700 

6705 

G710 

6715 

G720 

G725 

6730 

6735 

6740 

6745 

6750 

6755 

6760 

6765 

6770 

6775 

6780 

6785 

6790 

6795 

6800 

6805 

6810 

6815 

6820 

6825 

6830 

6835 

6840 

6845 

6850 

6855 

6860 

6865 

6870 

6875 

6880 

6885 

6890 

6895 

6900 

6905 

6910 

6915 

6920 

6925 



BEEP 

INPUT "ENTER THE STRING" , Labels 
PRINT "PR 2,0:SR 1 . 0 , 1 . 8 ; LB " ; Labels ; " " 

GOTO 6665 
END IF 
END IF 
BEEP 

INPUT "UANT TO COMBINE ANOTHER FILE? < 1 = Y , 0=N ) " , Icomb 
ASSIGN ©File TO * 

Xll=5 

Xul=45 

IF lcomb<>0 THEN 6220 
BEEP 

INPUT "WANT TO PLOT A LEAST-SQUARES LINE < 1=Y , 0=N ) " , 1 1 s 
IF Ils=l THEN 
BEEP 

INPUT "SELECT EXPONENT: 0=COMPUTE, 1=0.75 )", lexp 
BEEP 

INPUT “SELECT CURVE TYPE <0=SOLID, 1=DASHED )" , lit 

Ilt=Ilt+l 

PRINT "SM" 

IF lexp=0 THEN 

Bb = <Nru(n*Sxy-Sy»Sx )/< Nrun*Sx2-Sx''2 ) 

ELSE 

Bb=.75 
END IF 

Aa = (Sy-Bb*Sx )/Nrun 
Aa=EXP(Aa ) 

PRINTER IS 1 

PRINT USING "10X,""a = "",Z.4DE":Aa 
PRINT USING "10X,""n = "",Z.4DE";Bb 
PRINTER IS 705 
ln=0 

IF lht=0 THEN Xxstep=Xstep/40 
IF lht>0 THEN Xxstep=Xstep/ 10 
FOR Xa=Xll TO Xul STEP Xxstep 
IF Xa>.99*Xmex THEN 6995 
IF Iht=l THEN Ya=Aa*Xa''Bb 

IF lht=0 THEN Ya=Aa"< 1/Bb )*(Xa«l .E+6 )''< <Bb-l )/Bb ) 
IF Iht=2 THEN Ya=Aa*Xa''( Bb-1 ) 

IF lht=0 THEN 

Y=(Ye*l.E-3-Ymin)*Sfy 
X=(Xa-Xmin)»Sfx 
END IF 

IF Iht=l THEN 

Y=(Ya*l.E-6-Ymin)*Sfy 
X=<Xa-Xmin)*Sfx 
END IF 

IF Iht=2 THEN 
. Y=( Ya*l .E-3-Ymln >*Sfy 
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G930 

6935 

6940 

6945 

6950 

6955 

6960 

6965 

6970 

6975 

6980 

6985 

6990 

6995 

7000 

7005 

7010 

7015 

7020 

7025 

7030 

7035 

7040 

7045 

7050 

7055 

7060 

7065 

7070 

7075 

7080 

7085 

7090 

7095 

7100 

7105 

7110 

7115 

7120 

7125 

7130 

7135 

7140 

7145 

7150 

7155 

7160 

7165 

7170 

7175 



X = (Xa-Xmin >*Sfx 
END IF 

IF Y<0 THEN Y=0 
IF Y>100 THEN 6990 
IF 111=1 THEN 

PRINT "PA" ,X,Y, "PD" 

ELSE 

In=In+ 1 
Ir=In MOD 111 

IF Ir=l THEN PRINT "PA" , X , Y , “PD" 

IF lr=0 THEN PRINT "PA" , X , Y , "PU" 

END IF 
NEXT Xa 
PRINT "PU" 

END IF 
lcoi*ib=0 
60T0 6185 
END IF 

PRINT "PU SM" 

BEEP 

INPUT "UANT TO PLOT NUSSELT LINE < 1=Y , 0=N )? " , Inp 

IF lnp=0 THEN 7125 

BEEP 

INPUT "ENTER T5AT (DEFAULT=18 DE6 C>",Tsal 
HfQ=FNHfg<Tsat ) 

Xll=5 

Xul=45 

FOR Xe=Xll TO Xul STEP Xstep/50 
T f i ln=T sal -Xa* . 5 
Kf=FNK(Tf ilm) 

Rhof=FNRho<Tfilm) 

Muf=FNMu<Tf ilfi ) 

Ya=.728*(Kf "3*Rhof ''2*9.81*Hfg/(Muf*Do*Xa ) )".25 
X=(Xa-Xmin )*Sfx 
Y=<Ya*l.E-3-Ynin)*Sfy 
PRINT "PA",X,Y,"PD" 

NEXT Xa 

PRINT "PU PA 0,0" 

PRINT “PU PA 0,0 SP0" 

SUBEND 
SUB Plots 

COM /Drl/ Star , Sym , Icon 
COM /Fid/ Ift 
DIM C<9) 

Fw=l 

PRINTER IS 1 
BEEP 

PRINT USING "4X, ""Select Option X-Y Limits: 

PRINT USING "6X,""0 Use default values 

PRINT USING "GX,""1 Use new values 
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7180 

7185 

7190 

7195 

7200 

7205 

7210 

7215 

7220 

7225 

7230 

7235 

7240 

7245 

7250 

7255 

72B0 

7265 

7270 

7275 

7280 

7285 

7290 

7295 

7300 

7305 

7310 

7315 

7320 

7325 

7330 

7335 

7340 

7545 

7350 

7355 

7360 

7365 

7370 

7375 

7380 

7385 

7390 

7395 

7400 

7405 

7410 

7415 

7420 

7425 



INPUT Okd 
PRINTER IS 705 
IF Okd=0 THEN 
Xnin=0 
Ym in=0 
Xna?<=15 
Ynax=15 
Xstep=3 
Ystep=3 
ELSE 
BEEP 

INPUT "ENTER MINIMUM AND MAXIMUM X-UALUES" ,Xmin ,Xnsx 
BEEP 

INPUT "ENTER MINIMUM AND MAXIMUM Y-UALUES " , Ynin , Ynax 
BEEP 

INPUT "ENTER STEP SIZE FOR X-AXIS" ,Xstep 
BEEP 

INPUT "ENTER STEP SIZE FOR Y-AXIS" , Ystep 
END IF 
BEEP 

PRINT "IN;SP1:IP 2300,1800,8500,6800;" 

PRINT "SC 0, 100,0, 100;TL 2,0:" 

Sfx=100/(Xnax-Xnin ) 

Sfy=100/(Ynax-Ymin > 

BEEP 

lcg=0 

INPUT "LIKE TO BY-PASS CAGE (1=Y,0=N=DEFAULT )?" , leg 
IF Icg=l THEN 7620 
PRINT "PU 0,0 PD" 

FOR Xa=Xnin TO Xnax STEP Xstep 
X = ( Xa-Xnin )*Sfx 
PRINT "PA";X, ",0; XT; " 

NEXT Xa 

PRINT “PA 100,0;PU; " 

PRINT "PU PA 0,0 PD“ 

FOR Ya=Ynin TO Ynax STEP Ystep 
Y = ( Ya-Ynin )»Sf y 
PRINT "PA 0, ";Y, "YT" 

NEXT Ya 

PRINT "PA 0,100 TL 0 2" 

FOR Xa=Xnin TO Xnax STEP Xstep 
X=(Xa-Xnin )*Sfx 
PRINT "PA" :X , " , 100; XT" 

NEXT Xa 

PRINT "PA 100,100 PU PA 100,0 PD" 

FOR Ya=Ynin TO Ynax STEP Ystep 
Y=( Ya-Ynin )*Sfy 
PRINT "PD PA 100,",Y,"YT" 

NEXT Ya 

PRINT "PA 100,100 PU" 
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7430 

7435 

7440 

7445 

7450 

7455 

74B0 

74B5 

7470 

7475 

7480 

7485 

7490 

7495 

7500 

7505 

7510 

7515 

7520 

7525 

7530 

7535 

7540 

7545 

7550 

7555 

75B0 

7565 

7570 

7575 

7580 

7585 

7590 

7595 

7600 

7605 

7610 

7615 

7620 

7625 

7630 

7635 

7640 

7645 

7650 

7655 

7660 

7665 

7670 

7675 



PRINT "PA 0,-2 SR 1.5,2" 

FOR Xa=Xmn TO Xnax STEP Xstep 
X=< Xa-Xni n )*Sf >; 

PRINT "PA”;X,",0;" 

IF Xa<l AND Xa<>0 THEN PRINT "CP -1 . 5 , - 1 : LB0 ; PR -1 , 0 ; LB " : Xe ; " " 
IF Xa=0 THEN PRINT "CP -.5,-1 ;LB0" 

X 1 n=0 

IF Xa MOD 1=0 THEN Xin=l 

IF Xa>=10 THEN PRINT "CP -2 , -1 ; LB" : Xa : " " 

IF Xa>l AND Xa<10 AND Xin=l THEN PRINT "CP -1 .25,-1 ;LB" ; Xai " " 
IF Xa>l AND Xin=0 THEN PRINT "CP -2 , -1 ; LB" : Xa : " " 

IF Xa=l THEN PRINT "CP -1,-1:LB1.0" 

NEXT Xa 

PRINT "PU PA 0,0" 

Iht=2 ! MODIFY 
FOR Ya=Ymin TO Ynax STEP Ystep 
Y=( Ya-Ynin )*Sf y 
PRINT "PA 0,“;Y,"" 

IF lht=0 AND Ye>0 THEN PRINT "PR 2,0" 

IF Ya<l AND Ya<>0 THEN PRINT "CP -4 , - . 25 ; LB0 ; PR -2 , 0 ; LB " ; Ya : " " 
IF Ya=0 THEN PRINT "CP -2,-.25;LB0" 

IF Ya>l AND Iht<2 THEN PRINT "CP -5, - . 25 ; LB" ; Ya; " " 

IF Ya>9 AND Iht=2 THEN PRINT "CP -4 , - . 25 ; LB " ; Ya ; " " 

IF Ya>0 AND Ya<10 AND Iht=2 THEN PRINT "CP -3 , - . 25; LB" ; Ya ; " " 

IF Ya = l THEN PRINT "CP -4 , - . 25 : LBl . 0” 

NEXT Ya 
IF Okd=0 THEN 
XlabeI$="X" 

YIabeIf=“Y" 

ELSE 

BEEP 

INPUT "ENTER X-LABEL" ,Xlabel$ 

BEEP 

INPUT "ENTER Y-LABEL" , Ylabel$ 

END IF 

PRINT "SR 1.5,2 PU PA 50,-10 CP " : -LEN( Xlabel$ )/2 = "0 ; LB “ ; Xlabel$ ; " " 
PRINT "PA -11,50 CP 0, “;-LEN(Ylabel$)/2*5/6; "DI 0 , 1 ; LB" ; YlabelS : 
PRINT "CP 0,0 DI" 

Nrun=0 

BEEP 

INPUT "WANT TO PLOT DATA FROM A FILE (1 =Y , 0=N )? " , Ok 
Okp=0 

IF Ok=l THEN 
BEEP 

INPUT "ENTER THE NAME OF THE PLOT DATA FILE" ,D_f ile$ 

ASSIGN ©File TO D_file$ 

IF lconb<>0 THEN 7690 

Sx=0 

Sv=0 

S>.2=0 
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7G80 

7685 

7690 

7695 

7700 

7705 

7710 

7715 

7720 

7725 

7730 

7735 

7740 

7745 

7750 

7755 

7760 

7765 

7770 

7775 

7780 

7785 

7790 

7795 

7800 

7805 

7810 

7815 

7820 

7825 

7830 

7635 

7840 

7845 

7850 

7855 

7860 

7865 

7870 

7875 

7880 

7885 

7890 

7895 

7900 

7905 

7910 

7915 

7920 

7925 



Sxy=0 

Md=l 



BEEP 

INPUT "ENTER THE BEGINNING RUN NUMBER <DEF=l)",Md 
Npairs=9 



BEEP 

INPUT "ENTER THE NUMBER OF X-Y PAIRS STORED < DEF=9 ) " , Npa irs 
BEEP 

INPUT "SELECT TUBE NUMBER ( 0=TOP , l=SECOND )",Itube 

Nrun=Nrun+Npairs 



PRINTER IS 1 
BEEP 

PRINT USING "4X, ""Select a synbol = 



PRINT USING "4X,"" 1 
PRINT USING "4X,“" 3 
PRINT USING "4X,"" 5 
PRINT USING "4X,"" B 
PRINT USING "4X,"" 7 
INPUT Syn 
PRINTER IS 70S 
IF Syn=l THEN PRINT 
IF Sym=2 THEN PRINT 
IF Sytn=3 THEN PRINT 
Md=Itube»9 
IF Md>l THEN 

FOR 1=0 TO (Md-1 ) 



Star 2 Pius sign 

Circle 4 Square""" 
Rombus 

Right-side-up triangle" 
Up-side-down triangle"" 



SM«" 

SM+" 

SMo" 



II 



ENTER @File;X,Y 
NEXT I 
END IF 

FOR 1=1 TO Npairs 

ENTER @File;Xa,Ya 
Sx=Sx+Xa 



Sy=Sy+Ya 

Sx2=Sx2+Xa"2 



Sxy=Sxy+Xa«Ya 
X=( Xa-Xmin )*Sfx 



Y=<Ya-Yiiin)»Sfy 



IF Y>100 OR Y<0 THEN 7910 
IF Syn>3 THEN PRINT "SM" 

IF Syt«<4 THEN PRINT "SR 1.4, 2. 4" 

PRINT "PA",X,Y,"" 

IF Syn>3 THEN PRINT "SR 1.2, 1.6" 

IF Syw=4 THEN PRINT "UC2 , 4 , 99 , 0 , -8 , -4 , 0 , 0 , 8 , 4 , 0 : " 

IF Sym=5 THEN PRINT "UC3 , 0 , 99 , -3 , -6 , -3 , 6 , 3 , 6 , 3 , -6 ; " 
IF Sym=G THEN PRINT "UC0 ,5 . 3 , 99 , 3 ,-8 , -6 , 0 , 3 , 8 ; " 

IF Sym=7 THEN PRINT "UC0 ,-5 . 3 , 99 , -3, 8 , 6 , 0 , -3 , -8 ; " 
NEXT I 
BEEP 

INPUT "WANT TO LABEL < 1=Y , 0=N )? " , I Ibl 
IF Ilbl=l THEN 
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7930 

7935 

7940 

7945 

7950 

7955 

7960 

7965 

7970 

7975 

7980 

7985 

7990 

7995 

8000 

8005 

8010 

8015 

8020 

8025 

8030 

8035 

8040 

8045 

8050 

8055 

8060 

8065 

8070 

8075 

8080 

8085 

8090 

8095 

8100 

8105 

8110 

8115 

8120 

8125 

8130 

8135 

8140 

8145 

8150 

8155 

8160 

8165 

8170 

8175 



IF Syn>3 THEN PRINT "SM" 

IF Syn<4 THEN PRINT "SR 1.4, 2. 4" 

PRINT "PA" ,Xal , Yal , 

IF Syn>3 THEN PRINT "SR 1.2, 1.6" 

IF Sym=4 THEN PRINT "UC2 , 4 , 99 , 0 , -8 , -4 , 0 , 0 , 8 , 4 , 0 ; " 

IF Syrn=5 THEN PRINT "UC3 , 0 , 99 , -3 , -6 , -3 , 6 , 3 , 6 , 3 , -6 ; " 

IF Syn=6 THEN PRINT "UC0 , 5 . 3 , 99 , 3 , -8 , -6 , 0 , 3 , 8 ; " 

IF Syrn=7 THEN PRINT "UC0 , -5 . 3 , 99, -3 , 8 , 6 , 0 , -3 , -8 ; " 

PRINT "SM" 

IF Syrn<4 THEN PRINT "PR 2,0" 

PRINT "PR 2,-1.0;SR 1 . 0 , 1 . 8 : LB" : D_f i le$ ; " " 

YaI=Yal-5 

BEEP 

INPUT "WANT TO ADD ANOTHER STRING ( 1 = Y , 0=N )? " , las 
IF Ias=l THEN 
BEEP 

INPUT "ENTER THE STRING" , Label* 

PRINT "PR 2,0:SR 1 . 0 , 1 . 8 : LB " ; Label* ; " " 

GOTO 7990 
END IF 
END IF 
BEEP 

INPUT "WANT TO COMBINE ANOTHER FILE? ( 1=Y , 0=N ) " , Iconb 
ASSIGN ©File TO * 

IF lcomb<>0 THEN 7645 

Ils=l 

BEEP 

INPUT "WANT TO PLOT A LEAST-SQUARES LINE < 1=DEF=YES , 0=NO ) " , 1 1 s 
IF Ils=l THEN 
BEEP 

INPUT "SELECT CURVE TYPE ( 0=SOLID , 1=DASHED ) " >, 1 1 1 

Ilt=Ilt+l 

PRINT "SM" 

IF lexp=0 THEN 

Bb=( Nrun*Sxy-Sy*Sx )/(Nrun*Sx2-Sx'‘2 ) 

ELSE 

Bb=.75 
END IF 

Aa=( Sy-Bb*Sx )/Nrun 
PRINTER IS 1 

PRINT USING "10X,""a = " " , MZ . 3DE " : Aa 
PRINT USING "10X,""n = " " , MZ . 3DE" ; Bb 
PRINTER IS 705 
ln=0 

FOR Xa=Xnin TO Xnax STEP (X(nax-Xmn> 

Ya=Aa+Xa*Bb 
Y = ( Ya-Ynin )*Sfy 
X=<Xa-Xnin )*Sfx 
IF Y<0 THEN Y=0 
IF Y>100 THEN GOTO 8220 
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8180 

8185 

8190 

8195 

8200 

8205 

8210 

8215 

8220 

8225 

8230 

8255 

8240 

8245 

8250 

8255 

8260 

8270 

8280 

8290 

8300 

8320 

8330 

8340 

8350 

8370 

8380 

8390 

8400 

8420 

8430 

8440 

8450 

8460 

8470 

8480 

8490 



IF 111=1 THEN 

PRINT "PA" ,X,Y, "PD" 

ELSE 

In=In+l 
Ir=In MOD 111 

IF Ir=l THEN PRINT "PA" , X , Y , "PD" 
IF lr=0 THEN PRINT "PA" ,X , Y , "PU" 
END IF 
NEXT Xe 
PRINT "PU" 

END IF 
lcomb=0 
GOTO 7620 
END IF 

PRINT "PU PA 0,0" 

PRINT "PU PA 0,0 SP0" 

SUBEND 

DEF FNFncaK I ,T,Pc ) 

IF 1=0 THEN 

Mdot0=.00206*Pc-.00414 
Mdot24=.00225*Pc-. 00329 
END IF 
IF 1=1 THEN 

Mdol0=.00204*Pc-.00615 
Mdol24=.00216*Pc-. 00138 
END IF 
IF 1=2 THEN 

Mdot0=.00208*Pc-.0036 
Mdot24=.00223*Pc-. 00301 
END IF 
IF 1=3 THEN 

Mdot0= . 002 8* Pc- .00122 
Mdol24= . 002 9* Pc + . 0048 
END IF 

Mdt = (T*(Mdot24-Mdot0)/24 )+Mdot0 

RETURN Mdl 

FNEND 
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